地震勘探测量数据的快速处理和管理

目前山区地震测量施工主要的困难体现在山区放样难、内业数据处理复杂、测量成果无法进行异常检查等方面,其严重影响后续资料处理和解释质量。针对以上难点,利用AUTOCAD为平台,开发了一套物探测量数据处理系统。该系统具有展点、点内插、绘制测线剖面图、管理测点信息、打印成果报表等功能,实现了地震勘探测量数据的快速处理及有效检核．提高了野外测量放样的施工效率和数据的准确度。另外其图层管理及系统辅助功能,还可为地震勘探施工设计提供参考。     

Non-Linear Theory and Power-Law Models for Information Integration and Mineral Resources Quantitative Assessments

Singular physical or chemical processes may result in anomalous amounts of energy release or mass accumulation that, generally, are confined to narrow intervals in space or time. Singularity is a property of different types of non-linear natural processes including cloud formation, rainfall, hurricanes, flooding, landslides, earthquakes, wildfires, and mineralization. The end products of these non-linear processes can be modeled as fractals or multifractals. Hydrothermal processes in the Earth’s crust can result in ore deposits characterized by high concentrations of metals with fractal or multifractal properties. Here we show that the non-linear properties of the end products of singular mineralization processes can be applied for prediction of undiscovered mineral deposits and for quantitative mineral resource assessment, whether for mineral exploration or for regional, national and global planning for mineral resource utilization. In addition to the general theory and framework for the non-linear mineral resources assessment, this paper focuses on several power-law models proposed for characterizing non-linear properties of mineralization and for geoinformation extraction and integration. The theories, methods, and computer system discussed in this paper were validated using a case study dealing with hydrothermal Au mineral potential in southern Nova Scotia, Canada.     

(ϖμ/ϖT) P of the lower mantle

