Textural segmentation of digital rock images into bedding units using texture energy and cluster labels

The texture of digital rock images, as recorded, for instance, with borehole imaging devices, is shown to reflect different bedding types. Textural segmentation of borehole images, therefore, subdivides the recorded sequence into bedding units. We show that a textural segmentation algorithm based on the concept of texture energy achieves good results when compared with synthetic as well as real data in which petroleum geologists have performed zonations on cores. Texture energy involves filtering of the original image with a set of texture sensitive masks. The filtering is done as a finite convolution over the size of the masks. On the resulting images the variance is computed over a relatively large sliding window, which, in its practical implementation, covers the full width of the image. The resulting nine one-dimensional curves are then clustered hierarchically into a user-determined number of image texture or lithological bedding classes. Principal component analysis previous to clustering can be used to reduce redundancy in the data. A recurring and relatively ill-defined problem in this field are macro-textures, i.e., the cyclic interbedding of two or more bedding types. We show that sliding Fourier transforms and variable mask scale can successfully address the zonation of macro-textures. In general, the method gives best results with mask sizes equivalent to 2–4 centimeters, reflecting the length scale at which the investigated geological bedding seems to have its highest variation.     

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.     

“Tripoli”的含义及译名的讨论

“Tripoli”是一种新型的石英质工业矿物原料,在我国长期被误译为硅藻土,影响地质部门对这一新矿种的认识和寻找。本文追溯了“Tripoli”一名的由来,正确地阐述了“Tripoli”的本质特征,分析了在国内文献中长期被误译的根源,指明了其正确含义,并建议废除各种不正确的中译名,译为“粉石英”。     

Experimental investigation of the elastic modulus of a fractal system—A model of fractured rocks

The elastic properties of a physical model representing a damaged rock matrix were studied using a square lattice deformed under tensile stress. The elastic modulusM of such a system varies in agreement with percolation theory as|x–x _c | ^f , wherex is the damage parameter andx _c the threshold value of the damage parameter,f3.6. Atxx _c the scale dependence ofM can be expressed asML ^–f/v , whereL is the size of the sample andv the correlation exponent in percolation theory.The experimental results are of interest in assessing elastic properties in earthquake focal zones and fault zones in general.     

Na depletion in modern adakites via melt / rock reaction within the sub-arc mantle

Interaction between slab-derived melt and mantle peridotite and the role of slab melt as a metasomatizing agent in the sub-arc mantle is being increasingly recognized. Adakite, the slab melt erupted on the surface, usually exhibits anomalously high MgO, CaO, Cr and Ni contents that indicate interaction with mantle peridotitite. Here we note that Cenozoic adakites have Na₂O contents below 5.8 wt.% with ∼95% samples lower than 5.0 wt.%, and are generally depleted in this component relative to experimental basalt partial melts (mostly beyond 5.0 wt.% and up to 9.0 wt.% Na₂O) produced under 1.5-3.0 GPa conditions that are most relevant to adakite production. We interpret the adakite Na depletion to be also a consequence of the melt / rock reaction that takes place within the hot mantle wedge. During ascent and reaction with mantle peridotite, primary adakite melts gain mantle components MgO, CaO, Cr and Ni but lose Na₂O, SiO₂ and perhaps K₂O to the mantle, leading to Na-rich mantle metasomatism. Selective assimilation of predominately mantle clinopyroxene, some spinel and minor olivine at high T/P has been considered to be an important process in producing high-Mg adakites from primary low-Mg slab melts [Killian, R., Stern, C. R., 2002. Constraints on the interaction between slab melts and the mantle wedge from adakitic glass in peridotite xenoliths. Eur. J. Mineral. 14, 25-36]. In such a process, Na depletion in the assimilated melt is the result of dilution due to the increase in melt mass. Phase relationships in the reaction system siliceous melt + peridotite and quantitative calculation suggest that assimilation of mantle clinopyroxene, olivine and spinel and fractional crystallization of sodic amphibole and orthopyroxene, under conditions of moderate T/P and increasing melt mass, is also an important process that modifies the composition of adakites and causes the Na depletion.     

Cross-formational hydrocarbon fluid flows in the Tertiary deltaic system of the Beaufort–Mackenzie Basin

