Adjoint formulation and constraint handling for gradient-based optimization of compositional reservoir flow

An adjoint formulation for the gradient-based optimization of oil–gas compositional reservoir simulation problems is presented. The method is implemented within an automatic differentiation-based compositional flow simulator (Stanford’s Automatic Differentiation-based General Purpose Research Simulator, AD-GPRS). The development of adjoint procedures for general compositional problems is much more challenging than for oil–water problems due to the increased complexity of the code and the underlying physics. The treatment of nonlinear constraints, an example of which is a maximum gas rate specification in injection or production wells, when the control variables are well bottom-hole pressures, poses a particular challenge. Two approaches for handling these constraints are presented—a formal treatment within the optimizer and a simpler heuristic treatment in the forward model. The relationship between discrete and continuous adjoint formulations is also elucidated. Results for four example cases of increasing complexity are presented. Improvements in the objective function (cumulative oil produced) relative to reference solutions range from 4.2 to 11.6 %. The heuristic treatment of nonlinear constraints is shown to offer a cost-effective means for obtaining feasible solutions, which are, in some cases, better than those obtained using the formal constraint handling procedure.     

Spread-F during the magnetic storm of 22 January 2004 at low latitudes: Effect of IMF-Bz in relation to local sunset time

The paper describes the results of spread-F at low latitude stations around the world during the magnetic storm starting at 0130 UT on 22 January 2004. The storm can be divided into two phases, first phase up to 1000 UT when interplanetary magnetic field IMF-Bz was highly fluctuating around a small positive value and the second phase after a sudden large southward turning of IMF-Bz at 1030 UT. The first phase produced strong spread-F at Jicamarca, Sao Luis, and Ascension Island and caused complete inhibition of spread-F at Thumba and Waltair in India. It generated weak spread-F at Ho Chi Minh City in Vietnam and strong spread-F at Hainan and Chung Li. The strong spread-F at Hainan and Chung Li were caused by the positive IMF-Bz during the first phase of the storm and not by the negative pulse of IMF-Bz at 1000 UT.     

Microphysical properties of stratocumulus clouds during ACE‐2

Microphysical measurements performed during 8 flights of the CLOUDYCOLUMN component of ACE‐2, with the Meteo‐France Merlin‐IV, are analyzed in terms of droplet number concentration and size. The droplet concentration is dependent upon the aerosol properties within the boundary layer. Its mean value over a flight varies from 55 cm⁻³, for the cleanest conditions, to 244 cm⁻³, for the most polluted one. For each flight, the variability of the concentration, in selected cloud regions that are not affected by mixing with dry air or drizzle scavenging, ranges from 0.5 to 1.5 of the mean value. The mean volume diameter increases with altitude above cloud base according to the adiabatic cloud model. The frequency distribution of mean droplet volume normalized by the adiabatic value, for the selected regions, shows the same dispersion as the distribution of normalized concentration. The values of droplet concentration versus mean volume diameter are then examined in sub‐adiabatic samples to characterize the effects of mixing and drizzle scavenging. Finally, the ratio of mean volume diameter to effective diameter is analyzed and a simple relationship between these 2 crucial parameters is proposed.     

Glacial-interglacial mean sea level pressure change due to sea level, ice sheet and atmospheric mass changes

The change in the global mean atmospheric pressure between glacial and interglacial periods is evaluated at sea level. This change originates in a modification of topography and in a possible variation in the atmospheric mass. In this calculation the atmosphere is at hydrostatic equilibrium, and the parameters describing the glacial period are varied in a plausible range. The result, with constant atmospheric mass, is a mean sea level pressure decrease of 9–15 hPa linked with the deglaciation. The corresponding pressure change at the reference level corresponding to the present day sea level does not exceed one hPa. When considering only the change in the atmospheric mass, an increase which does not exceed 2 hPa is found, linked with the deglaciation.     

New 1982 — 1990 photometry of λ Andromedae and its 11-year cycle

We present photoelectric photometry of λ And never before published, obtained between February 1982 and December 1990 at 29 different observatories. Then we combine it with all other photometry available to us (previously published, contained in the I.A.U. Commission 27 Archives, and obtained with the Vanderbilt 16-inch automatic telescope but not yet published), to yield a 14.8-year data base. Analysis reveals a long-term cycle in mean brightness, with a full range of 0^m.15 and a period of 11.4 ± 0.4 years. Because most of our new photometry was concentrated in the 1983-84 observing season, we analyze that one well-defined light curve with a two-spot model. Spot A keeps a 0^m.04 amplitude throughout four rotation cycles whereas the amplitude of spot B diminishes from 0^m.09 down almost to 0^m.03. The spot rotation periods were 55^d.9 ± 0^d.6 and 520^d.8 ± 1^d.0, respectively.     

