Application of the EnKF and Localization to Automatic History Matching of Facies Distribution and Production Data

The performance of the ensemble Kalman filter (EnKF) for continuous updating of facies location and boundaries in a reservoir model based on production and facies data for a 3D synthetic problem is presented. The occurrence of the different facies types is treated as a random process and the initial distribution was obtained by truncating a bi-Gaussian random field. Because facies data are highly non-Gaussian, re-parameterization was necessary in order to use the EnKF algorithm for data assimilation; two Gaussian random fields are updated in lieu of the static facies parameters. The problem of history matching applied to facies is difficult due to (1) constraints to facies observations at wells are occasionally violated when productions data are assimilated; (2) excessive reduction of variance seems to be a bigger problem with facies than with Gaussian random permeability and porosity fields; and (3) the relationship between facies variables and data is so highly non-linear that the final facies field does not always honor early production data well. Consequently three issues are investigated in this work. Is it possible to iteratively enforce facies constraints when updates due to production data have caused them to be violated? Can localization of adjustments be used for facies to prevent collapse of the variance during the data-assimilation period? Is a forecast from the final state better than a forecast from time zero using the final parameter fields?To investigate these issues, a 3D reservoir simulation model is coupled with the EnKF technique for data assimilation. One approach to enforcing the facies constraint is continuous iteration on all available data, which may lead to inconsistent model states, incorrect weighting of the production data and incorrect adjustment of the state vector. A sequential EnKF where the dynamic and static data are assimilated sequentially is presented and this approach seems to have solved the highlighted problems above. When the ensemble size is small compared to the number of independent data, the localized adjustment of the state vector is a very important technique that may be used to mitigate loss of rank in the ensemble. Implementing a distance-based localization of the facies adjustment appears to mitigate the problem of variance deficiency in the ensembles by ensuring that sufficient variability in the ensemble is maintained throughout the data assimilation period. Finally, when data are assimilated without localization, the prediction results appear to be independent of the starting point. When localization is applied, it is better to predict from the start using the final parameter field rather than continue from the final state.     

Markov chain Monte Carlo methods for conditioning a permeability field to pressure data

Generating one realization of a random permeability field that is consistent with observed pressure data and a known variogram model is not a difficult problem. If, however, one wants to investigate the uncertainty of reservior behavior, one must generate a large number of realizations and ensure that the distribution of realizations properly reflects the uncertainty in reservoir properties. The most widely used method for conditioning permeability fields to production data has been the method of simulated annealing, in which practitioners attempt to minimize the difference between the ’ ’true and simulated production data, and “true” and simulated variograms. Unfortunately, the meaning of the resulting realization is not clear and the method can be extremely slow. In this paper, we present an alternative approach to generating realizations that are conditional to pressure data, focusing on the distribution of realizations and on the efficiency of the method. Under certain conditions that can be verified easily, the Markov chain Monte Carlo method is known to produce states whose frequencies of appearance correspond to a given probability distribution, so we use this method to generate the realizations. To make the method more efficient, we perturb the states in such a way that the variogram is satisfied automatically and the pressure data are approximately matched at every step. These perturbations make use of sensitivity coefficients calculated from the reservoir simulator.     

额尔古纳超大型火山热液型铀成矿带地质特征及找矿前景

本文介绍了在额尔古纳超大型火山热液型铀成矿带中发现的超大型火山热液脉型铀矿末有５个，大型火山热液脉型铀矿床有６个，总共有铀金属储量３３万ｔ，在世界热液脉型铀矿资源量中占有重要地位。经研究认为，该成矿带内铀矿就成矿作用在时间上具有突发性特征，成矿时代集中在早白垩世，受统一中生代大陆裂陷－火山岩浆活动控制；成矿介质中深源的火山岩水，大气水和变质水三者相混合的具中低温的火山热液，成矿物质主要来自中生代深     

On the generation of oceanic surface waves by underwater volcanic explosions

The oceanic surface waves generated by submarine volcanism are modeled as the response of a compressible ocean to an underwater point explosion at an arbitrary depth over an elastic half-space. We generalize the work of Kajiura and Miyoshi by allowing the elastic sea floor to influence both the sound and surface gravity waves. The classic work of Press and Ewing is extended to allow for the movement of the oceanic free surface. From numerical inversion techniques, we find that most of the energy is radiated away as surface gravity waves. A mixture of direct and reflected P- (compressional) and S-type (distortional) head waves may produce high-frequency, large-amplitude oscillations in the sea surface.     

Optimal Routefinding Across Landscapes Featuring High‐cost Linear Obstacles

