Estimation of statistical parameters for censored lognormal hydraulic conductivity measurements

Field tests of hydraulic conductivity (e.g., injection test, pumping test, etc.) in low permeability formations are subject to censoring due to the detection limit of the instruments used. An iterative method of estimating the mean and variance of hydraulic conductivity data with a presumed log-normal distribution function is presented. This method accounts for the data that are actually below the lower detection limit (called truncated data) and thus gives distribution parameters that are more representative for the underlying distribution. The proposed method is then tested on two simulated normally distributed random datasets having different variances. The results show that the means and variances estimated by the proposed method are very accurate. Finally, the method is used to estimate the mean and variance of hydraulic conductivity data from single hole water injection tests in a fractured geological formation.     

An integrable case of a rotational motion analogous to that of Lagrange and Poisson for a gyrostat in a Newtonian force field

The scope of the present paper is to provide analytic solutions to the problem of the attitude evolution of a symmetric gyrostat about a fixed point in a central Newtonian force field when the potential function isV ⁽²⁾.We assume that the center of mass and the gyrostatic moment are on the axis of symmetry and that the initial conditions are the following: (t ₀)=₀, (t ₀)=₀, (t ₀)=(t ₀)=₀, ₁(t ₀)=0, ₂(t ₀)=0 and ₃(t ₀)= ₃ ⁰ .The problem is integrated when the third component of the total angular momentum is different from zero (B ₁ 0). There now appear equilibrium solutions that did not exist in the caseB ₁=0, which can be determined in function of the value ofl ₃ ^r (the third component of the gyrostatic momentum).The possible types of solutions (elliptic, trigonometric, stationary) depend upon the nature of the roots of the functiong(u). The solutions for Euler angles are given in terms of functions of the timet. If we cancel the third component of the gyrostatic momentum (l ₃ ^r =0), the obtained solutions are valid for rigid bodies.     

矿产预测中成矿信息的提取和强化

成矿信息的提取是矿产预测的重要手段,其方法的有效性直接影响到预测成果的可靠性。在成矿规律研究的基础上,有意识地干预模型的构成,突出与成矿有关的信息;抑制某些成矿意义不明显和属干扰的信息,有目的地使模型向反映成矿信息的方向逼近,提高模型与矿床实际赋存地质环境的吻合程度,均是成矿信息提取和强化的内涵。本文从研究实践中总结了先验约束模型和非先验约束模型的强化方法,在新疆阿勒泰地区的地质－找矿工作中已取得     

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.     

极性转换期间地球磁场形态学研究

本文对采自中国黄土高原西峰（35.7°N,107.6°E）和段家坡（34.2°N,109.2°E）两个剖面中黄土层L8和古土壤层S8的1281块定向古地磁样品做了详细的岩石磁学和古地磁学研究.证实了Matuyama-Brunhes（M-B）极性转换带位于L8的中下部.提出了下列观点:1.M-B极性转换过程与地球磁场方向变化相联系的持续时间为3600-4500a,而与地球磁场强度变化相联系的持续时间则为8000-9000a,即强度变化存在“超前和滞后” 效应;2.M-B转换场的形态是由三次快速倒转和一次不成功的倒转构成,或者说,转换场具有快速变换极性的振荡特征;3.M-B转换过程中地球磁场并不是以轴对称的非偶极子场为主,而是偶极子场至少与非偶极子场相当;4.中国黄土-古土壤沉积物所含磁性矿物的主要成分是磁铁矿,它是研究极性转换期间地球磁场详细结构的良好物质.     

首都减灾圈的思考与建议

首都减灾圈,系首都减轻自然灾害预测防治圈。1991年12月20—21日在北京召开了《首都圈自然灾害及其减灾对策研讨会》。本文根据这次会议所提供的材料,在从整体上实现减灾的思想指导下,就首都减灾圈的“成灾背景”、“首都减灾圈的组成”、“首都减灾圈的灾害预测与防治状况”、“首都减灾圈的灾害关联性分析”、“首都减灾圈的减灾实效预估”和“首都减灾圈的减灾对策与实施”等六个带有共同性的问题,进行了讨论。供制定首都圈减灾方案时参考。     

新疆地震预报研究

新疆的地震预报是1970年开始的。二十多年来在“边观测、边研究、边预报”、“多路探索、多兵种联合作战”、“走综合预报之路”的方针指导下取得了长足的进展。建设了43个地震台站,投入189套仪器,建成了遍布全疆的地震无线通讯网,开展了历史地震调查和地震烈度区划工作,建立健全了地震会商预报制度,进行了多方面的地震预报研究。实际地震预报统计分析表明,扣除自然发震概率之后,趋势预报的成功率约0.3,短临预报的成功率约0.1。在前兆台网控制范围内取得了一些震例,说明地震确实是有前兆的,但又是非常复杂的。地震预报作为科学难题还有漫长的路要走。作者简要地讨论了地震观测的间接性和地震异常的离散性、难以区别的地震异常和地壳变动异常、建立在复杂现象基础上的地震前兆以及地震的混沌性对地震预报的影响等科学问题。     

新疆第四纪盐类矿产的形成和控制因素——兼论“高山深盆”成盐模式

根据袁见齐教授“高山深盆”成盐模式,探讨新疆天山对第四纪盐类矿床的形成和控制作用,阐述了天山地貌、气候、水文特征与盐类矿产的分布规律和特征。认为“高山深盆”并非一定是四周环山的深盆,可以是某一高山与其间深盆或两侧盆地的有机组合。高山的屏障作用造成了垂直的气候分带,在潮湿多雨的山区利于成盐组份的析出并迁移到干旱少雨的闭流深盆中,形成盐类矿床。     

大别杂岩中混合岩的矿物空间分布研究

介绍了矿物空间分布研究的基本原理及两种统计方法(接触频数法和线切法)。作者对混合岩矿物空间分布的研究表明:(1)前人提出的统计方法存在方法上的缺陷和应用上的局限性,作者推导出矿物接触类型的概率公式;(2)部分浅色体中矿物显示聚集分布的特征,而绝大部分的浅色体中矿物具有分散分布的特点。结合质量平衡和地球化学研究认为:大别杂岩中主体混合岩成因机制是重熔和交代作用。     

稀土矿床：基本特点与全球分布规律

本文介绍了全球稀土资源供需历史、现状和对未来的展望。从矿床成因视角切入，将稀土矿床分为内生和外生两大类型，其中内生稀土矿床包括碳酸岩型、碱性岩型、碱性岩型- 碳酸岩型、氧化铁铜金型、热液脉型，外生稀土矿床包括风化壳离子吸附型、沉积岩型、沉积矿产(煤矿、铝土矿和沉积磷矿)伴生型、砂矿和现代海洋底部含稀土的锰结核、结壳和软泥型。归纳总结了主要类型矿床的基本特点和时空分布；认为内生稀土矿床产出于四类构造环境，包括裂谷环境、碰撞后伸展环境、大陆碰撞环境和后俯冲伸展环境；从构造演化入手，探讨了在外生与内生地质过程中稀土元素的迁移和富集规律，建立了涵盖主要矿床类型的构造- 成矿模型。