Fast wavelet-based stochastic simulation using training images

Spatial uncertainty modelling is a complex and challenging job for orebody modelling in mining, reservoir characterization in petroleum, and contamination modelling in air and water. Stochastic simulation algorithms are popular methods for such modelling. In this paper, discrete wavelet transformation (DWT)-based multiple point simulation algorithm for continuous variable is proposed that handles multi-scale spatial characteristics in datasets and training images. The DWT of a training image provides multi-scale high-frequency wavelet images and one low-frequency scaling image at the coarsest scale. The simulation of the proposed approach is performed on the frequency (wavelet) domain where the scaling image and wavelet images across the scale are simulated jointly. The inverse DWT reconstructs simulated realizations of an attribute of interest in the space domain. An automatic scale-selection algorithm using dominant mode difference is applied for the selection of the optimal scale of wavelet decomposition. The proposed algorithm reduces the computational time required for simulating large domain as compared to spatial domain multi-point simulation algorithm. The algorithm is tested with an exhaustive dataset using conditional and unconditional simulation in two- and three-dimensional fluvial reservoir and mining blasted rock data. The realizations generated by the proposed algorithm perform well and reproduce the statistics of the training image. The study conducted comparing the spatial domain filtersim multiple-point simulation algorithm suggests that the proposed algorithm generates equally good realizations at lower computational cost.     

Invariant shape solutions of the spinning three craft Coulomb tether problem

In this paper we study shape-preserving formations of three spacecraft, where the formation keeping forces arise from the electric charges deposed on each craft. Inspired by Lagrange’s 3-body problem, the general conditions that guarantee preservation of the geometric shape of the electrically charged formation are derived. While the classical collinear configuration is a solution to the problem, the equilateral triangle configuration is found to only occur with unbounded relative motion. The three collinear spacecraft problem is analyzed and the possible solutions are categorized based on the spacecraft mass–charge ratio. Precise statements on the number of solutions associated with each category are provided. Finally, a methodology is proposed to study boundedness of the collinear solution that is inspired by past understanding and results for the 3-body problem. Given the initial position and the velocity vectors of each craft along with the charges, analytical solutions are provided describing the resulting relative motion.     

Effects of diagenesis on paleoclimate reconstructions from modern and young fossil corals

The integrity of coral-based reconstructions of past climate variability depends on a comprehensive knowledge of the effects of post-depositional alteration on coral skeletal geochemistry. Here we combine millimeter-scale and micro-scale coral Sr/Ca data, scanning electron microscopy (SEM) images, and X-ray diffraction with previously published δ¹⁸O records to investigate the effects of submarine and subaerial diagenesis on paleoclimate reconstructions in modern and young sub-fossil corals from the central tropical Pacific. In a 40-year-old modern coral, we find secondary aragonite is associated with relatively high coral δ¹⁸O and Sr/Ca, equivalent to sea-surface temperature (SST) artifacts as large as −3 and −5 °C, respectively. Secondary aragonite observed in a 350-year-old fossil coral is associated with relatively high δ¹⁸O and Sr/Ca, resulting in apparent paleo-SST offsets of up to −2 and −4 °C, respectively. Secondary Ion Mass Spectrometry (SIMS) analyses of secondary aragonite yield Sr/Ca ratios ranging from 10.78 to 12.39 mmol/mol, significantly higher compared to 9.15 ± 0.37 mmol/mol measured in more pristine sections of the same fossil coral. Widespread dissolution and secondary calcite observed in a 750-year-old fossil coral is associated with relatively low δ¹⁸O and Sr/Ca. SIMS Sr/Ca measurements of the secondary calcite (1.96-9.74 mmol/mol) are significantly lower and more variable than Sr/Ca values from more pristine portions of the same fossil coral (8.22 ± 0.13 mmol/mol). Our results indicate that while diagenesis has a much larger impact on Sr/Ca-based paleoclimate reconstructions than δ¹⁸O-based reconstructions at our site, SIMS analyses of relatively pristine skeletal elements in an altered coral may provide robust estimates of Sr/Ca which can be used to derive paleo-SSTs.     

Curie Point Depth Estimates from Aeromagnetic Data from Death Valley and Surrounding Regions, California

Aeromagnetic data were analyzed to determine the Curie point depth (CPD) by power density spectral and three-dimensional inversion methods within and surrounding Death Valley in southern California. We calculated the CPD for 0.5° regions using 2D power density spectral methods and found that the CPDs varied between 8 and 17 km. However, the 0.5° region may average areas that include shallow and deep CPDs, and because of this limitation, we used the 3D inversion method to determine if this method may provide better resolution of the CPDs. The final 3D model indicates that the depth to the bottom of the magnetic susceptible bodies varies between 5 and 23 km. Even though both methods produced roughly similar results, the 3D inversion method produced a higher lateral resolution of the CPDs. The shallowest CPDs occur within the central and southern Death Valley, Panamint Valley, Coso geothermal field and the Tecopa hot springs region. Deeper (>15 km) CPDs occur over outcropping granitic and Precambrian lithologies in the Panamint Range, Grapevine Mountains, Black Mountains and the Argus Range. The shallowest CPD occurs within the central Death Valley over a possible seismically imaged magma body and slightly deeper values occur within the Panamint Valley, southern Death Valley and Tecopa Hot Springs. The shallow CPD values suggest that partially molten material may also be found in these latter regions. The CPD computed heat flow values for the region suggest that the entire area has high heat flow values (>100 mW m^?2), on the other hand, locally extremely high values (>200 mW m^?2) occur within the Panamint Valley, the southern and central Death Valley and Tecopa Hot Springs region. These locally high heat flow values may be related to midcrustal magma bodies; but additional geophysical experiments are needed to determine if the magma bodies exist.     

