Estimation of earthquake hazard parameters for incomplete and uncertain data files

The maximum likelihood estimation of earthquake hazard parameters (maximum regional magnitudem _max, activity rate , and theb parameter in the Gutenberg-Richter distribution) is extended to the cases of incomplete and uncertain data. The method accepts mixed data containing only large (extreme) events and a variable quality of complete data with different threshold magnitude values. Uncertainty of earthquake magnitude is specified by two values, the lower and upper magnitude limits. It is assumed that such an interval contains the real unknown magnitude. The proposed approach allows the combination of different quality catalog parts, e.g. those where the assignment of magnitude is questionable and those with magnitudes precisely determined.As an illustration of the method, the seismic hazard analysis for western Norway and adjacent sea area (4–8°E, 58–64°N) is presented on the basis of the strongest earthquakes felt during the period 1831–1889 and three complete catalog parts, covering the period 1890–1987.Paper presented at the 21st General Assembly of the European Seismological Commission held in Sofia, 1988.On leave from Institute of Geophysics, Polish Academy of Sciences, Warsaw, Poland.     

稀土元素内潜同晶分馏模式及其意义

稀土元素所具有的内潜同晶现象,可将其14个元素则可视为一个元素,而其它稀土元素视为其同位素,其在地质作用过程中,将像同位素一样,存在分馏。按照同位素的分馏模式,以球粒陨石值作为标准,计算出每一个稀土元素的分馏值,其值按δREE=[(REE/ΣREE)样品/(REE/ΣREE)球粒陨石-1]×100计算,以此为纵坐标,以原子序数为横坐标,作一个模式图,这种模式图称为稀土元素内潜同晶分馏模式图。该模式图在岩石成因,特别是岩浆岩的成因、岩浆岩的演化等方面有其独特的意义。     

Statistical Analysis of the Completeness of a Seismic Catalogue

Among the numerous issues that the study of seismic events presents, theincompleteness of catalogues is certainly one of the most important. It is also one that only the contribution of many and different skills canprovide with a valid solution. In this paper the search for the complete part of a catalogue is expressed in terms of identification of the changepoint in a hierarchical Bayesian model. Stochastic simulation methods, recently presented in the literature, have enabled us to overcome the computational issues that previously made this approach prohibitive. We have applied the method on data, drawn from the Italian NT4.1.1 catalogue,related to some seismogenetic zones of ZS.4 zonation within which we assumespatial incompleteness to be homogeneous. The results obtained are given inthe concluding sections of the paper.     

X射线荧光光谱微区分析在铜矿物类质同象鉴定中的应用

铜矿物种类繁多,赋存状态各异,成分复杂,由于共生矿物和伴生矿物的存在,增加了传统岩矿鉴定工作的难度。特别是铜矿物间存在的类质同象现象(黝铜矿和砷黝铜矿)和铜矿物与其他矿物间(磁黄铁矿和方黄铜矿)存在镜下光学特征相似的现象,应用偏光或反光显微镜难以给铜矿物准确定名。本文应用具有微区分析功能的X射线荧光光谱仪扫描铜矿物标本谱图,通过选择测量条件、测量方式和干扰校正模型(谱线重叠和基体效应)建立了类质同象物相的鉴定方法。该方法具有较高的精密度和准确度,对黄铜矿标本进行12次测定,S、Fe、Cu的相对标准偏差(RSD)分别为0.74%、0.86%、0.18%;对铜矿物标准样品进行单次测定,待测元素的测量结果与推荐值的相对误差均小于6%。该方法可准确甄别出黝铜矿与含银砷黝铜矿、磁黄铁矿与方黄铜矿;测量效率高,能够提供主量元素和主要伴生元素的二维或三维图像,直观地呈现矿物的分布情况,为铜矿物的准确定名和矿物的综合利用提供了有价值的信息。     

利用ICCG迭代技术加快电阻率三维正演计算

一般而言,有限差分法求解点源三维地电场正问题所形成的大型稀疏线性方程组Ａｘ＝ｂ,直接解法的计算效率极低。本文从系数矩阵Ａ的不完全Ｃｈｏｌｅｓｋｙ分解及矩阵特征值的特点等角度,说明了不完全Ｃｈｏｌｅｓｋｙ共轭梯度（ＩＣＣＧ）迭代技术可大大提高电阻率三维正演速度的内在原因。结合矩阵Ａ的稀疏存储模式,使得内存需求也大大减少。     

