A mathematical programming approach for scheduling equipment in a surface coal mining operation

Summary A mathematical programming model for scheduling open pit mining was developed and validated using data from a surface mining operation. A two-phase solution procedure was used involving repeated evaluations of an integer scheduling model and a simple transportation model.     

Evaluation of biotite-garnet geothermometers: application to the English River subprovince,Ontario

In recent years, there have been numerous calibrations of the biotite-garnet Fe-Mg exchange geothermometer. The Eastern Lac Seul region of the English River subprovince, Ontario, provides an excellent field area in which to compare these calibrations.Trend surface analysis using the temperatures obtained from garnet cores and matrix biotites-showed almost identical trends in the eastern Lac Seul region regardless of the calibration used. The absolute temperatures and the precision of each calibration do, however, show large variation. Geothermometers based solely on lnK_D were found to give more precise results than the calibrations that attempt to incorporate non-Fe-Mg components. The Perchuk and Lavrent'eva (1983) thermometer yields the most precise and accurate results. If a sufficient number of samples are collected over a region, it can be used to estimate metamorphic temperature trends to ±30° C. Metamorphism and migmatization of the English River subprovince occurred during the Kenoran orogeny, 2.68 b.y. ago. Our results show that a thermal anticline has been preserved, with temperatures of 600° C at the north and south contacts with Uchi and Wabigoon Greenstone belts, increasing to 725 °C at the center of the subprovince. A garnet-cordierite in isograd occurs at 650° C and an orthopyroxene in isogradat 700° C.     

The control of maturity and kerogen type on quantitative analytical pyrolysis data

Quantitative pyrolysis-gas chromatography has been performed on 96 kerogen samples isolated from 17 wells on the Norwegian Continental shelf. Petrographic and bulk geochemical measurements were also performed on the samples, and a combined data set of 117 variables for each sample was analysed using principal components analysis (PCA). This approach provides an objective and reproducible means of kerogen characterisation, which can be easily automated. In addition to objective kerogen characterisation and facile visualisation of facies and maturity related chemical trends, the method has the potential to allow objective prediction of key geochemical parameters such as maturity level from pyrogram data.     

Chemical models of vitrinite reflectance evolution

A global kinetic model of vitrinite reflectance evolution is developed to illustrate a new approach to vitrinite reflectance modelling. The method utilizes quantitative pyrolysisgas Chromatographic data from isolated vitrinite kerogens to derive the concentrations of structurally specific moieties (alkylphenol precursors) in vitrinite kerogens as a function of rank. This data set is used to calibrate a gaussian distributed activation energy model for phenol-precursor loss from vitrinites using published burial histories. Phenol yield is then correlated with vitrinite reflectance. The predictive vitrinite reflectance model based on this phenol system is tested against a well from the US Gulf Coast area, the results being compared with a Lopatin model of the same well. This type of vitrinite reflectance model is numerically compatible with current generation kinetic models of hydrocarbon production from source rock kerogens, and has a better heating rate independence than Lopatin's method.

---

Zusammenfassung Um eine neue Annäherung an das Modellieren der Vitrinit-Reflexion zu demonstrieren, wurde ein global-kinetisches Modell der Entwicklung der Vitrinit-Reflexion erstellt. Die Methode verwendet Daten von quantitativer Pyrolyse-Gas Chromatographie isolierten Vitrinitkerogens, um die Konzentration von strukturspezifischen Moieten (Alkylophenol-Vorläufer) in Vitrinitkerogen als Funktion der Reihenfolge zu erhalten. Verwendet werden diese Daten, um ein nach Gauss verteiltes Aktivierungsenergiemodell zu kalibrieren, das mit Hilfe von veröffentlichten Ablagerungsgeschichten den Phenol-Vorläuferverlust von Vitriniten beschreibt. Die Phenol-Ausbeute wird dann mit der Vitrinit-Reflexion korreliert. Das auf dem Phenolsystem basierende Vitrinit-Reflexionsmodell wird in einem Test einem Bohrloch an der amerikanischen Golfküste gegenübergestellt. Die Ergebnisse werden mit dem Lopatinmodell des selben Bohrlochs verglichen. Diese Art eines Vitrinit-Reflexionsmodells ist numerisch wettbewerbsfähig mit der heutigen Generation kinetischer Modelle der Kohlenwasserstofferzeugung aus Trägergesteinskerogenen. Außerdem ist seine Unabhängigkeit von der Erhitzungsrate größer als mit der Lopatin-Methode.

---

