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.     

Estimating urban above ground biomass with multi-scale LiDAR

Background

Urban trees have long been valued for providing ecosystem services (mitigation of the “heat island” effect, suppression of air pollution, etc.); more recently the potential of urban forests to store significant above ground biomass (AGB) has also be recognised. However, urban areas pose particular challenges when assessing AGB due to plasticity of tree form, high species diversity as well as heterogeneous and complex land cover. Remote sensing, in particular light detection and ranging (LiDAR), provide a unique opportunity to assess urban AGB by directly measuring tree structure. In this study, terrestrial LiDAR measurements were used to derive new allometry for the London Borough of Camden, that incorporates the wide range of tree structures typical of an urban setting. Using a wall-to-wall airborne LiDAR dataset, individual trees were then identified across the Borough with a new individual tree detection (ITD) method. The new allometry was subsequently applied to the identified trees, generating a Borough-wide estimate of AGB.

Results

Camden has an estimated median AGB density of 51.6 Mg ha^–1 where maximum AGB density is found in pockets of woodland; terrestrial LiDAR-derived AGB estimates suggest these areas are comparable to temperate and tropical forest. Multiple linear regression of terrestrial LiDAR-derived maximum height and projected crown area explained 93% of variance in tree volume, highlighting the utility of these metrics to characterise diverse tree structure. Locally derived allometry provided accurate estimates of tree volume whereas a Borough-wide allometry tended to overestimate AGB in woodland areas. The new ITD method successfully identified individual trees; however, AGB was underestimated by ≤?25% when compared to terrestrial LiDAR, owing to the inability of ITD to resolve crown overlap. A Monte Carlo uncertainty analysis identified assigning wood density values as the largest source of uncertainty when estimating AGB.

Conclusion

Over the coming century global populations are predicted to become increasingly urbanised, leading to an unprecedented expansion of urban land cover. Urban areas will become more important as carbon sinks and effective tools to assess carbon densities in these areas are therefore required. Using multi-scale LiDAR presents an opportunity to achieve this, providing a spatially explicit map of urban forest structure and AGB.

     

An Automatic Well Planner for Complex Well Trajectories

Mathematical Geosciences - A data-driven automatic well planner procedure is implemented to develop complex well trajectories by efficiently adapting to near-well reservoir properties and geometry....     

Nonlinear influences on ocean waves observed by X-band radar

The purpose of this study is to discuss the influence of signal nonlinearity upon X-band radar observations. A method for estimating the degree of nonlinearity by bispectral analysis was applied and discussed. We found that bispectral analyses from spatial radar backscatter series are similar to results obtained from water level time series. In addition, the average nonlinear degree from radar backscatter is related to wind speed. The accuracy of wave observations derived by consideration of the nonlinear effect from radar backscatter was also investigated. The estimated error in wave height from the radar data is also related to the degree of nonlinearity. In order to improve accuracy, the modulation transfer function method was applied in order to eliminate the influence of nonlinearity.     

Isolation of protein from Chlorella sorokiniana CY1 using liquid biphasic flotation assisted with sonication through sugaring-out effect

Microalgae, a sustainable source of multi beneficial components has been discovered and could be utilised in pharmaceutical, bioenergy and food applications. This study aims to investigate the sugaring-out effect on the recovery of protein from wet green microalga, Chlorella sorokiniana CY1 which was assisted with sonication. A comparison of monosaccharides and disaccharides as one of the phaseforming constituents shows that the monosaccharides, glucose was the most suitable sugar in forming the phases with acetonitrile to enhance the production of protein(52% of protein). The protein productivity of microalgae was found to be significantly influenced by the volume ratio of both phases, as the yield of protein increased to 77%. The interval time between the sonication as well as the sonication modes were influencing the protein productivity as well. The optimum protein productivity was obtained with 10 s of resting time in between sonication. Pulse mode of sonication was suitable to break down the cell wall of microalgae compared to continuous mode as a lower protein yield was obtained with the application of continuous mode. The optimum condition for protein extraction were found as followed: 200 g/L glucose as bottom phase with volume ratio of I:1.25, 10 s of resting time for ultrasonication, 5 s of ultrasonication in pulse mode and 0.25 g of biomass weight. The high yield of protein about 81% could be obtained from microalgae which demonstrates the potential of this source and expected to play an important role in the future.     

Local isotropy and anisotropy in the sheared and heated atmospheric surface layer

Longitudinal velocity and temperature measurements above a uniform dry lakebed were used to investigate sources of eddy-motion anisotropy within the inertial subrange. Rather than simply test the adequacy of locally isotropic relations, we investigated directly the sources of anisotropy. These sources, in a daytime desert-like climate, include: (1) direct interaction between the large-scale and small-scale eddy motion, and (2) thermal effects on the small-scale eddy motion. In order to explore these two anisotropy sources, we developed statistical measures that are sensitive to such interactions. It was found that the large-scale/small-scale interaction was significant in the inertial subrange up to 3 decades below the production scale, thus reducing the validity of the local isotropy assumption. The anisotropy generated by thermal effects was also significant and comparable in magnitude to the former anisotropy source. However, this thermal anisotropy was opposite in sign and tended to counteract the anisotropy generated by the large-scale/smallscale interaction. The thermal anisotropy was attributed to organized ramp-like patterns in the temperature measurements. The impact of this anisotropy cancellation on the dynamics of inertial subrange eddy motion was also considered. For that purpose, the Kolmogorov-Obukhov structure function equation, as derived from the Navier-Stokes equations for locally isotropic turbulence, was employed. The Kolmogorov-Obukhov structure function equation in conjunction with Obukhov's constant skewness closure hypothesis reproduced the measured second- and third-order structure functions. Obukhov's constant skewness closure scheme, which is also based on the local isotropy assumption, was verified and was found to be in good agreement with the measurements. The accepted 0.4 constant skewness value derived from grid turbulence experiments overestimated our measurements. A suggested 0.26 constant skewness value, which we derived from Kolmogorov's constant, was found to be adequate.     

