Automated Segmentation of Resistivity Image Logs Using Wavelet Transform

We describe a wavelet-transform-based method for automated segmentation of resistivity image logs that takes into account the apparent dip in the data and addresses the problem of discriminating lithofacies boundaries from noise and intrafacies variations. Our method can be applied to borehole measurements in general, but might have an advantage when applied to resistivity image logs as it addresses explicitly the large variability in facies segments recorded with a high-resolution multiple-sensor tool. We have developed an algorithm based on this method that might outperform other existing segmentation methods in the cases of low to moderate dip. We made a detailed comparison of the segmentation from our method with the one done by a geologist to delineate different lithofacies blocks in a well drilled in a deepwater depositional environment. Our results show considerable success rates in reproducing the geologically defined lithofacies boundaries, and the generality of our procedure suggests it could also be applied to other depositional environments.     

向下延拓在深部矿产勘探中的应用——以青海某矿区为例

勘测深部盲矿体时,矿体距离观测点较远,磁异常的相互叠加常常在地面形成具有一个磁异常中心的宽缓磁异常,分辨不清矿致异常和背景异常。通常用向上延拓、方向导数等处理的效果差强人意。笔者采用向下延拓方法处理了青海某矿区的宽缓磁异常,向下延拓深度为30倍点距,获得了精细的各矿体异常特征。通过与大量的见矿、未见矿钻孔比较,向下延拓场反映的磁异常位置与已见矿钻孔有较好的对应关系,同时指示出全部的未见矿钻孔都处于向下延拓场反映的磁异常外围。基于当前向下延拓方法的稳定计算问题已经取得突破,笔者认为针对一些大埋藏深度的宽缓磁异常采用向下延拓处理,区分叠加异常,圈定异常位置,可以获得事半功倍的效果,为钻孔布置提供佐证。同时对勘探精度也作了一些有意义的讨论。     

An optimal extraction algorithm for imaging photometry

This paper is primarily an investigation of whether the 'optimal extraction' techniques used in CCD spectroscopy can be applied to imaging photometry. It is found that using such techniques provides a gain of around 10 per cent in signal-to-noise ratio over normal aperture photometry. Formally, it is shown to be equivalent to profile fitting, but offers advantages of robust error estimation, freedom from bias introduced by mis-estimating the point spread function, and convenience. In addition some other techniques are presented, which can be applied to profile fitting, aperture photometry and the 'optimal' photometry. Code implementing these algorithms is available at http://www.astro.keele.ac.uk/~timn/.     

基于LabVIEW虚拟仪器技术的海洋环境监测数据采集与处理系统

介绍一种基于LabVIEW虚拟仪器技术的数据采集与处理系统。该系统以单片机和计算机为硬件平台,采用高精度16位Σ-ΔA/D转换器AD7705对采集到的信号进行A/D转换,通过串口进行数据通信,实现了叶绿素浓度的计算机自动采集、实时数据显示、分析和数据管理等功能。该实现方案还可用于温度、盐度等参数的高精度自动数据采集系统,在海洋环境监测中具有推广应用的现实前景。     

2011年3月日本福岛核电站核泄漏在海洋中的传输

使用全球版本的迈阿密等密度海洋环流模式对2011年3月日本福岛核电站泄漏在海洋中的传输以及扩散进行了数值模拟。数值模式中核废料(示踪物)排放情景采取等通量连续排放,排放时间从3月25日开始,分别持续20 d以及1 a,两种情形分别积分20 a。为了减少海洋环流年际变化带来的数值模拟的的不确定性,20 a的模式积分分别用2010年、1991-2011年、1971-1991年以及1951-1971年4个不同时段的NCEP/NCAR逐日再分析资料作为大气强迫场,因此每种排放情形包含4个数值试验。模拟结果的分析表明,不同核废料排放情景及其在不同时段大气资料对海洋模式的驱动下,模拟的示踪物总体的传输扩散路径(包括表层以及次表层)、传输速率以及垂直扩展的范围没有显著的差异。集合平均数值模拟的结果显示:在两种排放情景下,日本福岛核泄漏在海洋的传输路径受北太平洋副热带涡旋洋流系统主导,其传输路径首先主要向东,到达东太平洋后,再向南向西扩散至西太平洋,可能在10~15 a左右影响到我国东部沿海海域,且海洋次表层的传输信号比表层信号早5 a左右。通过进一步分析模式积分过程中最大示踪物浓度随时间变化发现,在积分第20 a(2031年3月),海洋表层和次表层浓度的最高值分别只有模式积分第一年浓度的0.1%和1%。在积分的20 a里,排放的核废料主要滞留在北太平洋海域(超过86%±1.5%的核废料在积分结束时,滞留在北太平洋),而在积分的前10 a(2021年之前),几乎所有的核废料滞留在北太平洋;在核废料的垂直分布上,主要集中在海洋表层至600 m的深度,在积分的20 a时间里,没有核废料信号扩散至1 000 m以下的深度。数值模拟的结果也表明核废料浓度减弱的强度以及演变的时间特征主要受洋流系统的影响,与排放源的排放时间长短关系不大。值得指出的是,更加准确地评估一个真实的核泄漏事故对海洋环境所造成的可能影响,还需要考虑大气中的放射性物质的沉降以及海洋生态对核物质的响应。     

