Acoustic and petrophysical relationships in low-shale sandstone reservoir rocks

This paper describes the measurements of the acoustic and petrophysical properties of two suites of low‐shale sandstone samples from North Sea hydrocarbon reservoirs, under simulated reservoir conditions. The acoustic velocities and quality factors of the samples, saturated with different pore fluids (brine, dead oil and kerosene), were measured at a frequency of about 0.8 MHz and over a range of pressures from 5 MPa to 40 MPa. The compressional‐wave velocity is strongly correlated with the shear‐wave velocity in this suite of rocks. The ratio V_P/V_S varies significantly with change of both pore‐fluid type and differential pressure, confirming the usefulness of this parameter for seismic monitoring of producing reservoirs. The results of quality factor measurements were compared with predictions from Biot‐flow and squirt‐flow loss mechanisms. The results suggested that the dominating loss in these samples is due to squirt‐flow of fluid between the pores of various geometries. The contribution of the Biot‐flow loss mechanism to the total loss is negligible. The compressional‐wave quality factor was shown to be inversely correlated with rock permeability, suggesting the possibility of using attenuation as a permeability indicator tool in low‐shale, high‐porosity sandstone reservoirs.     

Seismoelectric measurements in a porous quartz‐sand sample with anisotropic permeability

Seismoelectric coupling coefficients are difficult to predict theoretically because they depend on a large numbers of rock properties, including porosity, permeability, tortuosity, etc. The dependence of the coupling coefficient on rock properties such as permeability requires experimental data. In this study, we carry out a set of laboratory measurements to determine the dependence of seismoelectric coupling coefficient on permeability. We use both an artificial porous “sandstone” sample, with cracks, built using quartz‐sand and Berea sandstone samples. The artificial sample is a cube with 39% porosity. Its permeability levels are anisotropic: 14.7 D, 13.8 D, and 8.3 D in the x‐, y‐, and z‐directions, respectively. Seismoelectric measurements are performed in a water tank in the frequency range of 20 kHz–90 kHz. A piezoelectric P‐wave source is used to generate an acoustic wave that propagates through the sample from the three different (x, y, and z) directions. The amplitudes of the seismoelectric signal induced by the acoustic waves vary with the direction. The highest signal is in the direction of the highest permeability, and the lowest signal is in the direction of the lowest permeability. Since the porosity of the sample is constant, the results directly show the dependence of seismoelectric coefficients on permeability. Seismoelectric measurements with natural rocks are performed using Berea sandstone 500 and 100 samples. Because the Berea samples are nearly isotropic in permeability, the amplitudes of the seismoelectric signals induced in the different directions are the same within the measurement error. Because the permeability of Berea 500 is higher than that of Berea 100, the amplitude of the seismoelectric signals induced in Berea 500 is higher than those in Berea 100. To determine the relative contributions of porosity and permeability on seismoelectric conversion, we carried out an analysis, using Pride (1994) formulation and Kozeny–Carman relationship; the normalized amplitudes of seismoelectric coupling coefficients in three directions are calculated and compared with the experimental results. The results show that the seismoelectric conversion is related to permeability in the frequency range of measurements. This is an encouraging result since it opens the possibility of determining the permeability of a formation from seismoelectric measurements.     

多孔可变临界孔隙度模型及储层孔隙结构表征

碳酸盐岩、致密砂岩和页岩等储层具有孔隙类型多样、孔隙结构复杂和非均质性强等特征,属于典型的多重孔隙储层,孔隙结构表征是多重孔隙储层预测和流体识别的关键.现有的孔隙结构表征方法大多利用孔隙纵横比或者构建一种新参数来描述孔隙结构.岩石临界孔隙度模型是一种常用的岩石物理模型,具有一定的物理意义和地质含义.本文推导了岩石临界孔隙度与岩石孔隙结构(孔隙纵横比)之间的关系,进而利用极化(形状)因子建立临界孔隙度与弹性参数之间的关系,构建了能够包含多种孔隙类型的多孔可变临界孔隙度模型.利用多孔可变临界孔隙度模型由储层的弹性参数反演不同孔隙类型的体积含量.实验室测量数据和实际测井数据表明,多孔可变临界孔隙度模型能够适用于多重孔隙储层岩石物理建模和孔隙结构表征.     

