Dynamic data integration for structural modeling: model screening approach using a distance-based model parameterization

This paper proposes a novel history-matching method where reservoir structure is inverted from dynamic fluid flow response. The proposed workflow consists of searching for models that match production history from a large set of prior structural model realizations. This prior set represents the reservoir structural uncertainty because of interpretation uncertainty on seismic sections. To make such a search effective, we introduce a parameter space defined with a “similarity distance” for accommodating this large set of realizations. The inverse solutions are found using a stochastic search method. Realistic reservoir examples are presented to prove the applicability of the proposed method.     

非均匀月壤介质的被动微波辐射传输模拟

基于非均匀月壤物理模型和辐射传输方程,模拟月壤介质中的微波辐射传输特性,探讨频率、月壤厚度等与月表亮温的关系。结果表明:在低频段,月壤微波辐射亮温的动态变化范围较大,可探测的月壤厚度大,3 GHz时的最大可探测月壤厚度达12.4 m;在高频段对应的可探测月壤厚度较小,特别是从50GHz往后的频率段内,最大可探测月壤厚度均小于2 m。不同频率的亮温-厚度变化曲线没有交叉点,且频率越高,所能探测的月壤厚度越小。根据模拟结果,建立了月壤厚度与亮温的查找表。基于查找表,利用单个波段的亮温数据即可得到月壤厚度信息。     

Multi-Scale Texture Modeling

Patterns of crystallographic preferred orientation are referred to as texture. The specific subject of texture analysis is the experimental determination and interpretation of the statistical distribution of orientations of crystals within a specimen of polycrystalline material, which could be metals or rocks. The objective is to relate an observed pattern of preferred orientation to its generating processes and vice versa. In geosciences, texture of minerals in rocks is used to infer constraints on their tectono-metamorphic history. Since most physical properties of crystals, such as elastic moduli, the coefficients of thermal expansion, or chemical resistance to etching depends on crystal symmetry and orientation, the presence of texture imparts directional properties to the polycrystalline material. A major issue of mathematical texture analysis is the resolution of the inverse problem to determine a reasonable orientation density function on SO(3) from measured pole intensities on , which relates to the inverse of the totally geodesic Radon transform. This communication introduces a wavelet approach into mathematical texture analysis. Wavelets on the two-dimensional sphere and on the rotational group SO(3) are discussed, and an algorithms for a wavelet decomposition on both domains following the ideas of Ta-Hsin Li is given. The relationship of these wavelets on both domains with respect to the totally geodesic Radon transform is investigated. In particular, it is shown that the Radon transform of these wavelets on SO(3) are again wavelets on . A novel algorithm for the inversion of experimental pole intensities to an orientation density function based on this relationship is developed.     

Actinic flux and photolysis frequency comparison computations using the model PHOTOGT

The accurate radiative transfer model GOMETRAN, initially designed to yield radiances at TOA in the wavelength range 240–790 nm, has been extended to allow for the computation of actinic fluxes down to 175 nm and for the calculation of photolysis frequencies in the atmosphere. The capability of the extended model PHOTOGT (PHOTOGOMETRAN) is demonstrated in a number of successful comparison studies both with recent experiments (ground-based, balloonborne, airborne) and model calculations of radiances, actinic fluxes and photolysis frequencies in the stratosphere and troposphere. In an atmospheric case study, the impact of new quantum yield data for the O₃ » O₂+O(¹ d) photodissociation channel on the photolytic production of O(¹ d) atoms in the lower atmosphere has been quantified.     

Evolution of cool skin and direct air-sea gas transfer coefficient during daytime

Absorption of solar radiation within the thermal molecular sublayer of the ocean can modify the temperature difference across the cool skin as well as the air-sea gas transfer. Our model of renewal type is based on the assumption that the thermal and diffusive molecular sublayers below the ocean surface undergo cyclic growth and destruction, the heat and gas transfer between the successive burst events are performed by molecular diffusion. The model has been upgraded to include heating due to solar radiation. The renewal time is parameterized as a function of the surface Richardson number and the Keulegan number. A Rayleigh number criterion characterizes the convective instability of the cool skin under solar heating. Under low wind speed conditions, the solar heating can damp the convective instability, strongly increasing the renewal time and correspondingly decreasing the interfacial gas exchange. In the ocean, an additional convective instability caused by salinity flux due to evaporation becomes of importance in such cases. The new parameterization is compared with the cool skin data obtained in the western equatorial Pacific during the Tropical Ocean Global Atmosphere Coupled Ocean Atmosphere Response Experiment in February 1993. In combination with a model of the diurnal thermocline it describes main features of the field data both in nighttime and daytime. Under low wind speed conditions (< 5 m s^-1) diurnal variations of the sea surface temperature due to the formation of a diurnal thermocline were substantially larger than those across the cool skin. Under wind speeds > 5 m s^-1, diurnal variations of the surface temperature due to the variations of the thermal molecular sublayer become more important.     

