中国数字地质调查系统的基本构架 及其核心技术的实现

数字地质调查系统是贯穿整个地质矿产资源调查全过程的软件,涵盖地质矿产调查、矿产资源勘查、矿体模拟、品住估计、资源量估算、矿山开采系统优化等内容。考虑数字地质调查系统的应用技术层面,从数据“层”模型、数据流“池”技术、不同阶段数据模型继承技术、数据互操作技术等几个方面讨论了核心技术及其实现,通过基于无缝一体化技术的数据采集、管理、综合处理与成果表达,在实体或矿块的矿床建模技术与品住估计、储量估算等方面展现了全地质调查过程的数字化,不但为地质人员应用高新技术降低了门坎,而且极大地提高了研究精度和效率,丰富了成果的表现形式和服务形式。随着数字地质调查系统的完善和应用水平的提高,数字地质调查系统将成为中国地质调查的主流软件体系。     

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.     

Oil and gas resource assessment: The linkage between geological analysis and discovery process model results

Two widely-used techniques to estimate the volume of remaining oil and gas resources are discovery process modeling and geologic assessment. Both were used in a recent national assessment of oil and gas resources of the United States. Parallel estimates were obtained for 27 provinces. Geological-based estimates can typically see into areas not available to discovery process models (that is areas with little or no exploration history) and thus, on average, yield higher estimates. However, a linear relation does exist between the mean estimates obtained from these two methods. In addition, other variables were found in a multiple regression model that explained much of the difference. Thus, it is possible to perform discovery process modeling and adjust the estimates to yield results that might be expected from geological-based assessments.     

Earthquake rupture complexity due to dynamic nucleation and interaction of subsidiary faults

We numerically study the dynamic interaction of propagating cracks. It is assumed that propagating cracks can nucleate and drive subsidiary cracks because of shear strain enhancement near the propagating crack tips. The critical strain fracture criterion is assumed in the analysis. Intense interaction is expected to occur among the cracks. All the cracks are assumed to be parallel and antiplane strain deformation is assumed in the computation.In the interaction of two non-coplanar cracks, a strain shadow is formed in the neighborhood of each crack because of the strain release by the introduction of the crack. The growth of each crack is accelerated when the propagating tips of each crack are outside of the strain shadow of the other crack. In general, the crack tips enter the strain shadow, and the crack tips decelerate. The calculation shows that only one of the two cracks can continue to grow, and the other's growth is decelerated and arrested. If we can assume that the suite of cracks interact in a pairwise manner only, then this may suggest that only a limited number of cracks can continue to grow during the final stage of the rupture process. Hence the crack interaction causes complexity in dynamic earthquake faulting. The concepts of barrier and asperity have been employed by many researchers for the interpretation of complex seismic wave data. However, the physical realities of such concepts are obscure. Our calculations show that dynamic crack interactions can produce barriers and asperities in some cases; the crack tip deceleration or arrest due to the interactions among non-coplanar cracks can be interpreted as being due to a barrier. The dynamic coalescence among the coplanar cracks can be regarded as an asperity.Umeda found a localized area that strongly radiates high-frequency seismic waves in the epicentral areas of some large shallow earthquakes. He defined this as an earthquake bright spot. Our analysis implies that only a limited number of cracks continue to grow when many interactive cracks nucleate, and that all other cracks stop extending soon after nucleation. Hence, if the nucleation and termination of several cracks occur in a localized area, it will be observed seismologically as an earthquake bright spot. This is because it is theoretically known that the sudden termination of crack growth and dynamic crack coalescence efficiently emits high-frequency elastic waves.     

Interpolation by regularized spline with tension: I. Theory and implementation

Bivariate and trivariate functions for interpolation from scattered data are derived. They are constructed by explicit minimization of a general smoothness functional, and they include a tension parameter that controls the character of the interpolation function (e.g., for bivariate case the surface can be tuned from a membrane to a thin steel plate), Tension can be applied also in a chosen direction, for modeling of phenomena with a simple type of anisotropy. The functions have regular derivatives of all orders everywhere. This makes them suitable for analysis of surface geometry and for direct application in models where derivatives are necessary. For processing of large datasets (thousands of data points), which are now common in geosciences, a segmentation algorithm with a flexible size of overlapping neighborhood is presented. Simple examples demonstrating flexibility and accuracy of the functions are presented.On leave from the Department of Physical Geography and Cartography, Comenius University, Mlynská dolina, Bratislava, Czechoslovakia.On leave from the Institute of Physics, Dúbravská cesta 9, Bratislava, Czechoslovakia.     

胶东海阳所超镁铁岩中强亲铁元素的地球化学特征及其构造环境探讨