Compressional- and shear-wave velocities as a function of confining stress in dry sandstones

A laboratory study was carried out to investigate the influence of confining stress on compressional- and shear-wave velocities for a set of rock samples from gas-producing sandstone reservoirs in the Cooper Basin, South Australia. The suite of samples consists of 22 consolidated sublitharenites with helium porosity ranging from 2.6% to 16.6%. We used a pulse-echo technique to measure compressional- and shear-wave velocities on dry samples (cylindrical 4.6 × 2 cm) at room temperature and at elevated confining stress (≤ 60 MPa). Compressional- and shear-wave velocities in samples increase non-linearly with confining stress. A regression equation of the form V = A ? Be^?DP gives a good fit to the measured velocities with improved prediction of velocities at high confining stresses compared with equations suggested by other studies. The predicted microcrack-closure stresses of the samples show values ranging from 70 MPa to 95 MPa and insignificant correlation with porosity, permeability or clay content. There is a positive correlation between change in velocity with core porosity and permeability, but this association is weak and diminishes with increasing confining stress. Experimental results show that pore geometry, grain-contact type, and distribution and location of clay particles may be more significant than total porosity and clay content in describing the stress sensitivity of sandstones at in situ reservoir effective stress. The stress dependence of Cooper Basin sandstones is very large compared with data from other studies. The implication of our study for hydrocarbon exploration is that where the in situ reservoir effective stress is much less than the microcrack-closure stress of the reservoir rocks, the variation of reservoir effective stress could cause significant changes in velocity of the reservoir rocks. The velocity changes induced by effective stress in highly stress-sensitive rocks can be detected at sonic-log and probably surface-seismic frequencies.     

TRANSIENT ELECTROMAGNETIC FIELD OF A POLARIZABLE HALF-SPACE DUE TO VARIOUS CURRENT PULSES *

In this paper, the expressions for the induced voltage function V(T) and mutual impedance function Z_m(T) have been derived for specific range of time T and specific values of induced polarization parameters a for ramp and saw-tooth type of current pulses. The computed results for various cases are also presented. The low values of induced polarization parameter a represent the medium possessing membrane polarization, whereas high values exhibit electrode polarization medium. The method has practical applicability and is best suited for the interpretation of transient electromagnetic fields over a polarizable half-space.     

Electrical Properties of Crustal and Mantle Rocks – A Review of Laboratory Measurements and their Explanation

The electrical properties of rocks and minerals are controlled by thermodynamic parameters like pressure and temperature and by the chemistry of the medium in which the charge carriers move. Four different charge transport processes can be distinguished. Electrolytic conduction in fluid saturated porous rocks depends on petrophysical properties, such as porosity, permeability and connectivity of the pore system, and on chemical parameters of the pore fluid like ion species, its concentration in the pore fluid and temperature. Additionally, electrochemical interactions between water dipoles or ions and the negatively charged mineral surface must be considered. In special geological settings electronic conduction can increase rock conductivities by several orders of magnitude if the highly conducting phases (graphite or ores) form an interconnected network. Electronic and electrolytic conduction depend moderately on pressure and temperature changes, while semiconduction in mineral phases forming the Earth’s mantle strongly depends on temperature and responds less significantly to pressure changes. Olivine exhibits thermally induced semiconduction under upper mantle conditions; if pressure and temperature exceed ~ 14 GPa and 1400 °C, the phase transition olivine into spinel will further enhance the conductivity due to structural changes from orthorhombic into cubic symmetry. The thermodynamic parameters (temperature, pressure) and oxygen fugacity control the formation, number and mobility of charge carriers. The conductivity temperature relation follows an Arrhenius behaviour, while oxygen fugacity controls the oxidation state of iron and thus the number of electrons acting as additional charge carriers. In volcanic areas rock conductivities may be enhanced by the formation of partial melts under the restriction that the molten phase is interconnected. These four charge transport mechanisms must be considered for the interpretation of geophysical field and borehole data. Laboratory data provide a reproducible and reliable database of electrical properties of homogenous mineral phases and heterogenous rock samples. The outcome of geoelectric models can thus be enhanced significantly. This review focuses on a compilation of fairly new advances in experimental laboratory work together with their explanation.     

