基于QAPM矿物模型遗传算法评价火成岩储层

在松辽盆地深层发现了含气火成岩储层。由于火成岩矿物组成复杂和含量的变化,使得选择用于测井评价的解释参数很困难。基于IUGS提出的QAPF分类方案,本文提出了采用遗传算法,利用测井数据确定火成岩矿物含量的方法。根据QAPF分类方案,将火成岩中的矿物分为五类:Q-石英;A-碱性长石;P-斜长石和方柱石;F-副长石(研究区未出现);M-铁镁矿物。本文提出用包括孔隙度在内的QAPM模型对储层进行分析。建立密度、视中子孔隙度、声波时差、自然伽玛和体积光电吸收截面指数的测井响应方程,各矿物参数从斯伦贝谢的矿物参数手册中得到。用遗传算法计算骨架中四种矿物的体积,根据四种矿物的体积含量,依据QAPF分类对火成岩命名。基于解释参数计算的孔隙度可与岩心分析的孔隙度相比,本文给出的火成岩命名与岩心化学分析的命名相一致。     

Evaluation of volcanic reservoirs with the “QAPM mineral model” using a genetic algorithm

Gas-bearing volcanic reservoirs have been found in the deep Songliao Basin, China. Choosing proper interpretation parameters for log evaluation is difficult due to complicated mineral compositions and variable mineral contents. Based on the QAPF classification scheme given by IUGS, we propose a method to determine the mineral contents of volcanic rocks using log data and a genetic algorithm. According to the QAPF scheme, minerals in volcanic rocks are divided into five groups： Q（quartz）, A （Alkaline feldspar）, P （plagioclase）, M （mafic） and F （feldspathoid）. We propose a model called QAPM including porosity for the volumetric analysis of reservoirs. The log response equations for density, apparent neutron porosity, transit time, gamma ray and volume photoelectrical cross section index were first established with the mineral parameters obtained from the Schlumberger handbook of log mineral parameters. Then the volumes of the four minerals in the matrix were calculated using the genetic algorithm （GA）. The calculated porosity, based on the interpretation parameters, can be compared with core porosity, and the rock names given in the paper based on QAPF classification according to the four mineral contents are compatible with those from the chemical analysis of the core samples.     

莫桑比克东部滨海砂矿矿物组成及特征

滨海砂矿是目前全球最大的潜在矿产资源之一。以莫桑比克的滨海砂矿作为研究对象,利用光学显微镜、X射线衍射、扫描电镜、等离子体质谱（ICP-MS）等手段对莫桑比克东部滨海砂矿的矿物组成及特征进行了研究。并用数理统计方法讨论了重砂矿的元素赋存状态及其载体矿物。研究结果表明：（1）砂矿的矿石矿物主要为钛铁矿、金红石和锆石,脉石矿物主要为石英和长石。（2）重砂矿的主要化学组分为Ti和zr;Hf及REE为次要组分;而主要有害元素组分Th和U,其总体含量低,故重砂矿对采选加工和产品利用均无明显环境危害。（3）重砂矿中元素含量与粒度之间呈一定相关性,显示zr、Hf、U、Cr、Ta的载体矿物的粒度相对较细;Ti、Fe、V、Nb的载体矿物的粒度相对较粗。     

Mineral classification from quantitative X-ray maps using neural network: Application to volcanic rocks

We developed a mineral classification technique of electron probe microanalyzer (EPMA) maps in order to reveal the mineral textures and compositions of volcanic rocks. In the case of lithologies such as basalt that include several kinds of minerals, X-ray intensities of several elements derived from EPMA must be considered simultaneously to determine the mineral map. In this research, we used a Kohonen self-organizing map (SOM) to classify minerals in the thin-sections from several X-ray intensity maps. The SOM is a type of artificial neural network that is trained using unsupervised training to produce a two-dimensional representation of multi-dimensional input data. The classified mineral maps of in situ oceanic basalts of the Juan de Fuca Plate allowed us to quantify mineralogical and textural differences among the marginal and central parts of the pillow basalts and the massive flow basalt. One advantage of mineral classification using a SOM is that relatively many minerals can be estimated from limited input elements. By applying our method to altered basalt which contains multiple minerals, we successfully classify eight minerals in thin-section.     