We estimate (/T) _P of the lower mantle at seismic frequencies using two distinct approaches by combining ambient laboratory measurements on lower mantle minerals with seismic data. In the first approach, an upper bound is estimated for |(/T) _P | by comparing the shear modulus () profile of PREM with laboratory room-temperature data of extrapolated to high pressures. The second approach employs a seismic tomography constraint ( lnV _S / lnV _P ) _P =1.8–2, which directly relates (/T) _P with (K _S /T) _P . An average (K _S /T) _P can be obtained by comparing the well-established room-temperature compression data for lower mantle minerals with theK _S profile of PREM along several possible adiabats. Both (K _S /T) and (/T) depend on silicon content [or (Mg+Fe)/Sil of the model. For various compositions, the two approaches predict rather distinct (/T) _P vs. (K _S /T) _P curves, which intersect at a composition similar to pyrolite with (/T) _P =–0.02 to –0.035 and (K _S /T) _P =–0.015 to –0.020 GPa/K. The pure perovskite model, on the other hand, yields grossly inconsistent results using the two approaches. We conclude that both vertical and lateral variations in seismic velocities are consistent with variation due to pressure, temperature, and phase transformations of a uniform composition. Additional physical properties of a pyrolite lower mantle are further predicted. Lateral temperature variations are predicted to be about 100–250 K, and the ratio of ( lnp/ lnV _S ) _P around 0.13 and 0.26. All of these parameters increase slightly with depth if the ratio of ( lnV _S / lnV _P ) _P remains constant throughout the lower mantle. These predicted values are in excellent agreement with geodynamic analyses, in which the ratios ( ln / lnV _S ) _P and ( / lnV _S ) _P are free parameters arbitrarily adjusted to fit the tomography and geoid data.     

Hierarchical decomposition of variance with applications in environmental mapping based on satellite images

A quadtree-based image segmentation procedure (HQ) is presented to map complex environmental conditions. It applies a hierarchical nested analysis of variance within the framework of multiresolution wavelet approximation. The procedure leads to an optimal solution for determining mapping units based on spatial variability with constraints on the arrangement and shape of the units. Linkages to geostatisiics are pointed out, but the HQ decomposition algorithm does not require any homogeneity criteria. The computer implementation can be parameterized by either the number of required mapping units or the maximum within-unit variance, or it can provide a spectrum of significances of nested ANOVA. The detailed mathematical background and methodology is illustrated by a salt-affected grassland mapping study (Hortobágy, Hungary), where heterogeneous environmental characteristics have been sampled and predicted based on remotely sensed images using these principles.     

The “Zero-order Model” revisited: Paul Schindler's influence on the development of trace metal scavenging models

In 1975 Paul Schindler produced the first oceanic trace metal scavenging model to explicitly include the role of surface chemistry as a control on trace metal water column residence times. The eighteen years that have elapsed since the publication of Schindler's seminal paper have seen the development of a variety of oceanic scavenging models; yet, the fundamental insight of his Zero-order Model remains the benchmark. This paper describes the role of Paul Schindler's work on surface chemistry in providing a framework for the current generation of trace element scavenging models.     

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.     

The development of radio astronomy at Hartebeesthoek

The development of radio astronomy at the Hartebeesthoek Radio Astronomy Observatory in South Africa is described. The Hartebeesthoek site was established originally by NASA as one of three Deep Space Stations equipped with 26-m parabolic reflector antennas. It was first used for radio astronomy by South Africa in terms of the NASA host nation agreement which allowed for its use at times when the facility was not needed for its primary purpose of tracking space probes. After NASA withdrew from South Africa in 1975, the South African Council for Scientific and Industrial Research took over the site and the 26-m parabolic reflector antenna, which NASA had abandoned in position, and established it as a national observatory. The development of the facility to the stage where it could support a variety of observing programmes such as continuum observations and mapping, spectroscopy and pulsar timing is described as well as the role played by the Observatory in global programmes of very long baseline interferometry.     

Using numerical modelling to evaluate tsunami hazard near the Kuril Island

A sequence of computer experiments is used to study questions concerning the tsunami problem as a quantitative estimate of tsunami danger, detailed geographical tsunami classification, determination of the parameters of critical tsunami waves, and the conditions of their development. We call a wave critical, if its impact on the coast is most hazardous.Using the Middle Kuril Island as an example, we present the results of a computer experiment which includes determining the wavefields on the shelf and estimating the effects connected with the deep-water Bussol and Diana Straits.Numerical simulation of tsunami waves of different sources permits the assessment of the extent of tsunami danger in different areas of the coastal zone of Simushir Island, depending on the location of the focus zone and their geometry.The major singularities of the wavefield arise in the zones of the deep-water straits. The distribution of the amplification factors is determined by both the global parameters of the wavefields and the local properties of individual harbours. The results obtained for a particular harbour in the northern part of Simushir Island, formed the basis for the quantitative estimate of tsunami danger for this area.     

Geochemical mapping using a geomorphologic approach based on catchments

The Extended Sample Catchment Basin (ESCB) mapping technique, discussed in this paper, can be used to display the spatial distribution of geochemical variables measured in stream sediments taking into consideration the geomorphologic settings and the hydrographic patterns of surveyed areas. This approach is based on the association of an area of statistical representativeness with each sample, and on the assumption that the concentrations measured in the stream sediments can be considered as average reference values for this area.ESCBs can be easily identified considering the position of the sampling points within the hydrographic network and using the confluences between the streams of highest rank as break points for representing changes of the geochemical background. This approach, different from the traditional geostatistical and deterministic ones, does not consider the Euclidean distance among the sampling points as a measure of geochemical similarity but only refers to their functional relationship along the streams (following the water and the sediment flow) to measure their proximity.ESCBs can be seen as a specific development of previous techniques based on catchments and proves to be especially useful for supporting land planning in a preliminary survey phase while it is not specifically suitable for the identification of point sources of geochemical anomalies.Due to the fractal nature of the hydrographic network, all the procedures can be driven in a GIS environment by using digital terrain models apart from their spatial resolution.     

塔里木河源区冰川系统变化趋势预测

塔里木河源区是我国冰川分布最集中的地区之一,总面积达17 745.51 km2,占全国冰川总面积的30%;同时本区又属于我国升温幅度最大的地区之一。应用冰川系统变化的功能模型,对塔里木河源区冰川系统在本世纪对气候变化的趋势进行预测。结果表明:到2050年,如气温比1961~1990年高出1.9~2.3℃,本区冰川面积将减少4%~6%,冰川径流将增加22%~34%,零平衡线将上升62~94 m;如此升温率持续到本世纪末,则本区冰川面积将减少10%~16%,冰川径流将会回落,但仍比本世纪初多11%~13%,零平衡线将上升156~233 m。