泥质含量及其分布形式对岩电关系影响的数字岩心仿真(英文)

由于泥质所造成的附加导电现象,泥质含量及其分布形式对电阻率增大系数I和含水饱和度Sw关系具有重要影响,由于岩石物理实验中岩心孔隙结构及其组分构成、分布的微观不可调性,因而泥质分布形式所造成的影响很难通过岩心实验来单独研究。基于数字岩心的格子气自动机方法是一种有效的微观数值模拟方法,本研究利用储层岩心薄片的骨架颗粒尺寸信息资料建立数字岩心模型,结合格子气自动机技术对数字岩心不同饱和流体情况下电的传输特性进行数值模拟研究,揭示了不同泥质含量和泥质分布形式对孔隙介质导电特性非阿尔奇现象产生的影响,建立饱和度指数和泥质含量之间的关系模型,其良好的吻合性表明该方法在岩石物理研究中是一种十分有效的研究方法,而新模型适于在非阿尔奇储层进行准确的饱和度评价。     

Effect of pore geometry on Gassmann fluid substitution

Although there is no assumption of pore geometry in derivation of Gassmann's equation, the pore geometry is in close relation with hygroscopic water content and pore fluid communication between the micropores and the macropores. The hygroscopic water content in common reservoir rocks is small, and its effect on elastic properties is ignored in the Gassmann theory. However, the volume of hygroscopic water can be significant in shaly rocks or rocks made of fine particles; therefore, its effect on the elastic properties may be important. If the pore fluids in microspores cannot reach pressure equilibrium with the macropore system, assumption of the Gassmann theory is violated. Therefore, due to pore structure complexity, there may be a significant part of the pore fluids that do not satisfy the assumption of the Gassmann theory. We recommend that this part of pore fluids be accounted for within the solid rock frame and effective porosity be used in Gassmann's equation for fluid substitution. Integrated study of ultrasonic laboratory measurement data, petrographic data, mercury injection capillary pressure data, and nuclear magnetic resonance T₂ data confirms rationality of using effective porosity for Gassmann fluid substitution. The effective porosity for Gassmann's equation should be frequency dependent. Knowing the pore geometry, if an empirical correlation between frequency and the threshold pore‐throat radius or nuclear magnetic resonance T₂ could be set up, Gassmann's equation can be applicable to data measured at different frequencies. Without information of the pore geometry, the irreducible water saturation can be used to estimate the effective porosity.     

Magnetic and mineralogical characteristics of rocks at the Mezmaiskaya cave Paleolithic site (Northern Caucasus)