基于逾渗模型的泥页岩导电机理模拟

泥页岩储层储集空间以纳米/微米级的微孔隙和微裂缝为主,结构复杂,加之干酪根和黏土矿物含量较高,并含有一定量的黄铁矿等导电矿物,造成储层导电机理异于常规储层,岩电实验I-S_w曲线呈现非线性特征,阿尔奇公式等传统评价模型适用性较差.针对上述问题,根据实际岩心实验资料,结合随机算法建立三维逾渗模型并通过数值模拟和超松弛迭代法进行求解,分析泥页岩储层非阿尔奇性产生原因以及泥页岩储层电性的影响因素及规律.模拟结果显示,岩石孔隙拓扑结构和形状尺寸、矿物组成以及地层水电阻率等因素均对泥页岩储层电阻率产生不同程度的影响,通过改变上述因素的设定值,可以建立储层电阻率与各因素的单相关关系,并据此建立修正模型,实现储层含水饱和度的计算.该模型在四川某页岩气产区取得了较好的效果,具有良好的应用前景,为利用逾渗模型模拟方法解决油田勘探开发中的复杂问题提供了新思路.     

Chemical divides and evaporite assemblages on Mars

Assemblages of evaporite minerals record detailed physical and chemical characteristics of ancient surficial environments. Accordingly, newly discovered regions of saline minerals on Mars are high priority targets for exploration. The chemical divide concept of evaporite mineral formation is used successfully to predict evaporite mineralogy and brine evolution on Earth. However, basaltic weathering largely controls fluid compositions on Mars and the robust predictive capabilities of terrestrial chemical divides cannot be used to interpret Martian evaporites. Here we present a new chemical divide system that predicts evaporite assemblages identified in SNC-type meteorites, ancient evaporites discovered on Meridiani Planum by the Opportunity rover, and Mars Express OMEGA data. We suggest that a common fluid type that has been buffered to different pH levels by basaltic weathering controls the variability among Martian evaporite assemblages and that evaporite mineralogy and brine evolution is essentially established by the initial composition of the dilute evaporating fluid.     

基于SVM测井数据的火山岩岩性识别——以辽河盆地东部坳陷为例

辽河盆地东部坳陷储集层由火山多期喷发形成,岩相岩性复杂,岩性以中、基性火山岩为主.本文将火山岩的岩心及岩矿鉴定资料与测井数据进行整合,应用测井数据建立支持向量机(SVM)两分类和多分类岩性识别模式.首先,深入研究支持向量机二分类及"一对一"、"一对多"和有向无环图三种经典多分类算法的基本原理及结构;然后,总结研究区域火山岩岩石特征,分析测井数据的测井响应组合特征,选择40口井中岩心分析和薄片鉴定资料完整、常规五种测井曲线(RLLD,CNL,DEN,AC,GR)齐全的1200个测井数据作为训练样本,构造三种支持向量机岩性识别模式;最后,对4测试井中800个测井数据进行岩性识别,识别结果与取心段岩心描述和岩心/岩屑薄片鉴定资料对比,实验结果表明有向无环图更适合辽河盆地火山岩的识别,识别正确率达到82.3%.     

On Earth’s Mantle Constitution and Structure from Joint Analysis of Geophysical and Laboratory-Based Data: An Example

