Estimation of covariance matrix from geochemical data with observations below detection limits     

Chang-Jo F. Chung 《Mathematical Geology》1993,25(7):851-865

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.  相似文献   

中国南北带地壳上地幔三维面波速度结构和各向异性   总被引：6，自引：1，他引：5       下载免费PDF全文

陈立华  庄真 《地球物理学报》1992,35(5):574-583

本文用长周期763地震仪面波群速度资料反演了中国南北带及邻区的三维速度结构.其中采集238条瑞利波和358条勒夫波混合频散曲线,使用均等显示滤波方法,并以4°×4°为一格将我国境内分为147格.用随机逆反演方法得到了研究区16格的纯路径频散.面波速度结构及演结果表明:1.莫霍界面深度一般在40-50km之间,最深达65km.总趋势是从东到西加深,且在南北带西侧南北两端向中部明显加深,东侧变化小.2.地幔顶部普遍出现很厚的低速层,上界面一般埋深60-80km.上地幔顶盖厚度一般为20-60km,速度为4.30-4.50km/s.3.研究区普遍存在各向异性,而且勒夫波和瑞利波速度的差值(V_SV-V_SH)的绝对值随深度有增大的特点,在南北带南部和西北部V_SV-V_SH各向异性现象更为明显.  相似文献   

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

罗庆生  韩宝玲 《武汉大学学报(信息科学版)》1990,(2)

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

The Tornquist Zone, a north east inclining lithospheric transition at the south western margin of the Baltic Shield: Revealed through a nonlinear teleseismic tomographic inversion     

P. Voss  K. Mosegaard  S. Gregersen  TOR Working Group 《Tectonophysics》2006,416(1-4):151

From July 1996 to August 1997 the TOR project operated 130 seismographs in North Germany, Denmark and South Sweden, with the aim of collecting signals from local, regional and teleseismic earthquakes. This data set is particularly interesting since the seismic antenna crosses the most significant geological boundary in Europe, the Tornquist Zone, which in the northern part is the border between the Baltic Shield and the younger European lithosphere. Previous studies have shown significant physical changes in the crust and upper mantle across this transition zone, including two independent teleseismic tomographic studies of the TOR data set. But these two studies disagree on the orientation of the slope of the transition. Both studies used an iterative linearized inversion method. We will in this work Preprint submitted to Elsevier Science 27 July 2005 present an inversion based on Bayesian statistics, where the solution space is examined in order to study a very large number of tomographic solutions and to examine the solution uniqueness and uncertainty. The method is applied to measurements of 3345 relative teleseismic P-phase travel times from 48 teleseismic earthquakes with good azimuthal coverage with respect to the great circle arc of the TOR array. We find the lithospheric transition to be a north east inclination of around 30° to 45° off vertical.  相似文献   

Time domain Gauss–Newton seismic waveform inversion in elastic media     

Dong-Hoon Sheen  Kagan Tuncay  Chang-Eob Baag  Peter J. Ortoleva 《Geophysical Journal International》2006,167(3):1373-1384

  相似文献   

New results on the resistivity structure of Merapi Volcano (Indonesia), derived from three-dimensional restricted inversion of long-offset transient electromagnetic data     

Michael Commer  Stefan L. Helwig  reas Hördt  Carsten Scholl  Bülent Tezkan 《Geophysical Journal International》2006,167(3):1172-1187

  相似文献   

Fault plane orientations of small earthquakes of the 1997 Umbria-Marche (Italy) seismic sequence from P-wave polarities and rise times     

Marilena Filippucci  Salvatore de Lorenzo  Enzo Boschi 《Geophysical Journal International》2006,166(1):322-338

  相似文献   

Strain partitioning and stress rotation at the North Anatolian fault zone from aftershock focal mechanisms of the 1999 Izmit M_w= 7.4 earthquake     

Marco Bohnhoff  Helmut Grosser  Georg Dresen 《Geophysical Journal International》2006,166(1):373-385

  相似文献   

Fault slip and rupture velocity inversion by isochrone backprojection     

Gaetano Festa  Aldo Zollo 《Geophysical Journal International》2006,166(2):745-756

  相似文献   

3D behaviour of contamination in landfill sites using 2D resistivity/IP imaging: case studies in Portugal   总被引：2，自引：0，他引：2  

E. Martinho  F. Almeida 《Environmental Geology》2006,49(7):1071-1078

Aiming at defining a valid spatial contamination model, resistivity and induced polarization (IP) measurements were used to investigate contamination plumes in the vicinity of two municipal landfills (Ovar and ílhavo). Previous geophysical surveys and underground water samples confirmed the contamination. However 2D resistivity/IP surveys enabled in obtaining a more accurate spatial model. The Ovar survey consisted of two profiles with nine Wenner soundings each; the ílhavo survey was carried out along two individual lines using a Wenner standard pseudo-section. In both situations, negative IP values were found associated with positive IP values, which can be explained mainly by 2D or 3D geometric effects caused by the presence of the conductive plumes. The data were modelled using a 2D inversion program (RES2DINV) and the resulting resistivity and chargeability distributions were displayed as pseudo-sections. The resistivity and chargeability pseudo-sections define the contamination plumes and the sedimentary structure. These case studies illustrate the advantages of 2D resistivity/IP surveys for the mapping of shape and dimension of contamination associated with landfills.  相似文献