Existing GIS‐based least‐cost routefinding approaches can in some cases be confounded by linear features in the landscape whose crossing costs are high relative to other traversing costs found throughout the landscape. Unfortunately, such high‐cost linear features are not uncommon; they frequently occur in the form of hydrologic features whose crossing costs are high relative to costs of traversing the surrounding dry land. This study (1) enumerates the situations where existing approaches can fail, (2) proposes a simple method for overcoming the limitations of the existing approaches, and (3) conducts an experiment to assess the impact of the weaknesses of existing approaches and the magnitude of the differences between the results produced by existing and revised approaches. Our results show that in mountainous terrain, linear hydrologic features with high crossing costs can have a profound impact on least‐cost routes, and the choice of solution method has equally profound impacts on the optimal route produced by the analysis.     

Mapping local regression for spatial object-pairs

Local regression methods such as geographically weighted regression (GWR) can provide specific information about individual locations (or places) in spatial analysis that is useful for mapping nonstationary covariate relationships. However, the distance-based weighting schemes used in GWR are only adaptable for spatial objects that are point or area features. In particular, spatial object-pairs pose a challenge for local analysis because they have a linear dimensionality rather than a point dimensionality. This paper proposes using an alternative local regression model – quantile regression (QR) – for investigating the stationarity of regression parameters with respect to these linear features as well as facilitating the visualization of the results. An empirical example of a gravity model analysis of trade patterns within Europe is used to illustrate the utility of the proposed method.     

Seasonal and regional patterns in egg production of Calanus glacialis/marshallae in the Chukchi and Beaufort Seas during spring and summer, 2002

Understanding the physical and biogeochemical processes that control the exchange of biogenic carbon within and between the arctic shelves, slopes, and deep basins is a key objective of the Western Arctic Shelf-Basin Interaction program (SBI). Here, egg production (EP) of the dominant copepod Calanus glacialis/marshallae was used as an indicator of food limitation for the mesozooplankton community in the Chukchi and Beaufort Seas in spring and summer, 2002. Both C. glacialis and C. marshallae may occur in this region but the two cannot easily be differentiated visually. Four oceanographic regions were objectively identified that roughly corresponded to the different pathways in circulation of nutrient-rich Pacific water. A ‘transition’ region characterized by ‘older’ Pacific water was located at the shelfbreak and separated the nutrient-rich shelf water and the low-nutrient waters of the deep basin. The observed spatial pattern in EP in C. glacialis/marshallae in spring and summer resulted both from the different water mass environments and from the reproductive cycle of the species. EP was greater on the shelf than in the basin, corresponding to differences in body size and nitrogen condition factor (NCF) in females, while the egg viability was generally high throughout the study area. EP showed no relationship with low-chlorophyll a biomass under heavy ice-cover in spring, while a significant relationship was observed in the more open water in summer. Adult female carbon condition factor (CCF) was much higher in summer, reflecting the accumulation of lipids during the growth season. Small animals with a markedly greater NCF dominated on the shelf. The shelfbreak region contained a mixture of females from the shelf and the basin with intermediate sizes, conditions, and EP rates. The occurrence of water typical of the ‘transition’ shelfbreak region and elevated EP in C. glacialis/marshallae offshore on the Barrow Canyon and East Barrow sections indicated offshore transport of productive shelf water and the associated plankton community. The input of nutrient-rich Pacific water and accompanying elevated production to the northern Chukchi Sea and the Chukchi-Beaufort shelfbreak region may contribute to the reproductive success of C. glacialis/marshallae in this region.     

The role of oxyradicals in intracellular proteolysis and toxicity in mussels

Studies were performed on the common mussel, M. edulis L., to determine whether copper (Cu) exposure can affect the extent to which digestive cell proteins are oxidised and whether such oxidative damage is mediated by free radicals. Three age groups of mussels were exposed for 6 ^-days to environmentally realistic concentrations of Cu and then digestive gland homogenates were examined for evidence of protein carbonyl formation. Significant increases in carbonyls relative to untreated control mussels were seen for the youngest (2–4 year-old) and oldest (≥ 10 year-old) mussels only after exposure for 6 days, followed by recovery from exposure for a further 6 days. Untreated mussels also showed an age-related difference in protein oxidation, with a significantly lower concentration in the youngest animals (2–4 year olds). Copper did not affect the levels of modified tryptophan or tyrosine residues or the extent of total lipid peroxidation in digestive gland homogenate. Significant depletion of total vitamin E (a-tocopherol) was seen only in young and medium-aged mussels following exposure for 6 days. The levels of protein carbonyl groups were increased in digestive cell cytosol and lighter lysosomes but not in heavier lysosomes or digestive gland microsomes following 5 days exposure to Cu. Dihydrohodamine-123 was converted to fluorescent rhodamine-123 following sequestration into digestive cell lysosomes. The results suggest a link between the lysosomal sequestration of copper, a concomitant increase in the production of oxyradicals and the potential for intracellular oxidative damage, as well as an increased capacity for oxidative damage in older animals.