Determining Earth’s structure is a fundamental goal of Earth science, and geophysical methods play a prominent role in investigating Earth’s interior. Geochemical, cosmochemical, and petrological analyses of terrestrial samples and meteoritic material provide equally important insights. Complementary information comes from high-pressure mineral physics and chemistry, i.e., use of sophisticated experimental techniques and numerical methods that are capable of attaining or simulating physical properties at very high pressures and temperatures, thereby allowing recovered samples from Earth’s crust and mantle to be analyzed in the laboratory or simulated computationally at the conditions that prevail in Earth’s mantle and core. This is particularly important given that the vast bulk of Earth’s interior is geochemically unsampled. This paper describes a quantitative approach that combines data and results from mineral physics, petrological analyses of mantle minerals, and geophysical inverse calculations, in order to map geophysical data directly for mantle composition (major element chemistry and water content) and thermal state. We illustrate the methodology by inverting a set of long-period electromagnetic response functions beneath six geomagnetic stations that cover a range of geological settings for major element chemistry, water content, and thermal state of the mantle. The results indicate that interior structure and constitution of the mantle can be well-retrieved given a specific set of measurements describing (1) the conductivity of mantle minerals, (2) the partitioning behavior of water between major upper mantle and transition-zone minerals, and (3) the ability of nominally anhydrous minerals to store water in their crystal structures. Specifically, upper mantle water contents determined here bracket the ranges obtained from analyses of natural samples, whereas transition-zone water concentration is an order-of-magnitude greater than that of the upper mantle and appears to vary laterally underneath the investigated locations.     

Correlation and petrogenesis of silicic pyroclastic rocks in the Northern Pannonian Basin,Eastern-Central Europe: In situ trace element data of glass shards and mineral chemical constraints

The Neogene volcanism of the Carpathian–Pannonian Region, Eastern-Central Europe started with repeated explosive eruption of silicic magmas. Volcanic products consist mostly of non-welded and partially to densely welded ignimbrites, which cover large areas in the Pannonian Basin. Since this volcanism occurred during a long time interval, from 21 Ma to 13.5 Ma ago, these pyroclastic deposits have great stratigraphic importance, as well as providing valuable information about petrogenetic processes during the formation of the back-arc basin area. In this study, we used in situ trace element data of glasses (glass shards, pumice and fiamme glass) determined by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS), combined with microprobe geochemical data of the main mineral phases (plagioclase, biotite), to correlate the scattered outcrops of the ignimbrite units. Based on these data, we distinguished four ignimbrite units in the Bükkalja Volcanic Field, Northern Pannonian Basin. Each of these units is characterized by specific geochemical fingerprints. Thorium, Nb, Y and the rare earth elements are the most effective discriminator elements in glasses. The modal composition of mineral phases (occurrence or lack of certain minerals) and chemistry of plagioclases and biotites are also good correlation tools, especially the Fe, Mg and Ti contents of biotites. We suggest that these correlation criteria, particularly the grain-specific LA-ICP-MS data, could be effectively used also in other ignimbrite fields and in wider tephrochronological studies. The in situ trace element composition of glasses, representing the liquid part of the erupted magma, can be also used to constrain the petrogenesis of the rhyolitic magmas. Trace element ratios such as La/Nb, La/Y and Th/Nb suggest the importance of minor (e.g., hornblende and ilmenite) and accessory (e.g., zircon, allanite) minerals controlling the composition of the erupted melt. Rhyolitic magmas probably evolved from metaluminous andesitic parental melts via fractional crystallization. Syn-eruptive magma mingling was detected in the genesis of the Middle Ignimbrite Unit, based on the strong intra-sample geochemical variation both in the glasses and in the phenocrysts.     

核磁共振测井孔隙度观测模式与处理方法研究

核磁共振（NMR）测井是确定地层孔隙度的有效方法,已经得到广泛应用.但是,在我国部分油田的应用中发现核磁测井孔隙度与常规孔隙度之间存在明显差异.本文对NMR测井孔隙度观测模式和处理方法进行研究,分析了NMR孔隙度误差的来源.针对富含蒙脱石、伊利石和伊/蒙混层粘土矿物的陆相沉积砂泥岩地层,根据这些粘土矿物的NMR横向弛豫时间T2值特征以及孔隙度处理结果得到的40 ms附近孔隙度分配量不合适,提出了新的观测模式和处理方法,并用实验验证了新方法的可行性.     

Calculation of oxygen isotope fractionation in hydroxyl-bearing silicates

The modified increment method has been applied to the calculation of oxygen isotope fractionation factors for hydroxyl-bearing silicate minerals. The order of ¹⁸O enrichment obtained in common rock-forming minerals is: pyrophyllite > kaolinite > tourmaline talc > prehnite topaz > illite > phengite > lepidolite muscovite staurolite > epidote > glaucophane > serpentine chlorite > tremolite > hornblende > phlogopite biotite > humite > norbergite > ilvaite. Hydroxyl-bearing silicates are enriched in ¹⁸O relative to hydroxyl groups but depleted in ¹⁸O relative to anhydrous counterparts.