Estimation from incomplete data files based on magnitude-frequency relation for earthquakes

The problem of incorporating the available seismological information provided by the major events of the historical catalog with those for the short period of instrumental data is investigated. Assuming that the frequency-magnitude law of Gutenberg and Richter is valid for both periods, an estimation procedure for the main parameter of this law and the rate of earthquake occurrence for historical period is presented. Application of the proposed method is demonstrated, using both real and simulated data. An extension to allow for variable quality of complete data with different magnitude values is also included.     

An enlarged and updated internally consistent thermodynamic dataset with uncertainties and correlations: the system K2O–Na2O–CaO–MgO–MnO–FeO–Fe2O3–Al2O3–TiO2–SiO2–C–H2–O2

We present, as a progress report, a revised and much enlarged version of the thermodynamic dataset given earlier (Holland & Powell, 1985). This new set includes data for 123 mineral and fluid end-members made consistent with over 200 P–T–X_CO2–f_O2 phase equilibrium experiments. Several improvements and advances have been made, in addition to the increased coverage of mineral phases: the data are now presented in three groups ranked according to reliability; a large number of iron-bearing phases has been included through experimental and, in some cases, natural Fe:Mg partitioning data; H₂O and CO₂ contents of cordierites are accounted for with the solution model of Kurepin (1985); simple Landau theory is used to model lambda anomalies in heat capacity and the Al/Si order–disorder behaviour in some silicates, and Tschermak-substituted end-members have been derived for iron and magnesium end-members of chlorite, talc, muscovite, biotite, pyroxene and amphibole. For the subset of data which overlap those of Berman (1988), it is encouraging to find both (1) very substantial agreement between the two sets of thermodynamic data and (2) that the two sets reproduce the phase equilibrium experimental brackets to a very similar degree of accuracy. The main differences in the two datasets involve size (123 as compared to 67 end-members), the methods used in data reduction (least squares as compared to linear programming), and the provision for estimation of uncertainties with this dataset. For calculations on mineral assemblages in rocks, we aim to maximize the information available from the dataset, by combining the equilibria from all the reactions which can be written between the end-members in the minerals. For phase diagram calculations, we calculate the compositions of complex solid solutions (together with P and T) involved in invariant, univariant and divariant assemblages. Moreover we strongly believe in attempting to assess the probable uncertainties in calculated equilibria and hence provide a framework for performing simple error propagation in all calculations in thermocalc, the computer program we offer for an effective use of the dataset and the calculation methods we advocate.     

The journal Communications in Numerical Methods in Engineering with Biomedical Applications becomes the International Journal for Numerical Methods in Biomedical Engineering (IJNMBE) from 1st January 2010

Communications in Numerical Methods in Engineering, founded by Roland W. Lewis in 1985, will change its title to the ‘International Journal for Numerical Methods in Biomedical Engineering’ and has a revised Aims and Scope. Copyright © 2009 John Wiley & Sons, Ltd.     

Experimental and theoretical transition probabilities in singly ionized gold

Absolute transition probabilities have been measured for lines originating from the 5d⁹6d and 5d⁹7s electronic configurations in the spectrum of singly ionized gold (Au  ii ). The laser-induced breakdown spectroscopy has been applied to free gold atoms and ions produced by laser ablation. Absolute transition probabilities have been determined using the branching fraction and the Boltzmann plot methods. Theoretical branching fractions as well as radiative lifetime values have also been obtained by a relativistic Hartree–Fock method taking core polarization and configuration interaction effects into account. The new results are compared with previous results when available.