Successful application of multiscale methods in a real reservoir simulator environment     

K.-A. Lie  O. Møyner  J. R. Natvig  A. Kozlova  K. Bratvedt  S. Watanabe  Z. Li 《Computational Geosciences》2017,21(5-6):981-998

For the past 10 years or so, a number of so-called multiscale methods have been developed as an alternative approach to upscaling and to accelerate reservoir simulation. The key idea of all these methods is to construct a set of prolongation operators that map between unknowns associated with cells in a fine grid holding the petrophysical properties of the geological reservoir model and unknowns on a coarser grid used for dynamic simulation. The prolongation operators are computed numerically by solving localized flow problems, much in the same way as for flow-based upscaling methods, and can be used to construct a reduced coarse-scale system of flow equations that describe the macro-scale displacement driven by global forces. Unlike effective parameters, the multiscale basis functions have subscale resolution, which ensures that fine-scale heterogeneity is correctly accounted for in a systematic manner. Among all multiscale formulations discussed in the literature, the multiscale restriction-smoothed basis (MsRSB) method has proved to be particularly promising. This method has been implemented in a commercially available simulator and has three main advantages. First, the input grid and its coarse partition can have general polyhedral geometry and unstructured topology. Secondly, MsRSB is accurate and robust when used as an approximate solver and converges relatively fast when used as an iterative fine-scale solver. Finally, the method is formulated on top of a cell-centered, conservative, finite-volume method and is applicable to any flow model for which one can isolate a pressure equation. We discuss numerical challenges posed by contemporary geomodels and report a number of validation cases showing that the MsRSB method is an efficient, robust, and versatile method for simulating complex models of real reservoirs.  相似文献   

Mapping vertical profile of discontinuous permafrost with ground penetrating radar at Nalaikh depression,Mongolia   总被引：3，自引：0，他引：3  

Tonghua Wu  Qinxue Wang  Masataka Watanabe  Ji Chen  Dorjgotov Battogtokh 《Environmental Geology》2009,56(8):1577-1583

Ground-penetrating radar (GPR) has become an important geophysical tool which can provide a wealth of interpretive information about the vertical profile of discontinuous permafrost. A GPR investigation was conducted in October 2006 at the Nalaikh site at the southern boundary of the Siberian discontinuous permafrost region in Mongolia. GPR data were collected along four 100-m-long profiles to identify the location of the permafrost body, which included an in situ drilling borehole and analysis of temperature observations and soil water content measurements from boreholes. The GPR interpretation results indicated that the thickness of discontinuous permafrost at the study site was only 1.9–3.0 m and the permafrost is vulnerable to climate change. The soil temperature and soil water content data demonstrate the precision of GPR image interpretation. This case demonstrated that GPR is well suited for mapping the internal structure of discontinuous permafrost with relatively low soil water content.  相似文献   

Spatial distribution of sister species of vesicomyid bivalves Calyptogena okutanii and Calyptogena soyoae along an environmental gradient in chemosynthetic biological communities in Japan     

Hiromi Watanabe  Eriko Seo  Yoshimi Takahashi  Takao Yoshida  Shigeaki Kojima  Katsunori Fujikura  Hiroshi Miyake 《Journal of Oceanography》2013,69(1):129-134

Vesicomyid bivalves have a substantial biomass in deep-sea chemosynthetic biological communities in the Pacific. Using a novel multiplex-PCR (mPCR) method to identify the co-occurring vesicomyids in Sagami Bay, we analyzed the distribution of Calyptogena okutanii and Calyptogena soyoae along environmental gradients. All the known distributions of C. okutanii indicated the different preferences in salinity and temperature to those of C. soyoae, and in Sagami Bay, depth seemed to be an important environmental factor, too. Although the concentration of hydrogen sulfide in sediment was not examined, our results showed that the distributions of these two Calyptogena clams were affected by salinity and temperature.  相似文献   

Holocene paleolimnological changes in Lake Skallen Oike in the Syowa Station area of Antarctica inferred from organic components in a sediment core (Sk4C-02)     

Genki I. Matsumoto  Yukinori Tani  Koji Seto  Tomoko Tazawa  Masumi Yamamuro  Takahiro Watanabe  Toshio Nakamura  Tetsuo Takemura  Satoshi Imura  Hiroshi Kanda 《Journal of Paleolimnology》2010,44(2):677-693

