Dependence on pressure of conduction by hopping of small polarons in minerals of the Earth's lower mantle

Electrical conductivity of the lower mantle-like assemblage (Mg,Fe)SiO₃ perovskite-(Mg,Fe)O magnesiowüstite is usually analyzed using the quasi-chemical Arrhenian approach of diffusion. The conductivity of this assemblage has often been attributed to hopping of small polarons, because of the low value of the activation energy and the small negative activation volume. However, the solid-state physics approach can provide more arguments, for or against conduction by polarons. We have tried to bridge the gap between the two approaches and identify the physical quantities entering the phenomenological activation parameters. In particular, we have investigated the pressure dependence of the activation energy, and the physical meaning of the activation volume. Hopping is controlled by the binding energy of the polaron and by the value of the exchange integral, which increases with pressure causing the observed decrease of the activation energy. From the physical theory and the results of experiments at pressures up to 40 GPa and temperatures up to 400 ^∘C, we have estimated the values of parameters characteristic of polarons: radius, mobility, time between jumps and adiabaticity. These values are compatible with conduction by small adiabatic polarons. The consequences for extrapolations to lower mantle conditions of the presence of a temperature dependent preexponential term in the expression for conductivity have been examined. It was found that the extrapolations are not significantly different from those using the Arrhenius equation. Received: 5 November 1998 / Revised, accepted: 4 May 1999     

Deep-Sea Sediment Compression Curves: Some Controlling Factors, Spurious Overconsolidation, Predictions, and Geophysical Reproduction

Oedometer tests have been carried out on 70 undisturbed surficial clays (at approximately 250 mm below the mudline), mostly collected by free-fall corers from sites widely scattered throughout the deep-sea North Atlantic. Acoustic measurements were also made, initially on contiguous samples and ultimately on the same sample using a geophysically instrumented oedometer which also collected electrical resistivity data. Apart from those quiescent areas below the carbonate compensation depth, such as north of the West Indies where very fine clays exist, most of the samples are silty clays whose geotechnical-geophysical properties are dependent on the type of clay minerals present (and their ability to take in moisture), the sand-size fraction, and the quantity of carbonate present. Thus the pure clays have high compressibilities which decrease on the addition of coarse particles, while the converse is true for the acoustic parameters, these increasing with the sand fraction. Using the notion of the intrinsic compression line for all samples, and comparison to it of the measured compression curves, it is clear that, contrary to some previously held ideas, most deep-sea clays are normally consolidated; the addition of carbonate has the effect of creating an open, stronger sediment skeleton. Interestingly, where information is available, the variation with depth of a sample's acoustic velocity follows the void ratio pressure relationship of the compression curve. This allows the construction of an in-situ sediment compression curve using the in-situ geophysical observations.     

面向目标自适应有限元法的带地形三维大地电磁各向异性正演模拟

传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础.     

基于MPI并行算法的电阻率法多种装置数据的三维联合反演

由于地表电阻率法受到浅地表局部异常体的干扰,反演精度受到影响,井中装置数据资料参与反演虽然可以减小浅地表局部异常体的干扰,但是由于钻井位置的局限性,数据量得不到保障,也会导致反演精度降低.为此,本文开发了一套结合地表、地-井、井-地和井-井装置数据的三维联合反演算法.首先,利用有限差分法实现正演模拟,采用非线性共轭梯度法（NLCG）恢复电阻率结构;其次,调用Message Passing Interface（MPI）函数库加速正演模拟和灵敏度矩阵运算,当开辟12个进程时,反演程序获得最大加速比4.51;最后,通过合成数据和实测数据算例证明该反演算法的有效性,也证实了多种装置组合数据体反演结果明显优于单一地表装置数据体反演结果.     

含碳结构对龙门山断层带电导率影响的实验探索

碳是影响岩石电导率大小的一个重要因素,可能是造成龙门山断层带电导率异常的重要原因之一.为了研究不同的碳含量、矿物颗粒粒径与碳晶体结构对断层带电导率的影响,在干燥、常温、0.2~300 MPa的压力条件下实验研究了人工模拟断层泥样品(石英粉末与含碳粉末混合的样品,简称模拟样品)和采自映秀-北川断层八角庙剖面的天然断层岩样品(简称天然样品)的电导率.实验结果显示,当模拟样品中的含碳粉末连通时,电导率与碳体积率的关系符合逾渗理论模型;而含碳粉末未连通时,电导率随总孔隙度降低而指数性升高.同时模拟样品的电导率也随石英颗粒粒径的变化而发生改变.相比于模拟样品中的含碳粉末主要分布于石英颗粒支撑的孔隙中,天然样品中的碳则主要以碳膜的形式赋存在颗粒边缘,导致碳体积率相同的条件下,模拟样品的电导率小于天然样品.此外,天然样品的电导率(9×10~(-4)S·m~(-1))也要小于野外大地电磁探测的结果(0.03~0.1 S·m~(-1)).在今后的实验中还需要考虑在动态摩擦条件下对含有完整含碳结构的天然样品进行电导率的实验研究.     