The results of magnetic and microprobe studies of the rock sequences in the Mezmaiskaya cave at the unique multilayer Paleolithic site are discussed. The magnetic properties of rocks are analyzed for 17 layers (upwards from layer 4 to layer 1\1) dated from over 73 ka ago to recent times. The rocks of layer 1C (Early Upper Paleolithic, ∼38 ka) are found to have the highest magnetic susceptibility (K) (up to 2500 × 10⁻⁶ SI), which is related to the intensive activity of [Homo sapiens]. The minimum K corresponds to the rocks of layer 2, which is overlain by layer 1D. The sizes of magnetic grains vary throughout the rock section. The largest grains are found in the middle part of the section in the Middle Paleolithic layers 2B3, 2B2, 2B1, 2A, and 2. The superparamagnetic fraction is identified in all layers. This fact supports the view that the cave was open as early as the formation of layer 4. According to the thermomagnetic data on the saturation magnetization and the temperature curves of magnetic susceptibility, magnetite is the main carrier of the rock magnetization; some samples contain iron hydroxides. Samples with iron sulfides (pyrite) are abundant. The study of the hysteresis parameters of rocks showed that the question on whether sulfide-bearing rocks are suitable for reliable paleomagnetic determinations requires further laboratory research into the origin of magnetite in the rocks. The chemical composition of rocks composing layer 2B3 and layers 1D (∼39 ka) and 2B1 (∼45 ka), in which the presence of volcanic ash has been previously established according to the presence of volcanic glass, was determined by detailed microprobe analysis. A wide variety of chemical elements (up to 18 items) was recognized in layers 1D and 2B1. The iron, titan, chrome, manganese content, and concentrations of other components vary from grain to grain. The microprobe analysis of samples from layers 1D and 2B1 revealed a set of magnetic particles with compositions characteristic of volcanic rocks, which supports the ash origin of these layers. Layer 2B3 is established not to be volcanic ash. The results on the volcanic glass in the rocks of layers 1D and 2B1 were published by Golovanova and her colleagues in Current Anthropology in October 2010 [Golovanova et al., 2010].     

Velocities of elastic waves andv P/v S ratios in dry and water-saturated rock samples at high pressure

Elastic wave velocities and the ratiov _P/v _S were studied for dry and initially saturated samples of carbonate and crystalline rocks at pressures to 2 kbar. In initially saturated samplesv _P increases in crystalline rock, whereas in sedimentary rock it may either increase or decrease with increasing pressure. Under the same conditionsv _S remains approximately constant in crystalline rocks and decreases in sedimentary samples. The ratiov _P/v _S as a function of pressure increases in dry rocks and decreases in saturated ones. Saturated samples always have higherv _P/v _S ratios than dry ones. It appears that the geometry of the pore space influences the acoustic properties of saturated crystalline rocks.     

Nuclear magnetic resonance <Emphasis Type="Italic">T</Emphasis><Subscript>2</Subscript> spectrum: multifractal characteristics and pore structure evaluation

Pore structure characteristics are important to oil and gas exploration in complex low-permeability reservoirs. Using multifractal theory and nuclear magnetic resonance (NMR), we studied the pore structure of low-permeability sandstone rocks from the 4th Member (E_S4) of the Shahejie Formation in the south slope of the Dongying Sag. We used the existing pore structure data from petrophysics, core slices, and mercury injection tests to classify the pore structure into three categories and five subcategories. Then, the T₂ spectra of samples with different pore structures were interpolated, and the one- and three-dimensional fractal dimensions and the multifractal spectrum were obtained. Parameters α (intensity of singularity) and f (α) (density of distribution) were extracted from the multifractal spectra. The differences in the three fractal dimensions suggest that the pore structure types correlate with α and f (α). The results calculated based on the multifractal spectrum is consistent with that of the core slices and mercury injection. Finally, the proposed method was applied to an actual logging profile to evaluate the pore structure of low-permeability sandstone reservoirs.     

Reservoir parameter classification of a Miocene formation using a fractal approach to well logging,porosimetry and nuclear magnetic resonance

A methodology for rock classification is presented that considers lithology and reservoir parameters on the basis of a combined fractal analysis of well logs and mercury porosimetry results with nuclear magnetic resonance outcomes. A sandy‐shaly thinly‐bedded Miocene gas bearing formation in the Carpathian Foredeep is investigated. Fractal correlation dimensions D₂ calculated for standard logs are used to distinguish sandstone as the most homogeneous lithological group with the highest porosity. The fractal analysis also confirmed observations of gas accumulations in sandstone, shaly sandstone and sandy claystone as rocks of high porosity. The results of two laboratory methods are combined to improve reservoir properties assessment and evaluate movable media in pore space; this technique was based on the similarity of mercury porosimetry results plotted as cumulative intrusion volume versus pressure or pore diameter and also the curves of cumulative porosity and transverse relaxation time distributions using nuclear magnetic resonance. Close values of porosity from logs, recorded in situ and in laboratory measurements, provide the link between fractal analysis and porosimetry and nuclear magnetic resonance measurements.     