Antarctic climate changes influence environmental changes at both regional and local scales. Here we report Holocene paleolimnological changes in lake sediment core Sk4C-02 (length 378.0 cm) from Lake Skallen Oike in the Soya Kaigan region of East Antarctica inferred from analyses of sedimentary facies, a range of organic components, isotope ratios of organic carbon and nitrogen, and carbon-14 dating by Tandetron accelerator mass spectrometry. The sediment core was composed of clayish mud (378.0–152.5 cm) overlain by organic sediments (152.5 cm-surface). The age of the surface and the core bottom were 150 (AD1950-1640) and ca. 7,030 ± 73 calibrated years before present (cal BP), respectively, and the mean sedimentation rate was estimated to be 0.55 mm/year. Multi-proxy analyses revealed that the principal environmental change in the core is a transition from marine to lacustrine environments which occurred at a depth of 152.5 cm (ca. 3,590 cal BP). This was caused by relative sea level change brought about by ongoing retreat of glaciers during the mid-Holocene warming of Antarctica, and ongoing isostatic uplift which outpaced changes in global (eustatic) sea level. The mean isostatic uplift rate was calculated to be 2.8 mm/year. The coastal marine period (378.0–152.5 cm, ca. 7,030–3,590 cal BP) was characterized by low biological production with the predominance of diatoms. During the transition period from marine to freshwater conditions (152.5-approximately 135 cm, ca. 3,590–3,290 cal BP) the lake was stratified with marine water overlain by freshwater, with a chemocline and an anoxic (sulfidic) layer in the bottom of the photic zone. Green sulfur bacteria and Cryptophyta were the major photosynthetic organisms. The Cryptophyta appeared to be tolerant of the moderate salinity and stratified water conditions. The lacustrine period (approximately 135 cm-surface, ca. 3,290 cal BP-present) was characterized by high biological production by green algae (e.g. Comarium clepsydra and Oedegonium spp.) with some contributions from cyanobacteria and diatoms. Biological production during this period was 8.7 times higher than during the coastal marine period.  相似文献   

The effect of a hydrous phase on P‐wave velocity anisotropy within a detachment shear zone in the slow‐spreading oceanic crust: A case study from the Godzilla Megamullion,Philippine Sea          下载免费PDF全文

Katsuyoshi Michibayashi  Tohru Watanabe  Yumiko Harigane  Yasuhiko Ohara 《Island Arc》2016,25(3):209-219

We studied the contributions of plagioclase, clinopyroxene, and amphibole to the P‐wave velocity properties of gabbroic mylonites of the Godzilla Megamullion (site KH07‐02‐D18) in the Parece Vela Rift of the central Parece Vela Basin, Philippine Sea, based on their crystal‐preferred orientations (CPOs), mineral modes, and elastic constants and densities of single crystals. The gabbroic mylonites have been classified into three types based on their microstructures and temperature conditions: HT1, HT2 and medium‐temperature (MT) mylonites. The P‐wave velocity properties of the HT1 mylonite are dominantly influenced by plagioclase CPOs. Secondary amphibole occurred after deformation in the HT1 mylonite, so that its effect on P‐wave velocity anisotropy is minimal due to weak CPOs. Although the HT2 mylonite developed deformation microstructures in the three minerals, the P‐wave velocity properties of the HT2 mylonite are essentially isotropic, resulting from the destructive interference of different P‐wave velocity anisotropy patterns produced by the distinct CPOs of the three constituent minerals (i.e., plagioclase, clinopyroxene, and amphibole). The P‐wave velocity properties of the MT mylonite are influenced mainly by amphibole CPOs, whereas the effect of plagioclase CPOs on P‐wave velocity anisotropy becomes very small with a decrease in the intensity of plagioclase CPOs. As a result, the gabbroic mylonites tend to have weak P‐wave velocity anisotropy in seismic velocity, although their constituent minerals show distinct CPOs. Such weakness in the whole‐rock P‐wave velocity anisotropy could result from the destructive contributions of the different mineral CPOs with respect to the structural framework (foliation and lineation). These results show that amphibole has a high potential for P‐wave velocity anisotropy by aligning both crystallographically and dimensionally during deformation in the hydrous oceanic crust. The results also suggest that the effect of a hydrous phase on P‐wave velocity anisotropy within the detachment shear zone in a slow‐spreading oceanic crust varies depending on the degree of deformation and on the timing of hydrothermal activity.  相似文献   

Distribution of local Cs anomalies on the seafloor near the Fukushima Dai-ichi Nuclear Power Plant     

Blair Thornton  Seiki Ohnishi  Tamaki Ura  Naoteru Odano  Shun Sasaki  Tsuneo Fujita  Tomowo Watanabe  Kaoru Nakata  Tsuneo Ono  Daisuke Ambe 《Marine pollution bulletin》2013

An estimated 3.5 ± 0.7 × 10¹⁵ Bq of ¹³⁷Cs is thought to have been discharged into the ocean following the melt down at Fukushima Dai-ichi Nuclear Power Plant (F1NPP). While efforts have been made to monitor seafloor radiation levels, the sampling techniques used cannot capture the continuous distribution of radionuclides. In this work, we apply in situ measurement techniques using a towed gamma ray spectrometer to map the continuous distribution of ¹³⁷Cs on the seafloor within 20 km of the F1NPP. The results reveal the existence of local ¹³⁷Cs anomalies, with levels of ¹³⁷Cs an order of magnitude higher than the surrounding seafloors. The sizes of the anomalies mapped in this work range from a few meters to a few hundreds of meters in length, and it is demonstrated that the distribution of these anomalies is strongly influenced by meter scale features of the terrain.  相似文献   