Uncertainty analysis and probabilistic segmentation of electrical resistivity images: the 2D inverse problem

In this paper, we present the uncertainty analysis of the 2D electrical tomography inverse problem using model reduction and performing the sampling via an explorative member of the Particle Swarm Optimization family, called the Regressive‐Regressive Particle Swarm Optimization. The procedure begins with a local inversion to find a good resistivity model located in the nonlinear equivalence region of the set of plausible solutions. The dimension of this geophysical model is then reduced using spectral decomposition, and the uncertainty space is explored via Particle Swarm Optimization. Using this approach, we show that it is possible to sample the uncertainty space of the electrical tomography inverse problem. We illustrate this methodology with the application to a synthetic and a real dataset coming from a karstic geological set‐up. By computing the uncertainty of the inverse solution, it is possible to perform the segmentation of the resistivity images issued from inversion. This segmentation is based on the set of equivalent models that have been sampled, and makes it possible to answer geophysical questions in a probabilistic way, performing risk analysis.     

Effects of fluid displacement pattern on complex electrical impedance in Berea sandstone over frequency range 104–106 Hz

To better understand the effect of fluid distribution on the electric response of rocks saturated with oil and brine, we conducted experimental studies on the complex electrical impedance in a Berea sandstone, together with in situ acquisitions of oil distribution images employing a high‐resolution medical X‐ray computed tomography. We performed two tests of brine displacement by oil under high (10 MPa) and low (5 MPa) pressures, which were accompanied by fingering and stable displacement patterns, respectively. The measured complex impedance data were fitted to the Cole model to obtain the resistance, capacitance, peak frequency of the imaginary impedance, and the exponent α of the rock–fluid system. With increasing oil saturation, the resistance showed an increasing trend, whereas the other three parameters decreased. The fingering displacement exhibited lower resistance and capacitance than the stable displacement. The analysis of the resistance changes using a simple parallel connection model indicates that there are more components of residual brine in parallel connections in the fingering pattern than in the stable displacement pattern at the same saturation. We also interpreted the normalised changes in the capacitance (or apparent dielectric constant) with respect to the oil saturation via an analysis of the shape factor of fluid distribution based on the Maxwell–Wagner–Brugermann–Hanai model. The changes in the shape factor suggest that the pinch‐off of the brine in parallel connection by the oil is a dominant mechanism reducing the capacitance. In the stable displacement, most of the connections in the brine phase are immediately pinched off by oil displacement front at a local oil saturation of 65%. Conversely, in the fingering displacement, there is a transition from the bulk or layered brine to the pinched‐off at a local oil saturation below 60%. The analyses indicate that the difference in the fluid distribution under different fluid conditions is responsible for the non‐Archie behaviour.     

一维连续导电剖面的反演问题

基于Schlumberger装置测得的视电阻率,在Coen理论的基础上对一维连续直流导电剖面的反演问题进行研究。分别给出了导电剖面是常函数、线性函数、幂函数、指数函数、分段函数情况下的精确数据和噪音数据反演结果。利用一个野外场地反演的例子,把利用该方法进行反演得到的连续模型和Inman得到的层状模型进行了相应的比较,发现该方法在处理实际场地数据也是有效的。     

Hydrogeophysical modelling of Hisarcik (Kütahya) geothermal field,western Turkey

Western Anatolia hosts many low-to-moderate and high-temperature geothermal sources in which active faults play a dominant role to control the recharge and the discharge of geothermal fluid. In this study, we used the two-dimensional geoelectric structure of Kütahya Hisarcık geothermal field, and created a conceptual hydrogeophysical model that includes faults, real topographical variations and geological units. The temperature distribution and fluid flow pattern are also investigated. The depth extension of Hisarcık Fault, electrical basement and low resistivity anomalies related to the presence of geothermal fluid are determined by using resistivity studies in the area. Numerical simulations suggest that Hisarcık fault functioning as a fluid conduit primarily enables hot fluid to be transported from depth to the surface. It is shown that the locations of predicted outflow vents coincide with those of hot springs in the area.     

Non-invasive lab test in the monitoring of vadose zone contaminated by light non-aqueous phase liquid

We discuss a strategy capable of a quantitative long-term monitoring of water saturation and volume variation of light non-aqueous phase liquids in the soil. The goal was reached monitoring a controlled sand cell contaminated with classical gasoline over 124 days, using geophysical methods such as electrical resistivity tomography, induced polarization and ground penetrating radar. We show that empirical relations, linking the water saturation to the physical parameters measured as resistivity from electrical resistivity tomography and travel time from georadar with advanced processing, are good tools for this purpose. The consistence of the proposed process is validated by both good overlap of results carried out from electrical resistivity tomography and georadar and theoretical models simulating the actual experiment.