Derivation Of Refractive Index And Temperature Gradients From Optical Scintillometry To Correct Atmospherically Induced Errors For Highly Precise Geodetic Measurements

Refraction effects of optical beams are generally caused by an inhomogeneous propagation medium and are a major source of systematic errors in the precise optical determination of angles and distances in the atmospheric surface layer. In this contribution a method for deriving vertical temperature and refractive index gradients from optical scintillation is presented. Knowledge of these gradients is required for the compensation of atmospherically induced errors for highly precise terrestrial geodetic measurements, like direct transfer and levelling. The advantage of the present optical method is, that temperature and refractive index gradients can be derived as line-averaged values over the propagation path, which is not possible by meteorological point measurements. Field observations have been carried out with a displaced-beam scintillometer over flat terrain and under different atmospheric conditions in order to verify this method. The experiments show, that this method allows to derive accurate correction values for precise terrestrial geodetic measurements.     

Killing of marine larvae by diesel oil

Fourteen species (five phyla) of pelagic larvae, when exposed to 0.5% of No. 2 diesel oil in sea water, survived from 3 to 72 hours excepting the larvae of the echinoderm Crossaster which lived for 8 days.     

Processing and quality control of flux data during LITFASS-2003

Different aspects of the quality assurance and quality control (QA/QC) of micrometeorological measurements were combined to create a comprehensive algorithm which was then applied to experimental data from LITFASS-2003 (Lindenberg Inhomogeneous Terrain—Fluxes between Atmosphere and Surface: a long term Study). Eddy-covariance measurements of the latent heat flux were the main focus of the QA/QC efforts. The results of a turbulence sensor intercomparison experiment showed deviations between the different eddy-covariance systems on the order of 15%, or less than 30 W m⁻², for the latent heat flux and 5%, or less than 10 W m⁻², for the sensible heat flux. In order to avoid uncertainties due to the post-processing of turbulence data, a comprehensive software package was used for the analysis of experimental data from LITFASS-2003, including all necessary procedures for corrections and quality control. An overview of the quality test results shows that for most of the days more than 80% of the available latent heat flux data are of high quality so long as there are no instrumental problems. The representativeness of a flux value for the target land-use type was analysed using a stochastic footprint model. Different methods to calculate soil heat fluxes at the surface are discussed and a sensitivity analysis is conducted to select the most robust method for LITFASS-2003. The lack of energy balance closure, which was found for LITFASS-2003, can probably be attributed to the presence of low-frequency flux contributions that cannot be resolved with an averaging time of 30 min. Though the QA/QC system has been developed for the requirements of LITFASS-2003, it can also be applied to other experiments dealing with similar objectives.     

The effect of mountainous topography on moisture exchange between the “surface” and the free atmosphere

Typical numerical weather and climate prediction models apply parameterizations to describe the subgrid-scale exchange of moisture, heat and momentum between the surface and the free atmosphere. To a large degree, the underlying assumptions are based on empirical knowledge obtained from measurements in the atmospheric boundary layer over flat and homogeneous topography. It is, however, still unclear what happens if the topography is complex and steep. Not only is the applicability of classical turbulence schemes questionable in principle over such terrain, but mountains additionally induce vertical fluxes on the meso-γ scale. Examples are thermally or mechanically driven valley winds, which are neither resolved nor parameterized by climate models but nevertheless contribute to vertical exchange. Attempts to quantify these processes and to evaluate their impact on climate simulations have so far been scarce. Here, results from a case study in the Riviera Valley in southern Switzerland are presented. In previous work, measurements from the MAP-Riviera field campaign have been used to evaluate and configure a high-resolution large-eddy simulation code (ARPS). This model is here applied with a horizontal grid spacing of 350 m to detect and quantify the relevant exchange processes between the valley atmosphere (i.e. the ground “surface” in a coarse model) and the free atmosphere aloft. As an example, vertical export of moisture is evaluated for three fair-weather summer days. The simulations show that moisture exchange with the free atmosphere is indeed no longer governed by turbulent motions alone. Other mechanisms become important, such as mass export due to topographic narrowing or the interaction of thermally driven cross-valley circulations. Under certain atmospheric conditions, these topographical-related mechanisms exceed the “classical” turbulent contributions a coarse model would see by several times. The study shows that conventional subgrid-scale parameterizations can indeed be far off from reality if applied over complex topography, and that large-eddy simulations could provide a helpful tool for their improvement.