Assessing image processing techniques for geological mapping: a case study in Eljufra,Libya     

N. M. Saadi  K. Watanabe 《国际地球制图》2013,28(3):241-253

Various image processing techniques were experimented with in this study to evaluate their efficiency for geological mapping in the Eljufra area of northwest Libya. Remote sensing data including multi-spectral optical Landsat Enhanced Thematic Mapper (ETM+), Synthetic Aperture Radar (ERS-2 SAR) and Digital Elevation Models (DEMs) extracted from the Shuttle Radar Topography Mission (SRTM) data were used to trace different lithological units as well as extracting geological lineaments in the study area. The study area is located in an arid environment mostly devoid of any vegetation. Most lithological and structural units are distinguishable based on their topographic form and spectral properties. Fusion of ETM+ and ERS-2 images was experimented with to further identify lithological units. Shaded relief techniques were implemented to enhance terrain perspective views and to extract geological lineaments. The results discriminated different rock units and modified formation boundaries and revealed new geological lineaments. Nine rock units were identified and plotted in the new geological map defined by the new boundaries. The dominant lineaments tend to run in the NNW-SSE and NNE-SSW directions. Analysis and interpretation of the lineaments provided information about the tectonic evolution of the study area.  相似文献   

Study on the selection of unsaturated flow model for the different types of soil and soft rock   总被引：1，自引：1，他引：0  

Sajeewani Rajika Amarasinghe  Kunio Watanabe  Koji Ishiyama 《Environmental Earth Sciences》2011,64(7):1795-1805

Precise estimation of unsaturated hydraulic properties of porous media is indispensable in various study areas, such as analyzing the moisture flow, the drying process occurring from the surface, and the pollutant migration beneath the ground surface. Although many empirical/theoretical models describing the unsaturated hydraulic properties have been proposed by several previous researchers, the best model for the different types of soil/rock may not be identical. Thus, the model selection process and the estimation technique of the parameters included in the models should be developed. In the present study, the inverse technique based on the transient evaporation change was investigated to select the model and estimate the model parameters. The experimental work was based on a relatively low permeable soft rock and a relatively high permeable sandy soil (Toyoura standard sand). Experimental equipment was developed to precisely measure the evaporation rate for the high permeable sandy soil. The Genetic Algorithm (GA) was adopted in the inverse technique as an optimization tool. In order to simplify the problem, only the drying process from the saturated condition was considered. It was established that the information concerning the transient evaporation change could be used for the model selection and parameter estimation. Further, the saturation distribution could be used for the selection of the models. The present study provides important information for the development of the model selection process.  相似文献   

Evaluation of image processing methods for geological interpretation in the different environments in Libya   总被引：1，自引：0，他引：1  

Nureddin M. Saadi  Koichiro Watanabe  Hideki Mizunaga 《Arabian Journal of Geosciences》2011,4(3-4):635-643

This study evaluates the use of image processing techniques and methodologies of digital integration of multi-disciplinary geoscientific data to reduce the ambiguity in geological interpretations in different geological environments in Libya. To realize this objective, three areas were selected for this study; Eljufra, an arid environment; the mountainous environment Tarhunah; and the Jifara Plain, a coastal plain environment. Two types of remote sensing data were used in this study: Landsat Enhanced Thematic Mapper Plus (ETM+) and European Remote Sensing Satellite. The digital elevation model extracted from Shuttle Radar Topography Mission and digital topographic maps scale 1:50,000 were used for remote sensing interpretation. GIS and remote sensing-based methods were used to process and integrate all raster and vector layers data. This study offered dramatic benefits for geological interpretations and provided new insights into the efficiency of image processing methods in different geological environments.  相似文献   

Phytoplankton productivity in the western subarctic gyre of the North Pacific in early summer 2006     

Tetsuichi Fujiki  Kazuhiko Matsumoto  Shuichi Watanabe  Takuji Hosaka  Toshiro Saino 《Journal of Oceanography》2011,67(3):295-303

We deployed a profiling buoy system incorporating a fast repetition rate fluorometer in the western subarctic Pacific and carried out time-series observations of phytoplankton productivity from 9 June to 15 July 2006. The chlorophyll a (Chl a) biomass integrated over the euphotic layer was as high as 45–50 mg Chl a m⁻² in the middle of June and remained in the 30–40 mg Chl a m⁻² range during the rest of observation period; day-to-day variation in Chl a biomass was relatively small. The daily net primary productivity integrated over the euphotic layer ranged from 144 to 919 mg C m⁻² day⁻¹ and varied greatly, depending more on insolation rather than Chl a biomass. In addition, we found that part of primary production was exported to a 150-m depth within 2 days, indicating that the variations in primary productivity quickly influenced the organic carbon flux from the upper ocean. Our results suggest that the short-term variability in primary productivity is one of the key factors controlling the carbon cycle in the surface ocean in the western subarctic Pacific.  相似文献