Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.     

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.     

鲁23井“动水位”观测的分析

本文介绍了鲁23井“动水位”观测情况。提出了一个分析有泄流的“动水位”的新方法,即把流量换算成水位,以等效静水位来表示“动水位”与流量。由于把“动水位”与流量用一个量来计算,从而为“动水位”观测定量的、统一的分析研究提供了新的方法。这对提高“动水位”井的监测预报地震能力起一定的作用。 本文还以鲁23井“动水位”观测资料为例,利用新方法进行处理的观测分析结果与理论研究结果是一致的。     

Paleohydrology and paleoclimate records in a saline prairie lake core: mineral,isotope and organic indicators

This pilot study examines the potential of obtaining a sedimentary record of paleoenvironmental/climatic/hydrologic conditions for saline Redberry Lake in southern Saskatchewan, Canada. The tools are mineralogy, stable isotopes and pigments. The upper meter of an offshore sediment core contains 10 to 20% by weight aragonite (CaCO₃), which apparently precipitated in the water column. The ¹⁸O and ¹³C of the bulk aragonite (corrected for content of detrital calcite) vary by 4 to 5. Enrichment in ¹⁸O in aragonite is significantly correlated with pigment concentrations (chlorophyll a, phaeophytin). The ¹⁸O and pigment data provide evidence for relatively dry and/or warm conditions and high limnetic productivity for the period 2500 to 1500 yrs B.P. After 1500 B.P., the climate was apparently similar to the present, with two episodes of relatively enhanced productivity, dryness and/or warmth, at around 1000 to 900 and 500 to 200 B.P. During the past century, Redberry Lake has decreased approximately 8 m in depth and its salinity has doubled. No clear sedimentary signal was observed in response to these recent hydrologic trends. These changes have not been associated with a significant climate trend in the region, but may have been induced by land use changes in the catchment.This publication is the third of a series of papers presented at the Conference on Sedimentary and Paleolimnological Records of Saline Lakes. This Conference was held August 13–16, 1991 at the University of Saskatchewan, Saskatoon, Canada. Dr. Evans is serving as Guest Editor for this series.     

珠穆朗玛峰高程与板块构造运动问题

喜马拉雅山是由印度板块和欧亚板块相碰撞而形成,这已被当代多数有关科学家所认识。一些测量数据也表明喜马拉雅山现在仍在快速隆升。而喜马拉雅山主峰——珠穆朗玛峰的高程的测定,也成为当代地学界的热点。本文结合最新资料对精确测定珠峰高程的一些问题,诸如水准基面、峰顶标志和积雪等进行了讨论。     

Sedimentary Features of Shallow Ancient River Channels on the Northern Shelf of the South China Sea

Quaternary buried ancient river channels are widespread in the shallow-level sediments of the northern shelf of the South China Sea. The sedimentary sequence mainly of fluvial deposits comprise an important component part of the low-stand system tract and transgressive system tract in the study region. The plannar variation and spatial association of the sedimentary features such as incised valley fillings, deltaic foreset wedges and block slides of shelf-marginal fans reflect the palaeogeographic environment during the fall of the regional sea level in the northern part of the South China Sea. Based on the high-resolution seismic reflection data and gelogical data from boreholes, the present paper makes an integrated interpretation of the Quaternary ancient river channels in the shallow sediments of the study area, studies the sedimentary features of the ancient channels such as their spatial distribution, seismic facies reflection indicators, sedimentary facies and sand -body types, and discusses thei     

使用价值分析优化方法在功能—技术矩阵求解中的应用

功能—技术矩阵,经众多工作者的努力研究,已渐趋成熟与实用。但如何利用功能—技术矩阵构造在广义现代设计领域中适用的科学模式,解决多学科设计方案的组成与评价问题,还存在薄弱环节。本文利用“使用价值分析优化”的概念,从设计方案技术先进性和经济合理性的协调统一的角度出发,较好地实现了功能—技术矩阵可行优化组合链的组成与评价,弥补了以往的缺陷。     

城乡一体化渐进式开发及发展模式研究--以成都市为例

城乡一体化成为当前研究的热点问题之一。成都市从2003年开始，在全市范围内开展城乡一体化试点，准确把握城乡一体化的发展模式，是顺利推进城乡一体化发展的重要一环。根据成都市城乡一体化发展的特点，采用回归方法预测了城市化发展水平，并借鉴国际城市和国内沿海城市的发展经验，提出成都市城乡一体化应采取的总体思路和发展模式，结合成都市城市化水平预测结果，提出其分阶段、分步骤实施的时间阶段。     

Testing earthquake prediction methods: «The West Pacific short-term forecast of earthquakes with magnitude MwHRV ≥ 5.8»

Analyzing the tables and probability maps posted by Yan Y. Kagan and David D. Jackson in April 2002–September 2004 at http://scec.ess.ucla.edu/~ykagan/predictions_index.html and the catalog of earthquakes for the same period, the conclusion is drawn that the underlying method could be used for prediction of aftershocks, while it does not outscore random guessing when main shocks are considered.