Well placement optimization subject to realistic field development constraints

This work considers the well placement problem in reservoir management and field development optimization. In particular, it emphasizes embedding realistic and practical constraints into a mathematical optimization formulation. Such constraints are a prerequisite for the wider use of mathematical optimization techniques in well placement problems, since constraints are a way to incorporate reservoir engineering knowledge into the problem formulation. There are important design limitations that are used by the field development team when treating the well placement problem, and these limitations need to be articulated and eventually formalized within the problem before conducting the search for optimal well placements. In addition, these design limitations may be explicit or implicit. In this work, various design limitations pertaining to well locations have been developed in close collaboration with a field operator on the Norwegian Continental Shelf. Moreover, this work focuses on developing constraint-handling capability to enforce these various considerations during optimization. In particular, the Particle Swarm Optimization (PSO) algorithm is applied to optimize for the well locations, and various practical well placement constraints are incorporated into the PSO algorithm using two different constraint-handling techniques: a decoder procedure and the penalty method. The decoder procedure maps the feasible search space onto a cube and has the advantage of not requiring parameter tuning. The penalty method converts the constrained optimization problem into an unconstrained one by introducing an additional term, which is called a penalty function, to the objective function. In contrast to the penalty method, only feasible solutions are evaluated in the decoder method. Through numerical simulations, a comparison between the penalty method and the decoder technique is performed for two cases. We show that the decoder technique can easily be implemented for the well placement problem, and furthermore, that it performs better than the penalty method in most of the cases.     

Recurrent mesoproterozoic continental magmatism in South-Central Norway

We report U–Pb dates and Lu–Hf isotope data, obtained by LAM-ICPMS, for zircons from metamorphic rocks of the Setesdalen valley, situated in the Telemark block south of the classic Telemark region of southern Norway. The samples include infracrustal rocks from the metamorphic basement, metaigneous rocks and metasediments from the Byglandsfjorden supracrustal cover sequence, and metaigneous rocks which intruded the whole succession. The main crustal evolution took place from 1,550–1,020 Ma, beginning with the emplacement of juvenile tonalitic melts; the contribution of older crustal material increased with time. Around 1,320 Ma, further addition of juvenile material occurred, involving both mafic and felsic melts, metamorphism and deformation. Acid magmas with high FeO*/MgO were intruded at 1,215 Ma, coinciding with underplating elsewhere in South Norway. The period starting at 1,215 Ma is represented by supracrustal rocks, principally metarhyolites with minor mafic material and immature sediments of the Byglandsfjorden Group. The crust generation processes ended with the intrusion of diorites and granodiorites at 1,030 Ma, late in the Sveconorwegian orogeny. Regional processes of metamorphism and deformation (around 1,290 and 1,000 Ma) can be related to the assembly of Rodinia. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A phenomenological model for disparitions brusques followed by flarelike chromospheric brightenings

Here I present a simple gravitational model for the flarelike brightenings of the chromosphere that follow most disparitions brusques (disappearing filaments). I assume the ascending prominence material is lifted out of the initially stable magnetic dips that characterize quiescent prominences and falls along the arched field lines into the chromosphere where the kinetic energy of fall is dissipated in the bright areas. The examination of prominence and chromospheric characteristics leads naturally to many predictions and relations during and after prominence eruptions. In general the predictions are specific, but the observations of necessary detail and quality are nonexistent; however, the predictions appear to agree with the data that are available. The model appears to explain all non-active-region brightenings of the chromosphere that follow disparitions brusques and an unknown fraction of active-region flares. The conclusion is that two-ribbon flares are due to the disparitions brusques chromospheric flarelike brightening mechanism. In this paper it will become clear that many specific observations in and out of active regions will be necessary to test the predictions of the model given here.     

Broad search for unstable resonant orbits in the planar circular restricted three-body problem

Unstable resonant orbits in the circular restricted three-body problem have increasingly been used for trajectory design using optimization and invariant manifold techniques. In this study, several methods for computing these unstable resonant orbits are explored including grid searches, flyby maps, and continuation. Families of orbits are computed focusing on orbits with multiple loops near the secondary in the Jupiter–Europa system, and their characteristics are explored. Different parameters such as period and stability are examined for each set of resonant orbits, and the continuation of several specific orbits is explored in more detail.     

Non-Linear Theory and Power-Law Models for Information Integration and Mineral Resources Quantitative Assessments

Singular physical or chemical processes may result in anomalous amounts of energy release or mass accumulation that, generally, are confined to narrow intervals in space or time. Singularity is a property of different types of non-linear natural processes including cloud formation, rainfall, hurricanes, flooding, landslides, earthquakes, wildfires, and mineralization. The end products of these non-linear processes can be modeled as fractals or multifractals. Hydrothermal processes in the Earth’s crust can result in ore deposits characterized by high concentrations of metals with fractal or multifractal properties. Here we show that the non-linear properties of the end products of singular mineralization processes can be applied for prediction of undiscovered mineral deposits and for quantitative mineral resource assessment, whether for mineral exploration or for regional, national and global planning for mineral resource utilization. In addition to the general theory and framework for the non-linear mineral resources assessment, this paper focuses on several power-law models proposed for characterizing non-linear properties of mineralization and for geoinformation extraction and integration. The theories, methods, and computer system discussed in this paper were validated using a case study dealing with hydrothermal Au mineral potential in southern Nova Scotia, Canada.     

A simplified method for analyzing the fundamental frequency of monopile supported offshore wind turbine system design