Molecular and isotopic compositions of crude oils in the Beaufort–Mackenzie Basin confirm three genetic end-member oil groups and suggest extensive cross-formational hydrocarbon fluid flows in the Tertiary deltaic system. Inter- and intra-fractional variations in the geochemistry of the Tertiary-reservoired oils indicate that the oil source/maturity signatures were substantially masked by biomarkers that were picked up along migration pathways. Thus, many of the previously recognized “immature non-marine oils” are in fact thermally mature, probably derived from unpenetrated deeper marine source rocks. Although the effective source rock volumes have not been evaluated and their exact stratigraphic levels remain unknown, the relative timing of oil generation versus trap formation, rather than poor source quality, may be the cause of under-filled traps in the offshore area.     

Biostratigraphy of Jurassic-Cretaceous boundary sediments in the Eastern Crimea

The Dvuyakornaya Formation section in the eastern Crimea is described and subdivided into biostratigraphic units based on ammonites, foraminifers, and ostracodes. The lower part of the formation contains first discovered ammonites of the upper Kimmeridgian (Lingulaticears cf. procurvum (Ziegler), Pseudowaagenia gemmellariana Olóriz, Euvirgalithacoceras cf. tantalus (Herbich), Subplanites sp.) and Tithonian (?(Lingulaticeras efimovi (Rogov), Phylloceras consaguineum Gemmellaro, Oloriziceras cf. schneidi Tavera, and Paraulacosphinctes cf. transitorius (Oppel)). Based on the assemblage of characteristic ammonite species, the upper part of the formation is attributed to the Berriasian Jacobi Zone. Five biostratigraphic units (zones and beds with fauna) distinguished based on foraminifers are the Epistomina ventriosa-Melathrokerion eospirialis Beds and Anchispirocyclina lusitanica-Melathrokerion spirialis Zone in the upper Kimmeridgian-Tithonian, the Protopeneroplis ultragranulatus-Siphoninella antiqua, Frondicularia cuspidiata-Saracenaria inflanta zones, and Textularia crimica Beds in the Berriasian. The Cyrherelloidea tortuosa-Palaeocytheridea grossi Beds of the Upper Jurassic and Raymoorea peculiaris-Eucytherura ardescae-Protocythere revili Beds of the Berriasian are defined based on ostracodes. A new biostratigraphic scale is proposed for the upper Kimmeridgian-Berriasian of the eastern Crimea. The Dvyyakornaya Formation sediments are considered as deepwater facies accumulated on the continental slope.     

Correlations in aftershock and seismicity patterns

Correlations in space and time play a fundamental role in earthquake processes. One direct manifestation of the effects of correlations is the occurrence of aftershocks due to the stress transfer in the vicinity of a main shock. Less obvious and more speculative changes in correlations may occur in the background seismicity before large earthquakes. Using statistical physics it is possible to introduce a measure of spatial correlations through a correlation length. This quantity characterizes how local fluctuations can influence the occurrence of earthquakes over distances comparable with the correlation length. In this work, the physical basis of spatial correlations of earthquakes is discussed in the context of critical phenomena and the percolation problem. The method of two-point correlation function is applied to the seismicity of California. Well defined variations in time of the correlation length are found for aftershock sequences and background seismicity. The scaling properties of our obtained distributions are analyzed with respect to changes in several scaling parameters such as lower magnitude cutoff of earthquakes, the maximum time interval between earthquakes, and the spatial size of the area considered. This scaling behavior can be described in a unified manner by utilizing the multifractal fit. Utilizing the percolation approach the time evolution of clusters of earthquakes is studied with the correlation length defined in terms of the radius of gyration of clusters. This method is applied to the seismicity of California.     

Évolution de l'onde semi-diurne M2 de la marée à Brest de 1846 à 2005

The work of searching, recovering and quality control of ancient sea-level measurements at Brest is presented. This work enables us to complete a study carried out by Cartwright in 1972, which showed a decrease in the tidal M2 semi-diurnal amplitude of 1% per century. After including these ancient data, as well as the last four decades of observations in the analysis, our results show an increase of the amplitude of M2 after 1960 and a decrease before 1880, suggesting a long-period oscillation rather than a steady secular trend. To cite this article: N. Pouvreau et al., C. R. Geoscience 338 (2006).     

冲击凿岩钎杆的瞬态动力学分析

以ANSYS为平台，按照真实工作状态对凿岩钎杆进行了瞬态动力学分析，并对不同活塞形状下分析结果数据进行了比较，为冲击凿岩系统设计提供了理论依据。