Effect of Stress on the Magnetic Memory of Induced Magnetic Anisotropy of Rocks and Its Mathematical Model

When re-heated to temperatures below the Curie temperature and subsequently cooled in a constant magnetic field (H _T), rock samples which contain magnetic minerals can acquire an induced magnetic anisotropy (IMA). As the result of acquiring the IMA, a constriction develops in the hysteresis loop of the magnetization of these rocks at the values of the magnetizing field close or equal to the HT. Thus the IMA is capable of retaining the information on the palaeointensity of the geomagnetic field, i.e., if IMA was created in a rock in the geomagnetic field in a past geological epoch, it preserves the information on the intensity of that field. Investigations have shown, that when IMA is created in a rock under external stress, the stress has an impact on the magnetic memory. Here we also deal with the issue of how stress affects the magnetic memory of IMA. A mathematical model for the effect of stress on magnetic memory phenomena related to induced magnetic anisotropy in rocks containing multidomain magnetite and titanomagnetite grains is proposed herewith. The effect of temperature on the magnetic memory of rocks is discussed also.     

Permeability Prediction for Low Porosity Rocks by Mobile NMR

Estimating permeability from NMR well logs or mobile NMR core scanner data is an attractive method as the measurements can be performed directly in the formation or on fresh cores right after drilling. Furthermore, the method is fast and non-destructive. Compared to T ₁ relaxation times, commonly measured T ₂ distributions are influenced by external and internal magnetic field gradients. We performed two-dimensional T ₁ and T ₂ relaxation experiments on samples of Rhaetian sandstone, a rock with low porosity and small pore radii, using a mobile NMR core scanner which operates within a nearly homogeneous static magnetic field. Because small pore sizes are associated with high internal magnetic field gradients, standard methods from NMR logging in the oil industry cannot be applied for accurate permeability prediction. Therefore, a new model theory was developed, which describes the pore radius dependence of the surface relaxivity ρ ₂ by both an analytical and a more practical empirical equation. Using corrected ρ ₂ values, permeability can be predicted accurately from the logarithmic mean of the T ₂ distribution and the physically based Kozeny-Carman equation. Additional core plug measurements of structural parameters such as porosity, permeability, specific inner surface area and pore radius distributions supported the NMR results.     

Estimates of possible diffusional effects in trace element separations

Igneous material dredged from the Rio Grande rise, South Atlantic Ocean, includes basaltic rocks, some having mafic nodules and megacrysts, and volcanic breccias composed largely of basaltic fragments. These samples represent the only volcanic rocks recovered from this aseismic rise. Bulk compositions show alkalic basalt, trachybasalt, and trachyandesite; the rock types are similar to those of nearby Tristan da Cunha, Gough, and the Walvis ridge. Microprobe analyses show basaltic groundmass to have olivine, Fo₈₅, pyroxene, Fs₁₃Wo₄₆, feldspar, An₇₁, plus interstitial alkali feldspar. Mafic nodules and megacrysts have olivine, Fo_86–90 and pyroxene Fs_6–7.5Wo_45–46; Al₂O₃ 2.5–4 wt.%.The Rio Grande rise rocks have compositional characteristics of an alkalic basaltic suite, and not of mid-ocean ridge tholeiite. Based on mineral compositions, nodules and megacrysts in basalt are interpreted as cognate inclusions. Because oceanic alkalic basaltic rocks are almost invariably associated with islands and seamounts, the Rio Grande rise probably represents a series of alkalic-basalt islands that formed and eventually subsided during rifting of the South Atlantic; the dredged volcanic breccias are probably slump deposits from those volcanoes. This interpretation lends support to the Rio Grande rise having formed at a hot spot, but the possibility of alkalic rocks having formed along fracture zones should not be discounted.     

Petrogenesis and tectonic setting of bimodal volcanism in the Sakoli Mobile Belt, Central Indian shield

