Generation of a stochastic precipitation model for the tropical climate

Theoretical and Applied Climatology - A tropical country like Malaysia is characterized by intense localized precipitation with temperatures remaining relatively constant throughout the year. A...     

Kernel Principal Component Analysis for Efficient,Differentiable Parameterization of Multipoint Geostatistics

This paper describes a novel approach for creating an efficient, general, and differentiable parameterization of large-scale non-Gaussian, non-stationary random fields (represented by multipoint geostatistics) that is capable of reproducing complex geological structures such as channels. Such parameterizations are appropriate for use with gradient-based algorithms applied to, for example, history-matching or uncertainty propagation. It is known that the standard Karhunen–Loeve (K–L) expansion, also called linear principal component analysis or PCA, can be used as a differentiable parameterization of input random fields defining the geological model. The standard K–L model is, however, limited in two respects. It requires an eigen-decomposition of the covariance matrix of the random field, which is prohibitively expensive for large models. In addition, it preserves only the two-point statistics of a random field, which is insufficient for reproducing complex structures. In this work, kernel PCA is applied to address the limitations associated with the standard K–L expansion. Although widely used in machine learning applications, it does not appear to have found any application for geological model parameterization. With kernel PCA, an eigen-decomposition of a small matrix called the kernel matrix is performed instead of the full covariance matrix. The method is much more efficient than the standard K–L procedure. Through use of higher order polynomial kernels, which implicitly define a high-dimensionality feature space, kernel PCA further enables the preservation of high-order statistics of the random field, instead of just two-point statistics as in the K–L method. The kernel PCA eigen-decomposition proceeds using a set of realizations created by geostatistical simulation (honoring two-point or multipoint statistics) rather than the analytical covariance function. We demonstrate that kernel PCA is capable of generating differentiable parameterizations that reproduce the essential features of complex geological structures represented by multipoint geostatistics. The kernel PCA representation is then applied to history match a water flooding problem. This example demonstrates that kernel PCA can be used with gradient-based history matching to provide models that match production history while maintaining multipoint geostatistics consistent with the underlying training image.     

Monazite chemistry: evidence for a prehistoric journey

The chemistry of monazite from the black sand deposits of northern Sinai beach suggests that pegmatites and granites of the Eastern Desert of Egypt are the most likely source rocks. The floods associated with the pluvial periods prevailed in Egypt during the Pleistocene were able to erode the source rocks and liberate heavy minerals including monazite. The mineral grains were moved through several wadis and tributaries to the River Nile. At the confluence sites, these heavy minerals were mixed with the Ethiopian and central African heavy mineral assemblages. The grains continued to move together downriver until being deposited in their current locations.     

Integration of remote sensing data with the field and laboratory investigation for lithological mapping of granitic phases: Kadabora pluton, Eastern Desert, Egypt

In the current study, an integration of Enhanced Thematic Mapper Plus (ETM+), field, and laboratory data have been used for lithological mapping of different granitic phases in the Kadabora area, Eastern Desert, Egypt. Application of enhancement techniques, including a new proposed band ratio combination (ratio 5/3, 3/1, 7/5 in RGB, respectively) and supervised classification images are used in discriminating different granitic phases in the Kadabora pluton from each other and from their environs. The data have been proved with the help of field and geochemical investigations. The results revealed that: (1) the Kadabora granitic pluton could be distinguished into three phases that recognized by field and laboratory investigation including granodiorite (phase I), monzogranite (phase II), and syeno-alkali feldspar granite (phase III). These phases are arranged according to their relative ages while the country rocks include ophiolitic mélange and metagabbro–diorite complex. It is also confirmed that the granitic pluton is invaded by dyke swarms which is trending in N–S direction. Geochemically, results show that the granodiorite is calc-alkaline, I-type and formed under subduction tectonic regime. Monzogranite falls within the alkaline and highly fractionated calc-alkaline granites, whereas syeno-alkali feldspar granite extends into proper alkaline granitoids field. Monzogranite and syeno-alkali feldspar granite belong to the A₂-subtype granite. This A₂-subtype granite was probably formed in an extensional regime, subsequent to subduction which can lead to tensional break-up of the crust (i.e., post-collisional, post-orogenic granites). The monzogranite and the syeno-alkali feldspar granite were probably formed by partial melting of relatively anhydrous lower crust source and/or tonalite to granodiorite is viable alternative to the granulite source.     

Gradient-based Pareto optimal history matching for noisy data of multiple types