Three sets of self-consistent fractionation factors: between quartz and the hydroxyl-bearing silicate minerals, between calcite and the silicate minerals, and between the silicate minerals and water, have been calculated for a temperature range of 0–1200°C. The fractionation factors calculated for mineral pairs are applicable to isotopic geothermometry in igneous, metamorphic and sedimentary petrology. They can be used as a test of isotopic equilibrium or disequilibrium in natural mineral assemblages over all temperature ranges of geological interest. The difference in oxygen isotope composition between the hydroxyl-bearing mineral and the OH group is quantitatively demonstrated to be temperature dependent and, therefore, can be used as a single-mineral geothermometer.     

High precision in-situ Pb isotopic analysis of sulfide minerals by femtosecond laser ablation multi-collector inductively coupled plasma mass spectrometry

An in-situ microanalysis of Pb isotopic compositions in sulfide minerals is carried out by using femtosecond laser-ablation multi-collector inductively coupled plasma mass spectrometry(fs LA-MC-ICP-MS).High-temperature-activated carbon was used to filter Hg contained in the carrier gas,which reduced the Hg background signal by 48%and also lowered the detection limit of the analysis.Fractionation and mass discrimination effects existing in the ICP-MS analytical processes were corrected using an internal reference Tl in conjunction with an external reference NIST SRM 610.The proposed method was used to analyze the Pb isotopic compositions of chalcopyrite,pyrite,and sphalerite from the Dulong Sn-Zn-In polymetallic ore district.The results showed that in this ore district,the sulfide minerals and different grains of the same sulfide mineral show a large variation in Pb content up to 1000-fold.The studied pyrites show relatively higher Pb contents and homogeneous Pb isotopic compositions,whereas the sphalerites have low Pb contents but most variable Pb isotopic compositions.It is suggested that the large variation of Pb isotopic composition may reflect a late hydrothermal superimposition on the primary sulfide formation.In addition,radiogenic Pb accumulated by radioactive decay of trace amounts of U over time in the host minerals may also be one of the causes for the large variation range of Pb content and Pb isotopic composition of those low-Pb sphalerites.Chalcopyrite and sphalerite grains with Pb content greater than 10 ppm presented a consistent Pb isotopic distribution,whereas all the sulfide grains with Pb content greater than 100 ppm had consistent Pb isotopic composition within 2s measurement uncertainties.The in-situ analysis of Pb isotopic composition agreed well with the results obtained by conventional chemical methods within2s measurement uncertainties,indicating that the data obtained by fs LA-MC-ICP-MS are reliable.Additionally,this study indicates that the Pb isotopic composition could truthfully record the source of ore-forming minerals only for sulfide minerals with high Pb content.On the contrary,the Pb isotopic composition of low-Pb sulfide minerals may be affected by trace amounts of U in the host minerals that may lead to a highly radiogenic Pb isotope ratio.Alternatively,it is also possible that late fluid metasomatic overprinting may alter the Pb isotopic compositions.     

DETECTION OF THIN BEDS WITH THE PSEUDO-ACOUSTIC IMPEDANCE SECTION*

The conventional seismic response of a thin bed approximates the time derivative of the incident wavelet, whereas the pseudo-impedance response approximates the incident wavelet. Consequently the pseudo-impedance response of a geological sequence composed of thin beds is simpler and easier to interpret than the conventional response. By calibrating the sonic log data with check-shot data and performing zero-phase seismic processing, the fit of the sonic log and pseudo-velocity section is improved. Discrepancies in amplitude and phase, however, generally remain. A five-step processing and interpretation procedure, which benefits from multichannel interpretation along the model seismic section generated from the sonic logs, is described. The method has been tested with field data. In the test the detection of thin beds and the estimation of the natural gas content was more reliable with the proposed procedure than with the conventional method.     

On the Causes of Electrical Conductivity Anomalies in Tectonically Stable Lithosphere