The Sakoli Mobile Belt comprises bimodal volcanic rocks that include metabasalt, rhyolite, tuffs, and epiclastic rocks with metapelites, quartzite, arkose, conglomerate, and banded iron formation (BIF). Mafic volcanic rocks are tholeiitic to quartz‐tholeiitic with normative quartz and hypersthene. SiO₂ shows a large compositional gap between the basic and acidic volcanics, depicting their bimodal nature. Both the volcanics have distinct geochemical trends but display some similarity in terms of enriched light rare earth element–large ion lithophile element characteristics with positive anomalies for U, Pb, and Th and distinct negative anomalies for Nb, P, and Ti. These characteristics are typical of continental rift volcanism. Both the volcanic rocks show strong negative Sr and Eu anomalies indicating fractionation of plagioclases and K‐feldspars, respectively. The high Fe/Mg ratios for the basic rocks indicate their evolved nature. Whole rock Sm–Nd isochrons for the acidic volcanic rocks indicate an age of crystallization for these volcanic rocks at about 1675 ± 180 Ma (initial ¹⁴³Nd/¹⁴⁴Nd = 0.51017 ± 0.00017, mean square weighted deviate [MSWD] = 1.6). The εNd_t (t = 2000 Ma) varies between ?0.19 and +2.22 for the basic volcanic rock and between ?2.85 and ?4.29 for the acidic volcanic rocks. Depleted mantle model ages vary from 2000 to 2275 Ma for the basic and from 2426 to 2777 Ma for the acidic volcanic rocks, respectively. These model ages indicate that protoliths for the acidic volcanic rocks probably had a much longer crustal residence time. Predominantly basaltic magma erupted during the deposition of the Dhabetekri Formation and part of it pooled at crustal or shallower subcrustal levels that probably triggered partial melting to generate the acidic magma. The influence of basic magma on the genesis of acidic magma is indicated by the higher Ni and Cr abundance at the observed silica levels of the acidic magma. A subsequent pulse of basic magma, which became crustally contaminated, erupted as minor component along with the dominantly acidic volcanics during the deposition of the Bhiwapur Formation.     

Magnetic properties of high-Ti basaltic rocks from the Krušné hory/Erzgebirge MTS. (Bohemia/Saxony), and their relation to mineral chemistry

This study provides new thermomagnetic and petrographic data on specific basaltic rock association from the broader vicinity of the Lou?ná-Oberwiesenthal volcanic centre, western Bohemia/Saxony. Two types of volcanic rocks were recognized there: (i) high-Ti types (3.5–5.2 wt% TiO₂) represented by (mela)nephelinite s.s., and sporadically present (ii) medium-Ti types (2.5–3.5 wt% TiO₂) of olivine nephelinite, nepheline basanite and phonotephrite compositions. In order to examine the rock-magnetic behaviour, they were studied for their variations in the Curie temperature (T_C) and field-dependent susceptibility, spinel group minerals, chemistry and petrology. Magnetic susceptibility of ulvöspinel-rich titanomagnetite, as a dominant magnetic carrier, depends on the amplitude of measured magnetic field, whereas pure magnetite is field-independent. Field dependence parameter ^kHD of the studied basaltic rocks ranges from 0.8 to 18.7%, TiO₂ contents in titanomagnetite range from 12.7 to 20.1 wt.%. TiO₂ content in titanomagnetite does not correlate with whole-rock TiO₂ content (2.8 to 5.6 wt.%). The content of substituted titanium in the sublattice of magnetite is also sensitively reflected in the Curie temperature, ranging from 200 to 580°C. The spinel group minerals are designated as titanomagnetite with the dominance of ulvöspinel, magnetite and magnesioferrite components, or titanomagnetite with the magnetite, ulvöspinel and magnesioferrite components. Only two samples are characterized by a significant presence of Cr-spinel and magnesiochromite components forming cores of titanomagnetites. The titanomagnetite-bearing rocks in the studied area, likewise the low- to medium-Ti basaltic rocks from the ?eské st?edoho?í Mts., provide similar thermomagnetic curves.     

