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.     

用褶积滤波处理井水位对固体潮响应的滞后

本文给出了一种用褶积滤波处理井水位对固体潮响应滞后的新方法。鲁07井水位观测资料用这种方法处理比一般回归方法的平均中误差少4.6mm。文中讨论了褶积滤波积分区间长度与水井含水层渗透系数的关系及褶积滤波比别尔采夫滤波的优越性。     

The almost constant-speed two-body problem with resistance

The almost constant-speed motion of a mass acted upon by a Newtonian attraction and a resisting force is treated. The equation of orbit is derived for a specific type of resistance which covers the familiar case of Danby's drag(=r ^–2) whilst the vector invariants are obtained by direct operation on the vector form of the equation of motion.     

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

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

单频GPS接收机天线扼流圈的研制与测试

讨论了多路径效应对GPS共视比对的影响及削弱其影响的方法。结合综合原子时项目的需求,研制了适用于单频GPS接收机的天线扼流圈,并用零基线比对方法进行了测试。结果表明,研制的天线扼流圈可有效地削弱多路径效应的影响,并可明显地改善单频GPS接收机NTSCGPS-1的抗干扰能力。     

Mais comment s'écoule donc un glacier ? Aperçu historique

Ice and snow have often helped physicists understand the world. On the contrary it has taken them a very long time to understand the flow of the glaciers. Naturalists only began to take an interest in glaciers at the beginning of the 19th century during the last phase of glacier advances. When the glacier flow from the upslope direction became obvious, it was then necessary to understand how it flowed. It was only in 1840, the year of the Antarctica ice sheet discovery by Dumont d'Urville, that two books laid the basis for the future field of glaciology: one by Agassiz on the ice age and glaciers, the other one by canon Rendu on glacier theory. During the 19th century, ice flow theories, adopted by most of the leading scientists, were based on melting/refreezing processes. Even though the word ‘fluid’ was first used in 1773 to describe ice, more the 130 years would have to go by before the laws of fluid mechanics were applied to ice. Even now, the parameter of Glen's law, which is used by glaciologists to model ice deformation, can take a very wide range of values, so that no unique ice flow law has yet been defined. To cite this article: F. Rémy, L. Testut, C. R. Geoscience 338 (2006).     

三江平原旱田耕作对湿地土壤理化性质的累积影响初探

在中国科学院三江平原湿地生态试验站综合实验场,对不同耕作年限的湿地土壤(0~20 cm)进行环境累积效应分析。结果表明,随着开垦时间的增加,土壤的理化性质发生渐变,物理性质方面,土壤容重和比重逐渐增大,而孔隙度和田间持水量逐年减少;化学性质方面,土壤pH值随开垦时间的增加而增加,有机质和其他养分则随开垦时间的增加而逐年降低。弃耕后土壤性质有所恢复。土壤性质在开垦初期变化较明显,而后逐渐变缓。