Magnetotelluric (MT) data can image the electrical resistivity of the entire lithospheric column and are therefore one of the most important data sources for understanding the structure, composition and evolution of the lithosphere. However, interpretations of MT data from stable lithosphere are often ambiguous. Recent results from mineral physics studies show that, from the mid-crust to the base of the lithosphere, temperature and the hydrogen content of nominally anhydrous minerals are the two most important controls on electrical conductivity. Graphite films on mineral grain boundaries also enhance conductivity but are stable only to the uppermost mantle. The thermal profile of most stable lithosphere can be well constrained, so the two important unknowns that can affect the conductivity of a lithospheric section are hydrogen content and graphite films. The presence of both of these factors is controlled by the geological history of the lithosphere. Hydrogen in nominally anhydrous minerals behaves as an incompatible element and is preferentially removed during melting or high-temperature tectonothermal events. Grain-boundary graphite films are only stable to ~900 °C so they are also destroyed by high-temperature events. Conversely, tectonic events that enrich the lithosphere in incompatible elements, such as interaction with fluids from a subducting slab or a plume, can introduce both hydrogen and carbon into the lithosphere and therefore increase its electrical conductivity. Case studies of MT results from central Australia and the Slave Craton in Canada suggest that electrical conductivity can act as a proxy for the level of enrichment in incompatible elements of the lithosphere.     

First-principles calculations of elasticity of minerals at high temperature and pressure

The elasticity of minerals at high temperature and pressure(PT) is critical for constraining the composition and temperature of the Earth's interior and understand better the deep water cycle and the dynamic Earth. First-principles calculations without introducing any adjustable parameters, whose results can be comparable to experimental data, play a more and more important role in investigating the elasticity of minerals at high PT mainly because of(1) the quick increasing of computational powers and(2) advances in method. For example, the new method reduces the computation loads to one-tenth of the traditional method with the comparable precise as the traditional method. This is extraordinarily helpful because first-principles calculations of the elasticity of minerals at high PT are extremely time-consuming. So far the elasticity of most of lower mantle minerals has been investigated in detail. We have good idea on the effect of temperature, pressure, and iron concentration on elasticity of main minerals of the lower mantle and the unusual softening in bulk modulus by the spin crossover of iron in ferropericlase. With these elastic data the lower mantle has been constrained to have 10–15 wt% ferropericlase, which is sufficient to generate some visible effects of spin crossover in seismic tomography. For example, the spin crossover causes that the temperature sensitivity of P wave at the depth of ~1700 km is only a fraction of that at the depth of ~2300 km. The disruptions of global P wave structure and of P wave image below hotspots such as Hawaii and Iceland at similar depth are in consistence with the spin crossover effect of iron in ferropericlase. The spin crossover, which causes anomalous thermodynamic properties of ferropericlase, has also been found to play a control role for the two features of the large low shear velocity provinces(LLSVPs): the sharp edge and high elevation up to 1000 km above core-mantle boundary. All these results clearly suggest the spin crossover of iron in the lower mantle. The theoretical investigations for the elasticity of minerals at the upper mantle and water effect on elasticity of minerals at the mantle transition zone and subducting slab have also been conducted extensively. These researches are critical for understanding better the composition of the upper mantle and water distribution and transport in the Earth's mantle. Most of these were static calculations, which did not include the vibrational(temperature) effect on elasticity, although temperature effect on elasticity is basic because of high temperature at the Earth's interior and huge temperature difference between the ambient mantle and the subducting slab. Including temperature effect on elasticity of minerals should be important future work. New method developed is helpful for these directions. The elasticity of iron and iron-alloy with various light elements has also been calculated extensively. However, more work is necessary in order to meet the demand for constraining the types and amount of light elements at the Earth's core.     

Systematics in uranium-ionium dating of sediments

The mineral separates of a core have been systematically dated by the Io/Th chronology method. The results show that authigenic minerals, though not always completely authigenic, are suitable for dating.We give some evidence against the hypothesis of Io/²³²Th isotopic equilibrium in marine conditions and suggest that the minerals are mixtures of two components.     