聚落形态演变与防洪功能变化的关系——以广东高要地区为例

高要地区特殊的地理环境与气候条件造成了历来严重的洪涝灾害。其历史洪水淹浸区和历史堤防范围等历史地理环境与该地区的“八卦”形态聚落分布存在着特殊的“同构”现象。通过图解分析与逻辑推理,解读了“同构关系”的内在成因：“八卦”形态是聚落防洪的最优形态;由于民国以前的高要地区区域防洪未完善,因此位于西江南岸缺乏堤防保护的洪涝区域内大量聚落采用了“八卦”形态的聚落被动防洪方式。高要地区“八卦”形态聚落分布与历代成堤区域的同构关系说明了区域防洪方略与聚落被动防洪方略互为补充、此消彼长的辩证关系。在区域防洪设施未完善的时代,聚落被动防洪方略作为防洪方略的重要补充,导致了聚落的特殊形态;而在区域防洪设施基本完善、有效保护聚落免受洪涝灾害影响的今天,聚落被动防洪方略的存在意义大为减弱,“八卦”形态逐渐瓦解改变。这是聚落形态演变的一个重要动因。区域防洪问题解决后,大量聚落不断地向低洼地区延伸,与当初的被动防洪方略已经互相矛盾。这种扩张是非理性、无序型的扩张。当前对聚落防灾的轻视心态,埋下了洪涝灾害的隐患。     

一种原状欠压密土的力学特性实验研究

原状土的结构性对它的力学性质影响显著。通过对原状土与重塑土多应力路径的实验结果比较,来揭示原状土体从小应变到大应变的特殊力学特性及应力路径的影响。从等压固结对比实验可看出,该原状样明显的结构性与欠压密性。对于原状欠压密土,体应变对土体结构性的破坏要大于剪应变。欠压密原状土的结构性对土体强度的影响主要表现在内摩擦角的增大。土体小应变时表现出明显的高刚度与非线性。通过不同方向切样的原状样实验结果对比,反映了原状土体存在明显的各向异性,结构性引起的各向异性在应变越小时越明显。     

Structural identification with incomplete output‐only data and independence of measured information for shear‐type systems

Structural identification is the inverse problem of estimating physical parameters of a structural system from its vibration response measurements. Incomplete instrumentation and ambient vibration testing generally result in incomplete and arbitrarily normalized measured modal information, often leading to an ill‐conditioned inverse problem and non‐unique identification results. The identifiability of any parameter set of interest depends on the amount of independent available information. In this paper, we consider the identifiability of the mass and stiffness parameters of shear‐type systems in output‐only situations with incomplete instrumentation. A mode shape expansion‐cum‐mass normalization approach is presented to obtain the complete mass normalized mode shape matrix, starting from the incomplete non‐normalized modes identified using any operational modal analysis technique. An analysis is presented to determine the minimum independent information carried by any given sensor set‐up. This is used to determine the minimum necessary number and location of sensors from the point of view of minimum necessary information for identification. The different theoretical discussions are illustrated using numerical simulations and shake table experiments. It is shown that the proposed identification algorithm is able to obtain reliably accurate physical parameter estimates under the constraints of minimal instrumentation, minimal a priori information, and unmeasured input. The sensor placement rules can be used in experiment design to determine the necessary number and location of sensors on the monitored system. John Wiley & Sons, Ltd.     

钴的赋存状态及其在找矿和资源评估中的意义

通过分析含钴矿物类型的多样性 ,得出地球系统中钴呈硫化物、类硫化物、铁氧化物和呈类质同像代替铁多种赋存状态及其受晶体化学因素控制特点。基于钴赋存状态的不同 ,提出当前钴的地球化学找矿和矿产资源评估中仅以钴量化丰度值作为化探异常和矿石品位判据的局限性。指出元素状态丰度提供资源信息的重要意义     

高精度曲面建模的一种快速算法

 高精度曲面建模(HASM)是一种全新的曲面建模方法,其整个过程可分为偏微分方程的离散、采样方程建立和代数方程组求解3个阶段。目前所采用的求解对称正定方程组的方法主要是共轭梯度法。为了解决HASM的计算速度问题,本文给出了2种新的预处理共轭梯度算法,分别为不完全Cholesky分解共轭梯度法和对称逐步超松弛预处理共轭梯度法。实验表明,不完全Cholesky分解共轭梯度法收敛速度最快,且这2种预处理方法均比其他方法具有更快的收敛速度。