地震综合效应场函数及其在地震预报中的应用──（一）地震综合效应场函数及其计算方法

在系统地分析了目前各种测震学地震预报方法科学思路的基础上,认为测震学地震预报方法基本上可以分为两大类。一类是以已经发生的一些地震作为未来可能发生的地震的“因”,即由于已经发生的地震对区域应力场的影响,导致未来发生较强地震。这一类包括的预报方法较多,如空区、条带、ｂ值、地震迁移、相关地震等等及其由此衍生出来的各种方法。另一类是把已经发生的一些地震作为区域应力场增强的“果”,即已经发生的地震是区域应力场增强过程中的一种反映,而未来地震不一定是已经发生的地震所导致的结果。这一类包括“地震窗口”、小震群活动等方法。针对第一类方法,各种预报方法都是力图从地震三要素中提取未来地震的信息,而具体作法又都是利用地震三要素这个多维空间的某个剖面。为了从地震活动诸要素的多维空间提取综合信息,我们对每个地震加入了破裂面方位,构成了地震第四要素,并依据地震４要素建立了地震综合效应场函数。地震综合效应场函数概括了多种测震学地震预报方法的科学思路和预报经验,从而可以形成测震学的综合预报方法。     

Modelling of the vertical fluxes of nitric acid,ammonia, and ammonium nitrate

Results from numerical investigations regarding the exchange of HNO₃, NH₃, and NH₄NO₃ between the atmosphere and the biosphere are presented. The investigations were performed with a modified inferential method which is based on the generally accepted micrometeorological ideas of the transfer of momentum, sensible heat and matter near the Earth's surface and the chemical reactions among these nitrogen compounds. This modified inferential method calculates the micrometeorological quantities (such as the friction velocity and the fluxes of sensible and latent heat), the height-invariant fluxes of the composed chemically conservative trace species with group concentrationsc ₁=[HNO₃]+[NH₄NO₃] (total nitrate),c ₂=[NH₃]+[NH₄NO₃] (total ammonia), andc ₃=[HNO₃]-[NH₃] as well as the fluxes of the individual nitrogen compounds. The parameterization of the fluxes is based on the flux-gradient relationships in the turbulent region of the atmospheric surface layer. The modified inferential method requires only the data of wind velocity, temperature, humidity and concentrations (HNO₃, NH₃, and NH₄NO₃) measured at a reference height by stations of a monitoring network.     

Estimation of covariance matrix from geochemical data with observations below detection limits

Multivariate statistical analyses have been extensively applied to geochemical measurements to analyze and aid interpretation of the data. Estimation of the covariance matrix of multivariate observations is the first task in multivariate analysis. However, geochemical data for the rare elements, especially Ag, Au, and platinum-group elements, usually contain observations the below detection limits. In particular, Instrumental Neutron Activation Analysis (INAA) for the rare elements produces multilevel and possibly extremely high detection limits depending on the sample weight. Traditionally, in applying multivariate analysis to such incomplete data, the observations below detection limits are first substituted, for example, each observation below the detection limit is replaced by a certain percentage of that limit, and then the standard statistical computer packages or techniques are used to obtain the analysis of the data. If a number of samples with observations below detection limits is small, or the detection limits are relatively near zero, the results may be reasonable and most geological interpretations or conclusions are probably valid. In this paper, a new method is proposed to estimate the covariance matrix from a dataset containing observations below multilevel detection limits by using the marginal maximum likelihood estimation (MMLE) method. For each pair of variables, sayY andZ whose observations containing below detection limits, the proposed method consists of three steps: (i) for each variable separately obtaining the marginal MLE for the means and the variances, , , , and forY andZ: (ii) defining new variables by and and lettingA=C+D andB=C–D, and obtaining MLE for variances, and forA andB; (iii) estimating the correlation coefficient _YZ by and the covariance _YZ by . The procedure is illustrated by using a precious metal geochemical data set from the Fox River Sill, Manitoba, Canada.