Rare earth abundances in chondritic phosphates and their implications for early stage chronologies

The abundances of nine rare earth elements (REE) in phosphate separates from three ordinary chondrites, Saint Séverin (LL6), Bruderheim (L6) and Richardton (H5), were measured by instrumental neutron activation analysis. All REE except europium are enriched in the phosphate minerals (merrillite and chlorapatite) by factor of 200–300 relative to the chondritic average, whereas Eu is enriched by a factor of 40–50. Electron microprobe analysis showed no significant differences in phosphate mineral composition among the three chondrites studied, though the relative proportions of two minerals varied.According to our data, REE are enriched by almost the same factor in merrillite and chlorapatite in the Bruderheim and, with less certainty, in the other two chondrites. This behavior of REE contrast with that of the actinoid elements, Th, U and Pu, which are also enriched in phosphate but are fractionated between merrillite and chlorapatite. Since Pu and REE show different fractionation behavior in chondritic phosphates, it may be difficult to use REE as stand-ins for Pu in²⁴⁴Pu chronology.     

低阻油层双侧向测井的反演研究

对于双侧向测井反演，本文将广义逆反演方法用于消除侵入影响的三参数(侵入半径、侵入带电阻率、原状地层电阻率)反演中。通过实际多口井连续曲线的反演处理表明，本文提出的反演方法令理而且有效。同时，广义逆方法提高了反演计算的稳定性和精度。对低阻油层实际资料的处理结果与试油结论一致。     

含裂缝地层的过套管电阻率测井响应特征研究

由于石油压裂开采等开发措施的实施,会使在套管周围存在着许多被油气等高阻流体填充的微裂缝.尽管这些裂缝或孔隙可能很小,但对套管井电阻率测量会有非常大的影响,甚至会改变原地层的电阻率测井特征,因此裂缝测井响应的计算及考察对过套管电阻率测井十分重要.为解决裂缝测井响应的计算问题,本文提出了计算等效电阻的电流通量管模型,利用该电流通量管模型给出了含垂直裂缝地层横向电阻的计算方法,基于传输线方程法及地层电阻的过套管测量方法实现了含裂缝地层的过套管电阻率测井响应的数值计算,通过计算实例考察了地层裂缝对过套管电阻率测井响应的影响.算例表明:裂缝中的高阻流体对地层视电阻率测量结果会有较大的影响;环形裂隙比垂直裂缝有更大的电阻率测井响应.本文的研究为解决微裂缝过套管电阻率测井响应的计算这一关键技术问题提供了一种可行的计算与考察方法.     

Alteration and mass transfer inferred from the Hirabayashi GSJ drill penetrating the Nojima Fault, Japan

Abstract Mineralogical and geochemical studies on the fault rocks from the Nojima–Hirabayashi borehole, south-west Japan, are performed to clarify the alteration and mass transfer in the Nojima Fault Zone at shallow depths. A complete sequence from the hornblende–biotite granodiorite protolith to the fault core can be observed without serious disorganization by surface weathering. The parts deeper than 426.2 m are in the fault zone where rocks have suffered fault-related deformation and alteration. Characteristic alteration minerals in the fault zone are smectite, zeolites (laumontite, stilbite), and carbonate minerals (calcite and siderite). It is inferred that laumontite veins formed at temperatures higher than approximately 100°C during the fault activity. A reverse component in the movement of the Nojima Fault influences the distribution of zeolites. Zeolite is the main sealing mineral in relatively deep parts, whereas carbonate is the main sealing mineral at shallower depths. Several shear zones are recognized in the fault zone. Intense alteration is localized in the gouge zones. Rock chemistry changes in a different manner between different shear zones in the fault zone. The main shear zone (MSZ), which corresponds to the core of the Nojima Fault, shows increased concentration of most elements except Si, Al, Na, and K. However, a lower shear zone (LSZ-2), which is characterized by intense alteration rather than cataclastic deformation, shows a decreased concentration of most elements including Ti and Zr. A simple volume change analysis based on Ti and Zr immobility, commonly used to examine the changes in fault rock chemistry, cannot account fully for the different behaviors of Ti and Zr among the two gouge zones.