强亲铁元素与亲石元素具有不同的地球化学行为，因此能够从不同的角度为造山带中超镁铁岩的成因及演化提供重要信息。位于苏鲁造山带东北端的胶东海阳所超镁铁岩主要由橄榄岩和辉石岩组成，它们常以团块状赋存于花岗质片麻岩中。虽然前人对这些超镁铁岩已经开展大量岩石学研究，但关于其成因及构造属性仍存在较大争议。本文开展了海阳所超镁铁岩的全岩主微量元素、强亲铁元素及Re-Os同位素的分析工作，结果显示蛇纹石化橄榄岩具有较高的MgO和Fe₂O₃^T含量，较低的Al₂O₃、TiO₂和CaO含量，明显富集流体迁移元素（U、Pb），亏损高场强元素（Zr、Hf），强亲铁元素没有发生明显分异，但Ru显示正异常，表明海阳所蛇纹石化橄榄岩是经历了低-中等程度部分熔融及熔/流体交代作用影响的残余地幔橄榄岩。海阳所辉石岩的主量元素表现出明显的结晶分异特征，稀土元素较原始地幔富集，铂族元素（PGEs）含量较低且发生了明显的分异，表明辉石岩的地幔源区经历过高程度的部分熔融和硫化物的分离。海阳所蛇纹石化橄榄岩的Os同位素地球化学特征表现出大洋亲和性，与辉石岩不具有熔体-残留体的关系。由于该地区发育较深层次的韧性剪切带，蛇纹石化橄榄岩中的橄榄石与辉石表现出韧性变形的特征，同时有辉石岩侵入到橄榄岩的现象，表明该地区的蛇纹石化地幔橄榄岩与辉石岩既不同时，也不同源，因此，暗示了该套岩石组合可能形成于大洋核杂岩（OCC）与洋脊型蛇绿岩（MOR）堆晶岩交互发育环境。     

内蒙古中西部宝音图群石榴角闪岩的地球化学、变质作用及年代学特征

宝音图群分布于内蒙古中西部狼山、图古日格和达茂旗一带，主要由云母片岩、石英片岩、石英岩、角闪岩和大理岩组成。在狼山格尔敖包沟和图古日格西南地区发育十分典型的石榴角闪岩，均以似层状或透镜状产自云母片岩中，但两地具有不同的矿物组合，前者主要以石榴石+角闪石+斜长石+含钛矿物（金红石和钛铁矿）为特征，后者主要以石榴石+角闪石+绿帘石+含钛矿物（金红石和钛铁矿）为特征。地球化学研究显示两个地区的石榴角闪岩原岩具有相似的化学组成，均为亚碱性玄武岩的拉斑系列。宝音图群的石榴角闪岩具有右倾的稀土配分模式，REE总量为83.31×10^-6~125.9×10^-6，（La/Yb）_N比值为2.17~6.48，δEu=0.87~0.98。Ta、Nb、Ti没有明显的负异常，其配分型式类似E-MORB特点。微量元素构造判别图解表明这些石榴角闪岩的原岩产于板内拉张环境中。本文通过模拟格尔敖包沟样品LS01和图古日格西南地区样品LS35的P-T视剖面图研究其变质作用及矿物演化过程，结果显示这两个样品均经历了早期进变质，峰期以及峰后近等温减压的顺时针型P-T轨迹。样品LS01的峰期温压条件为~11kbar/~735℃，模拟得到的峰期矿物组合可能为石榴石+角闪石+斜长石+黑云母+透辉石+金红石+熔体，P-M（H₂O）视剖面图显示，由于峰后演化过程中存在饱和流体渗透，导致峰期透辉石无法保留。样品LS35的峰期温度压力条件为~8kbar/675℃，其峰期组合为石榴石+角闪石+绿帘石+金红石。两个地区的石榴角闪岩的温度压力条件明显不同，推测宝音图群中呈透镜体或似层状产出的变质基性岩，存在着递增型的变质作用。锆石LA-ICP-MS U-Pb年代学的研究结果显示石榴角闪岩样品LS01的变质锆石年龄为394±8Ma。由此推测这些石榴角闪岩的原岩可能形成于华北北缘在新元古代发育的裂陷盆地，在泥盆纪中期，宝音图群作为华北克拉通的一部分，被卷入到与古亚洲洋闭合有关的造山过程，发育中压型递增变质作用。     

Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

Seasonal variation of atmospheric dimethylsulfide at Amsterdam Island in the southern Indian Ocean

Daily measurements of atmospheric concentrations of dimethylsulfide (DMS) were carried out for two years in a marine site at remote area: the Amsterdam Island (37°50S–77°31E) located in the southern Indian Ocean. DMS concentrations were also measured in seawater. A seasonal variation is observed for both DMS in the atmosphere and in the sea-surface. The monthly averages of DMS concentrations in the surface coastal seawater and in the atmosphere ranged, respectively, from 0.3 to 2.0 nmol l^-1 and from 1.4 to 11.3 nmol m^-3 (34 to 274 pptv), with the highest values in summer. The monthly variation of sea-to-air flux of DMS from the southern Indian Ocean ranges from 0.7 to 4.4 mol m^-2 d^-1. A factor of 2.3 is observed between summer and winter with mean DMS fluxes of 3.0 and 1.3 mol m^-2 d^-1, respectively.