http://www.sciencedirect.com/science/article/pii/S167498711400108X

We present three 3D numerical models of deep subduction where buoyant material from an oceanic plateau and a plume interact with the overriding plate to assess the influence on subduction dynamics,trench geometry,and mechanisms for plateau accretion and continental growth.Transient instabilities of the convergent margin are produced,resulting in:contorted trench geometry;trench migration parallel with the plate margin;folding of the subducting slab and orocline development at the convergent margin;and transfer of the plateau to the overriding plate.The presence of plume material beneath the oceanic plateau causes flat subduction above the plume,resulting in a "bowed" shaped subducting slab.In plateau-only models,plateau accretion at the edge of the overriding plate results in trench migration around the edge of the plateau before subduction is re-established directly behind the trailing edge of the plateau.The plateau shortens and some plateau material subducts.The presence of buoyant plume material beneath the oceanic plateau has a profound influence on the behaviour of the convergent margin.In the plateau + plume model,plateau accretion causes rapid trench advance.Plate convergence is accommodated by shearing at the base of the plateau and shortening in the overriding plate.The trench migrates around the edge of the plateau and subduction is re-established well behind the trailing edge of the plateau,effectively embedding the plateau into the overriding plate.A slab window forms beneath the accreted plateau and plume material is transferred from the subducting plate to the overriding plate through the window.In all of the models,the subduction zone maintains a relatively stable configuration away from the buoyancy anomalies within the downgoing plate.The models provide a dynamic context for plateau and plume accretion in Phanerozoic accretionary orogenic systems such as the East China Orogen and the Central Asian Orogen(Altiads),which are characterised by accreted ophiolite complexes with diverse geochemical affinities,and a protracted evolution of accretion of exotic terranes including oceanic plateau and terranes with plume origins.     

NIGHTGLOW: an instrument to measure the Earth’s nighttime ultraviolet glow—results from the first engineering flight

We have designed and built an instrument to measure and monitor the “nightglow” of the Earth’s atmosphere in the near ultraviolet (NUV). In this paper we describe the design of this instrument, called NIGHTGLOW. NIGHTGLOW is designed to be flown from a high altitude research balloon, and circumnavigate the globe. NIGHTGLOW is a NASA, University of Utah, and New Mexico State University project. A test flight took place from Palestine, Texas on July 5, 2000, lasting about 8 h. The instrument performed well and landed safely in Stiles, Texas with little damage. The resulting measurements of the NUV nightglow are compared with previous measurements from sounding rockets and balloons.     

Fossil hyrax dung and evidence of Late Pleistocene and Holocene vegetation types in the Namib Desert

Pollen was derived from fossil dung of herbivorous hyraxes, deposited in a rock shelter on the highest mountain in Namibia, Dâures or Brandberg, situated on the Namib Desert margin. Radiocarbon dating ranging in age between modern times and 30 000 yr BP showed it represents the first empirical pollen evidence of continental palaeovegetation during the Late Pleistocene along the western escarpment of southern Africa. The initial results indicate Last Glacial Maximum vegetation differed totally from the current pattern as vegetation types were dominated by small Asteraceae shrubs, in contrast to those of the Holocene and modern times which show more succulents, grass and woody elements (arboreal pollen). The results suggest that Cape floral communities did not reach into the tropics along the western escarpment of Africa, despite such pollen types occurring in marine cores. Copyright © 2004 John Wiley & Sons, Ltd.     

THE FREE-WILSON PARADIGM REDUX:SIGNIFICANCE OF THE FREE-WILSON COEFFICIENTS,INSIGNIFICANCE OF COEFFICIENT ‘UNCERTAINTIES’ AND STATISTICAL SINS

The Free-Wilson paradigm is an established and powerful tool for quantitatively relating activity withchemical structure.Current implementations of the paradigm,however,are flawed both conceptually andin execution.As part of an attempt to more fully realize the promise of the paradigm,it was necessaryto examine these limitations in detail.This report introduces a robust,theory-founded Free-Wilson implementation:stepwise principalcomponents regression analysis(SPCRA).SPCRA is computationally superior to previousimplementations but does not in itself correct their conceptual flaws.The development of SPCRA did,however,facilitate derivation of a simple and chemically significantinterpretation of the Free-Wilson structure-activity model.A number of statistical aspects of this modelcommonly misused in previous applications are discussed at length.These discussions provide criticalbackground for the development of an alternative implementation of the Free-Wilson paradigm.