The advantages of the simultaneous integration of production and time-lapse seismic data for history matching have been demonstrated in a number of previous studies. Production data provide accurate observations at particular spatial locations (wells), while seismic data enable global, though filtered/noisy, estimates of state variables. In this work, we present an efficient computational tool for bi-objective history matching, in which data misfits for both production and seismic measurements are minimized using an adjoint-gradient approach. This enables us to obtain a set of Pareto optimal solutions defining the optimal trade-off between production and seismic data misfits (which are, to some extent, conflicting). The impact of noise structure and noise level on Pareto optimal solutions is investigated in detail. We discuss the existence of the “best” trade-off solution, or least-conflicting posterior model, which corresponds to the history-matched model that is expected to provide the least-conflicting forecast of future reservoir performance. The overall framework is successfully applied in 2D and 3D compositional simulation problems to provide a single least-conflicting posterior model and, for the 2D case, multiple samples from the posterior distribution using the randomized maximum likelihood method.     

Relativistic and the first sectorial harmonics corrections in the critical inclination

The problem of the critical inclination is treated in the Hamiltonian framework taking into consideration post-Newtonian corrections as well as the main correction term of sectorial harmonics for an earth-like planet. The Hamiltonian is expressed in terms of Delaunay canonical variables. A canonical transformation is applied to eliminate short period terms. A modified critical inclination is obtained due to relativistic and the first sectorial harmonics corrections.     

消除滑动扫描地震数据中谐波干扰的反相关方法(英文)

滑动扫描技术是高效、高保真、环保的可控震源勘探技术之一,是下一组震源不必等待上一组震源震动结束即可开始震动的高效采集方法。该技术由于缩短了相邻两炮的等待时间,使得生产效率得到显著提高。但是后一炮的谐波畸变与前一炮的基波信号混叠在一起,不易分离,在相关后的地震记录上形成了严重的谐波干扰,降低了地震资料的质量。本文提出一种反相关方法来压制滑动扫描地震数据中的谐波干扰。该方法首先把地面力信号分解为基波和各阶谐波分量;然后将后一炮的相关前数据分别与各分量相关,只选取正时间轴中对应分量的自相关部分,利用各分量的反相关算子提取各阶谐波信息;最后从前一炮数据中减去提取出的高阶谐波,得到压制谐波后的地震记录。该方法对有效信号影响小,可同时处理相关前和相关后数据,而且算法简单稳定,计算效率高。本文分别对理论模型和实际数据进行处理,验证了该方法消除谐波干扰的有效性。     

Effect of Ozone and Ozone/Fenton in the Advanced Oxidation Process on Biodegradable Characteristics of Semi‐aerobic Stabilized Leachate

This study investigated the effects of O₃ and O₃/H₂O₂/Fe²⁺ in the advanced oxidation processes (AOPs) on the biodegradable and soluble characteristics of semi‐aerobic stabilized solid waste leachate. The biodegradability (BOD₅/chemical oxygen demand, COD) ratio improved from 0.034 to 0.05 and 0.1 following O₃ and O₃/H₂O₂/Fe²⁺, respectively. Fractions of biodegradable COD_(bi) (24%), non‐biodegradable COD_(ubi) (76%), soluble COD_(s) (59%), biodegradable soluble COD_(bsi) (38%), non‐biodegradable soluble COD_(ubsi) (62%), and particulate COD (PCOD) (41%) in stabilized leachate were also investigated. The fraction of COD_(bi) increased to 28 and 36% after applying O₃ and O₃/AOPs, respectively. COD_(S) increased to 59% after O₃ and to 72% after O₃/AOPs, whereas COD_(bsi) increased to 38 and 51% after O₃ and O₃/AOPs, respectively. The removal efficiency of COD_(S) was obtained at 5% after O₃ alone and improved to 51% following ozone‐based AOPs, whereas the removal efficiency of PCOD improved from 25% after O₃ to 71% after ozone‐based AOPs.     

Preliminary evidence of a link between surface soil properties and the arsenic content of shallow groundwater in Bangladesh

The extremely heterogeneous distribution of As in Bangladesh groundwater has hampered efforts to identify with certainty the mechanisms that lead to extensive mobilization of this metalloid in reducing aquifers. We show here on the basis of a high-resolution transect of soil and aquifer properties collected in Araihazar, Bangladesh, that revealing tractable associations between As concentrations in shallow (< 20 m) groundwater with other geological, hydrological, and geochemical features requires a lateral sampling resolution of 10–100 m. Variations in the electromagnetic conductivity of surface soils (5–40 mS/m) within a 500 m × 200 m area are documented with 560 EM31 measurements. The results are compared with a detailed section of groundwater As concentrations (5–150 μg/L) and other aquifer properties obtained with a simple sampling device, “the needle-sampler”, that builds on the local drilling technology. By invoking complementary observations obtained in the same area and in other regions of Bangladesh, we postulate that local groundwater recharge throughout permeable sandy soils plays a major role in regulating the As content of shallow aquifers by diluting the flux of As released from reducing sediments.