基于等效介质理论的孔隙纵横比分布反演

孔隙纵横比是描述多孔岩石微观孔隙结构特征的重要参数,目前用于获取岩石完整孔隙纵横比分布的经典模型为David-Zimmerman(D-Z)孔隙结构模型,该模型假设岩石由固体矿物基质、一组纵横比相等的硬孔隙以及多组纵横比不等的微裂隙构成,并认为固体矿物基质和硬孔隙均不受压力影响,在此基础上,利用超声纵横波速度的压力依赖性反演岩石硬孔隙和各组微裂隙的孔隙纵横比及孔隙度.该方法的关键点在于以累积裂隙密度为桥梁,借助等效介质理论建立了岩石弹性模量和孔隙纵横比之间的内在联系.但在D-Z模型中,多重孔隙岩石累积裂隙密度的计算直接由单重孔隙裂隙密度公式实现,这种近似导致该模型在许多情况下难以获得良好的反演精度.为了完善经典D-Z模型,本文提出了一种基于虚拟降压的孔隙纵横比分布反演策略,通过多个假想降压过程实现累积裂隙密度的准确计算,并将基于DEM和MT的经典D-Z模型推广到KT和SCA中,结合四种等效介质理论建立了一套完整的反演流程.采用一系列砂岩和碳酸盐岩样品,测试了反演流程在实际岩芯孔隙纵横比提取中的应用效果,研究结果表明:与D-Z模型相比,本文方法的模拟结果与实际数据吻合更好,并同时适用于砂岩和碳酸盐岩;此外,通过分析四种等效介质理论的模拟结果发现,本文方法并不十分依赖于等效介质理论的选择,这些理论获得的孔隙结构参数随压力的变化趋势基本一致,数值上仅存在略微差异,且这种差异随着压力的增大逐渐消失.本文方法是经典D-Z孔隙结构模型的重要补充,对岩石孔隙结构表征、流体饱和岩石速度预测以及孔间喷射流效应的模拟具有十分重要的意义.     

Flow and fracture maps for basaltic rock deformation at high temperatures

Understanding how the strength of basaltic rock varies with the extrinsic conditions of stress state, pressure and temperature, and the intrinsic rock physical properties is fundamental to understanding the dynamics of volcanic systems. In particular it is essential to understand how rock strength at high temperatures is limited by fracture. We have collated and analysed laboratory data for basaltic rocks from over 500 rock deformation experiments and plotted these on principal stress failure maps. We have fitted an empirical flow law (Norton’s law) and a theoretical fracture criterion to these data. The principal stress failure map is a graphical representation of ductile and brittle experimental data together with flow and fracture envelopes under varying strain rate, temperature and pressure. We have used these maps to re-interpret the ductile–brittle transition in basaltic rocks at high temperatures and show, conceptually, how these failure maps can be applied to volcanic systems, using lava flows as an example.     

Early Jurassic intra‐oceanic arc system of the Neotethys Ocean: Constraints from andesites in the Gangdese magmatic belt,south Tibet

The Gangdese magmatic belt is located in the southern margin of the Lhasa terrane, south Tibet. Here zircon U–Pb ages and Hf isotopic data, as well as whole‐rock geochemistry and Sr–Nd isotopes on andesites from the Bima Formation with a view to evaluating the history of the Gangdese magmatism and the evolution of the Neotethys Ocean. Zircon U–Pb dating yields an age of ca 170 Ma from six samples, representing the eruptive time of these volcanic rocks. Zircon Hf isotopes show highly positive ε_Hf(t) values of +13 to +16 with a mean of +15.2. Whole‐rock geochemical and Sr–Nd isotopic results suggest that the magma source of these andesites was controlled by partial melting of a depleted mantle source with addition of continental‐derived sediments, similar to those in the southern arcs of the Lesser Antilles arc belt. In combination with published data, the volcanic rocks of the Bima Formation are proposed to have been generated in an intra‐oceanic arc system, closely associated with northward subduction of the Neotethyan oceanic lithosphere.