Résumé L'auteur présente un modèle cinétique global de l'évaluation du pouvoir réflecteur de la vitrinite, afin d'illuster une approche nouvelle dans la modélisation de ce pourvoir réflecteur. Cette méthode utilise des données quantitatives fournies par la pyrolyse combinée à la Chromatographie en phase gazeuse et relatives à des kérogènes de vitrinite. A partir de ces données, on exprime les concentrations des fractions à structure spécifique dans les kérogènes vitrinitiques en fonction du degré d'évolution. A partir de cet ensemble de données, on établit un modèle de l'énergie d'activation selon une distribution gaussienne, applicable à la perte de phénol-précurseur dans les vitrinites, en se référant à cet effet aux déroulements historiques de l'enfouissement fournis par la littérature. La production de phénol est alors mise en relation avec le pourvoir réflecteur de la vitrinite. Le modèle du pouvoir réflecteur prévisible, basé sur ce système des phénols a été contrôlé dans un forage de la côte du golfe du Mexique (USA), et les résultats comparés à ceux du modèle de Lopatin appliqué au même forage. Le modèle présenté est numériquement compatible avec les modèles cinétiques habituels qui rendent compte de la production d'hydrocarbure à partir des kérogènes des roches-mères; il présente une meilleure indépendance que la méthode de Lopatin par rapport au taux d'échauffement.

---

, , , .. , , - , . , . . , . ; ; .

     

Cloud detection and analysis: A review of recent progress

The major types of cloud retrieval algorithms are reviewed with special emphasis being placed upon recent (i.e., post 1981/1982) developments and novel techniques. Satellite-based retrieval algorithms can be grouped into three classes: threshold methods, statistical procedures, and radiative transfer techniques, although each algorithm depends upon implicit, if not overt, inversion of the radiative transfer equation. The fourth type of retrieval differs very considerably from the satellite-based techniques as it depends upon surface-based measurements which are generally, but not always, human rather than computer based. There is a tendency to assume that surface-based observations, especially of total cloud amount, are “correct” but that they may differ from satellite-based retrievals because of the differences in viewing geometry. Actually the literature reveals surprisingly few intercomparison studies. None of the satellite-based techniques have yet been well-validated in a variety of situations, while surface-based observations are made in all terrain and climate regimes. Overall there seems to be good reason to believe that surface-based observations of cloud amount and especially of low cloud amount, character, and base height, can add significantly to satellite-based global nephanalyses now operational or currently being planned.     

Inversion for seismic anisotropy using genetic algorithms1

A general inversion scheme based on a genetic algorithm is developed to invert seismic observations for anisotropic parameters. The technique is applied to the inversion of shear-wave observations from two azimuthal VSP data sets from the Conoco test site in Oklahoma. Horizontal polarizations and time-delays are inverted for hexagonal and orthorhombic symmetries. The model solutions are consistent with previous studies using trial and error matching of full waveform synthetics. The shear-wave splitting observations suggest the presence of a shear-wave line singularity and are consistent with a dipping fracture system which is known to exist at the test site. Application of the inversion scheme prior to full waveform modelling demonstrates that a considerable saving in time is possible whilst retaining the same degree of accuracy.     

A three-dimensional,iterative mapping procedure for the implementation of an ionosphere-magnetosphere anisotropic Ohm’s law boundary condition in global magnetohydrodynamic simulations

The mathematical formulation of an iterative procedure for the numerical implementation of an ionosphere-magnetosphere (IM) anisotropic Ohm’s law boundary condition is presented. The procedure may be used in global magnetohydrodynamic (MHD) simulations of the magnetosphere. The basic form of the boundary condition is well known, but a well-defined, simple, explicit method for implementing it in an MHD code has not been presented previously. The boundary condition relates the ionospheric electric field to the magnetic field-aligned current density driven through the ionosphere by the magnetospheric convection electric field, which is orthogonal to the magnetic field B, and maps down into the ionosphere along equipotential magnetic field lines. The source of this electric field is the flow of the solar wind orthogonal to B. The electric field and current density in the ionosphere are connected through an anisotropic conductivity tensor which involves the Hall, Pedersen, and parallel conductivities. Only the height-integrated Hall and Pedersen conductivities (conductances) appear in the final form of the boundary condition, and are assumed to be known functions of position on the spherical surface R=R₁ representing the boundary between the ionosphere and magnetosphere. The implementation presented consists of an iterative mapping of the electrostatic potential , the gradient of which gives the electric field, and the field-aligned current density between the IM boundary at R=R₁ and the inner boundary of an MHD code which is taken to be at R₂>R₁. Given the field-aligned current density on R=R₂, as computed by the MHD simulation, it is mapped down to R=R₁ where it is used to compute by solving the equation that is the IM Ohm’s law boundary condition. Then is mapped out to R=R₂, where it is used to update the electric field and the component of velocity perpendicular to B. The updated electric field and perpendicular velocity serve as new boundary conditions for the MHD simulation which is then used to compute a new field-aligned current density. This process is iterated at each time step. The required Hall and Pedersen conductances may be determined by any method of choice, and may be specified anew at each time step. In this sense the coupling between the ionosphere and magnetosphere may be taken into account in a self-consistent manner.