Preliminary design of offshore wind turbines requires high precision simplified methods for the analysis of the system fundamental frequency. Based on the Rayleigh method and Lagrange's Equation, this study establishes a simple formula for the analysis of system fundamental frequency in the preliminary design of an offshore wind turbine with a monopile foundation. This method takes into consideration the variation of cross-section geometry of the wind turbine tower along its length, with the inertia moment and distributed mass both changing with diameter. Also the rotational flexibility of the monopile foundation is mainly considered. The rigid pile and elastic middle long pile are calculated separately. The method is validated against both FEM analysis cases and field measurements, showing good agreement. The method is then used in a parametric study, showing that the tower length Lt, tower base diameter d0, tower wall thickness δt, pile diameter db and pile length Lb are the major factors influencing the fundamental frequency of the offshore wind turbine system. In the design of offshore wind turbine systems, these five parameters should be adjusted comprehensively. The seabed soil condition also needs to be carefully considered for soft clay and loose sand.     

A Note on Reviving the Goddard Satellite-Based Surface Turbulent Fluxes （GSSTF） Dataset

Accurate sea surface flux measurements are crucial for understanding the global water and energy cycles. The oceanic evaporation, which is a major component of the global oceanic fresh water flux, is useful for predicting oceanic circulation and transport. The global Goddard Satellite-based Surface Turbulent Fluxes Version-2 (GSSTF2; July 1987--December 2000) dateset that was officially released in 2001 has been widely used by scientific community for global energy and water cycle research, and regional and short period data analyses. We have recently been funded by NASA to resume processing the GSSTF dataset with an objective of continually producing a uniform dataset of sea surface turbulent fluxes, derived from remote sensing data. The dataset is to be reprocessed and brought up-to-date (GSSTF2b) using improved input datasets such as a recently upgraded NCEP/DOE sea surface temperature reanalysis, and an upgraded surface wind and microwave brightness temperature V6 dataset (Version 6) from the Special Sensor Microwave Imager (SSM/I) produced by Remote Sensing Systems (RSS). A second new product (GSSTF3) is further proposed with a finer temporal (12-h) and spatial (0.25^ox0.25^o) resolution. GSSTF2b (July 1987--December 2008) and GSSTF3 (July 1999--December 2009) will be released for the research community to use by late 2009 and early 2011, respectively.     

Bulk Transfer Relations for the Roughness Sublayer

In the roughness sublayer (RSL), Monin–Obukhov surface layer similarity theory fails. This is problematic for atmospheric modelling applications over domains that include rough terrain such as forests or cities, since in these situations numerical models often have the lowest model level located within the RSL. Based on empirical RSL profile functions for momentum and scalar quantities, and scaling the height with the RSL height z _*, we derive a simple bulk transfer relation that accounts for RSL effects. To verify the validity of our approach, these relations are employed together with wind speed and temperature profiles measured over boreal forest during the BOREAS experimental campaign to estimate momentum and heat fluxes. It is demonstrated that, when compared with observed flux values, the inclusion of RSL effects in the transfer relations yields a considerable improvement in the estimated fluxes.     

Errors and correction of precipitation measurements in China

In order to discover the range of various errors in Chinese precipitation measurements and seek a correction method, 30 precipitation evaluation stations were set up countrywide before 1993. All the stations are reference stations in China. To seek a correction method for wind-induced error, a precipitation correction instrument called the ＂horizontal precipitation gauge＂ was devised beforehand. Field intercomparison observations regarding 29,000 precipitation events have been conducted using one pit gauge, two elevated operational gauges and one horizontal gauge at the above 30 stations. The range of precipitation measurement errors in China is obtained by analysis of intercomparison measurement results. The distribution of random errors and systematic errors in precipitation measurements are studied in this paper. A correction method, especially for wind-induced errors, is developed. The results prove that a correlation of power function exists between the precipitation amount caught by the horizontal gauge and the absolute difference of observations implemented by the operational gauge and pit gauge. The correlation coefficient is 0.99. For operational observations, precipitation correction can be carried out only by parallel observation with a horizontal precipitation gauge. The precipitation accuracy after correction approaches that of the pit gauge. The correction method developed is simple and feasible.     

Surface and groundwater quality characterization of Deoria District,Ganga Plain,India

A water quality investigation was carried out in the Deoria district, Ganga plain, to assess the suitability of surface and groundwaters for domestic, agricultural, and industrial purposes. As much as 50 representative samples from river and groundwater were collected from various stations to monitor the water chemistry of various ions, comprising Ca²⁺, Mg²⁺, Na⁺, K⁺, HCO₃ ⁻, SO₄ ²⁻, NO₃ ⁻, Cl⁻, F⁻, and trace metals, such as Fe, Cu, Mn, Zn, Cd, and Pb. The results showed that electrical conductance (EC), total dissolved solids (TDS), HCO₃ ⁻, Mg²⁺, Na⁺, and total hardness (TH) are above the maximum desirable limit, and apart from Fe and Mn all other trace metals are within the maximum permissible limit for drinking water. The calculated values for sodium absorption ratio (SAR), salinity, residual sodium carbonate (RSC), and permeability index (PI) indicate well to permissible use of water for irrigation. High values of Na%, RSC, and Mg-hazard (MH) at some stations restrict its use for agricultural purpose. Anthropogenic activities affect the spatial variation of water quality. Economic and social developments of the study area is closely associated with the characteristics of the hydrological network.