Grain Dust as an Energy and Food Resource

Grain dust is a by-product produced in large quantities from grain storage and milling facilities but as of yet is treated as waste locally and regionally. Based on the conviction that material properties determination is a necessary first step for the efficient handling, processing, and use of materials, this work was initiated with the purpose of characterizing grain dust in terms of its physicochemical and thermal properties. Grain dust as a special waste material is a complex, dry, grey powder, composed mainly of grain chips, soil, whole wheat and barley seeds, and straw. The material was found to have a bulk density of 362 kg/m³, specific heat, thermal conductivity, and gross calorific value were 1540 J/kg K, 0.27 W/m K, and 4120 kJ/kg, respectively. Nutritional composition showed low concentrations of protein, and ether extract, and high concentrations of carbohydrates, ash, and crude fiber. Mineral analysis showed that Ca, Fe, and K were available in significant amounts ranging from 30 to 135 ppm. The as-received material is dry enough to be stored very safely and burned for energy production in either its loose form or as pellets. Moreover, adjustments to the ash content and C:N ratio of the material deems it a promising candidate as an animal feed and as a substrate for composting and biogasification.     

Seismic hazard assessment and rheological implications: a case study selected for cities of Saudi Arabia along the eastern coast of Red Sea

A probabilistic approach is used to evaluate the seismic hazard for 12 strategic cities in Saudi Arabia along the eastern coast of Red Sea. The focal depth variations controlled by rheological characteristics are taken into account for hazard calculations, and its creditability is tested through sensitivity analysis for hazard results. This study presents a neo-probabilistic seismic hazard assessment methodology in which the focal depth distribution of earthquakes within seismogenic layer is divided into three depth slices. These depth slices are based upon rheological characteristic of seismogenic layer. The hazard results are obtained using this depth-slice methodology and conventional approach in which uniform distribution of seismicity within seismogenic layer is assumed. The sensitivity analysis culminated in underestimation of hazard values in higher frequencies for uniform distribution of seismicity within seismogenic layer. Foregoing the observations recorded above, it can be concluded that the exploitation of depth-slices biased by the rheology to calculate hazard is relatively preferable in the situations demanding safety measures.     

Imaging fracture distributions of the Al-Khuff Formation outcrops using GPR and ERT geophysical techniques,Al-Qassim area,Saudi Arabia

Khuff Formation is of utmost importance in Saudi Arabia for oil and gas reservoir although it is composed mainly of limestone. This reason refers to the existence of intensive fractures which play a vital role in the increase in porosity and permeability of this formation. The fracture pattern in the study area was verified through 2D and 3D ground penetrating radar (GPR)-defined and electrical resistivity tomography (ERT)-defined surveys. In this respect, ten of 2D GPR surveys were collected along an intersected grid of profiles covering the study area while ERT data were collected along three profiles of the GPR grid. The results were interpreted in light of the field-based structural and stratigraphic assessment of the outcropping rocks. The analysis of the inverted ERT and filtered GPR sections revealed the presence of fractures. Several resistivity and electromagnetic reflection anomalies were laterally and vertically identified across the measured sections clarifying fractures that extend to a depth of 24 m in the limestone. Most fractures are oriented vertical to sub-vertical dipping both east-west and north-south.     

Groundwater development,Kharga Oases,Western Desert of Egypt: A long-term environmental concern

The groundwater reserves in Kharga Oases have been studied for the long-term socioeconomic development in the area. The Nubian Sandstone, which consists of a thick sequence of coarse clastic sediments of sandstone, sandy clay interbedded with shale, and clay beds, forms a complex aquifer system. The Nubian Aquifer has been providing water to artesian wells and springs in the Kharga Oases for several thousand years. Groundwater in the Kharga Oases is withdrawn from springs and shallow and deep artesian wells Nearly all the wells originally flowed, but with the exploitation of ground-water from deep wells for irrigation beginning about 1959. the natural flows declined as more and more closely spaced deep wells were drilled By 1975 many deep wells had ceased to flow The water demand in the area has been met by pumping both shallow and deep wells The total annual extraction from deep wells has fluctuated over the year, however, the annual withdrawal from deep wells has exceeded extraction from shallow wells About 17 billion m³ of water was withdrawn from the combination of shallow and deep wells during the period 1960–1980 The Nubian complex aquifer in the Kharga Oases has a very large groundwater potential that could be exploited and beneficially used for a long-term agricultural development in the area, provided proper well spacing and management are implemented Other major environmental considerations for which precise hydrogeologic data are needed include

1. Determination of the long-term yield available from properly constructed and producing artesian wells that will support a planned migration of population from the overcrowded Nile delta and flood plain areas
2. Development of an effective management program and adequate staff to maintain groundwater production over an extended period of years
3. The impact on climate caused by extensive irrigation in the oases of the Western Desert of Egypt
4. Protection against water logging of soils from irrigation practices
5. Protection against salinization of soils from irrigation practices
6. Development of effective surface and subsurface drainage practices
7. The impact of farming and pest control practices on the shallow groundwater of the oases
8. Determination of the long-term development of the artesian water on the quality of the water from the aquiter systems in the Western Desert

This paper addresses items 1, 2 and 8.