可用于带耗散力卫星轨道长间隔积分的数值积分方法

本文分析了对称方法不适合带耗散力的卫星轨道长间隔积分的缺陷和本质原因,并针对这个问题,采用文所介绍的方法,构造并推荐了一组积分公式。     

利用ARC/INFO进行三维地形模型的建立、处理和动态模拟

介绍了利用 ARC/ INFO建立三维地形模型的方法 ,以及利用其新版本的部分新功能处理三维地形模型、进行动态模拟的方法 ,与以往处理三维模型必须借助外部三维图形开发软件包的方法相比 ,具有高速、灵活、可视性强、支持强大分析功能等优点。形成的三维模型可随时与 ARC/ INFO的其它数据模型交换信息 ,与外部进程实时通讯 ,为使用者提供了二维分析工具所不能提供的信息 ,为开发视景仿真系统等与表面分析、立体显示有关的应用系统提供了强有力的支持。     

矿床统计预测单元划分的方法与程序

在矿床统计预测中,如何确定最佳的统计预测单元,目前尚无通用的准则和算法,它的划分往往取决于地质学家对研究区控矿条件的认识,矿化的实际分布殂态及统计预测所采用的评价模型等多种因素。系统功能齐全,基本上满足了资源评价的各种需求,实现了预测单元划分的自动化和智能化。     

An Improved Perturbation Mechanism for Simulated Annealing Simulation

Simulated annealing (SA) is being increasingly used for the generation of stochastic models of spatial phenomena because of its flexibility to integrate data of diverse types and scales. The major shortcoming of SA is the extensive CPU requirements. We present a perturbation mechanism that significantly improves the CPU speed. Two conventional perturbation mechanisms are to (1) randomly select two locations and swap their attribute values, or (2) visit a randomly selected location and draw a new value from the global histogram. The proposed perturbation mechanism is a modification of option 2: each candidate value is drawn from a local conditional distribution built with a template of kriging weights rather than from the global distribution. This results in accepting more perturbations and in perturbations that improve the variogram reproduction for short scale lags. We document the new method, the increased convergence speed, and the improved variogram reproduction. Implementation details of the method such as the size of the local neighborhood are considered.     

Accounting for Estimation Optimality Criteria in Simulated Annealing

This paper presents both estimation and simulation as optimization problems that differ in the optimization criteria, minimization of a local expected loss for estimation and reproduction of global statistics (semivariogram, histogram) for simulation. An intermediate approach is proposed whereby an initial random image is gradually modified using simulated annealing so as to better match both local and global constraints. The relative weights of the different constraints in the objective function allow the user to strike a balance between smoothness of the estimated map and reproduction of spatial variability by simulated maps. The procedure is illustrated using a synthetic dataset. The proposed approach is shown to enhance the influence of observations on neighboring simulated values, hence the final realizations appear to be better conditioned to the sample information. It also produces maps that are more accurate (smaller prediction error) than stochastic simulation ignoring local constraints, but not as accurate as E-type estimation. Flow simulation results show that accounting for local constraints yields, on average, smaller errors in production forecast than a smooth estimated map or a simulated map that reproduces only the histogram and semivariogram. The approach thus reduces the risk associated with the use of a single realization for forecasting and planning.     

Flow Simulations To Evaluate Upscaling of Permeability

We study upscaling of the permeability for porous media flow on a grid with one million blocks. The purpose is to illustrate how flow simulations can be used to evaluate upscaling methods.     

Closed Form Solutions of the Two-Dimensional Turning Bands Equation

The turning bands method generates realizations of isotropic Gaussian random fields by means of appropriately summed line processes. For two-dimensional simulations the relation between the isotropic correlation function of the random field and the correlation function to be simulated along the lines is given by an integral equation of Abel type. We present closed form solutions of this integral equation for almost all two-dimensional correlation models encountered in practice and discuss their numerical implementation. As an additional benefit, our tables and illustrations serve as a concise guide to correlation models useful in geostatistics.     

Laboratory sulfurisation of the marine microalga Nannochloropsis salina

To understand more fully the mode of preservation of organic matter in marine sediments, laboratory sulfurisation of intact cells of the cultured microalga Nannochloropsis salina was performed using inorganic polysulfides in seawater at 50°C. Solvent extractable and non-extractable materials were analysed by GC–MS and Py–GC–MS, respectively, to study the incorporation of sulfur into the microalgal organic matter. No GC-amenable sulfur-containing compounds were found in the extracts apart from some minor thiophenes with a phytanyl carbon skeleton. The residue after extraction and hydrolysis contained abundant macromolecular sulfur-containing moieties as revealed by the presence of dominant C₂₈–C₃₂ thiols, thiophenes, thianes and thiolanes in the flash pyrolysates. These products are thought to be formed from moieties derived from sulfurisation of C₂₈–C₃₂ diols and alkenols, characteristic lipids of N. salina. C₁–C₂ alkylated thiophenes were also found in the pyrolysates and probably result from moieties formed upon sulfurisation of carbohydrates. The highly resistant biomacromolecule (algaenan) synthesised by N. salina remains unaffected by sulfurisation. The non-hydrolysable residue isolated from the sulfurised N. salina thus comprises algaenan and (poly)sulfide-bound long alkyl chains. The sulfurisation experiments show that both selective preservation of algaenans and lipid and carbohydrate “vulcanisation” can be involved in the preservation of algal organic matter in marine environments.     

Some applications of humidity profiles estimated from INSAT infra red digital cloud imagery data

Moisture profiles have been estimated over the region bounded by the latitudes 40°N and 40°S and longitudes 30°E to 130°E using INSAT digital infra red cloud imagery data. The representativeness of these profiles in representing moisture field associated with the development and movement of synoptic scale systems during the period September 15th, 1996 to March 31st, 1997 has been examined. It has been shown that the changes in the moisture field associated with the withdrawal of the southwest and northeast monsoons from the Indian sub-continent, development and movement of synoptic scale sytems (depressions, tropical cyclones and waves in easterlies) and equatorial troughs in the Indian Ocean could be clearly seen in humidity profiles. The initial development of tropical systems is first seen in the humidity field in the upper troposphere. These profiles could be used in monitoring the initial development and subsequent movement of tropical systems. Further the data on moisture distribution from the data gap regions of the Indian Ocean could be used as an additional source of moisture in numerical analysis and prognosis.     

A comparison of estimated and calculated effective porosity

 Effective porosity in solute-transport analyses is usually estimated rather than calculated from tracer tests in the field or laboratory. Calculated values of effective porosity in the laboratory on three different textured samples were compared to estimates derived from particle-size distributions and soil–water characteristic curves. The agreement was poor and it seems that no clear relationships exist between effective porosity calculated from laboratory tracer tests and effective porosity estimated from particle-size distributions and soil–water characteristic curves. A field tracer test in a sand-and-gravel aquifer produced a calculated effective porosity of approximately 0.17. By comparison, estimates of effective porosity from textural data, moisture retention, and published values were approximately 50–90% greater than the field calibrated value. Thus, estimation of effective porosity for chemical transport is highly dependent on the chosen transport model and is best obtained by laboratory or field tracer tests. Received, March 1997 · Revised, August 1997 · Accepted, August 1997