高压下单晶橄榄石的电导率

田浩然; 徐良旭; 李娜娜; 张茜; 林俊孚; 刘锦

doi:10.11858/gywlxb.20190775

高压下单晶橄榄石的电导率

doi: 10.11858/gywlxb.20190775

田浩然^1,,
徐良旭¹,
李娜娜¹,
张茜¹,
林俊孚²,
刘锦^1,*, ,

1.
北京高压科学研究中心，北京　100193
2.
德克萨斯大学奥斯汀分校，美国德克萨斯 78701

基金项目: 国家自然科学基金（U1530402）

详细信息

作者简介:
田浩然（1993－），男，硕士研究生，主要从事矿物物理学研究. E-mail：haoran.tian@hpstar.ac.cn

通讯作者:
刘　锦（1984－），男，博士，研究员，主要从事地球深部结构与物质组成研究. E-mail：jin.liu@hpstar.ac.cn

中图分类号: O732; P574.1
计量
- 文章访问数: 288
- HTML全文浏览量: 126
- PDF下载量: 6
出版历程
- 收稿日期: 2019-05-13
- 修回日期: 2019-06-18
- 刊出日期: 2019-12-01

High-Pressure Electrical Conductivity of Single-Crystal Olivine

TIAN Haoran^1
,,
XU Liangxu¹,
LI Nana¹,
ZHANG Qian¹,
LIN Junfu²,
LIU Jin^{1,*
, ,}

1.
Center for High Pressure Science & Technology Advanced Research, Beijing 100193, China
2.
The University of Texas at Austin, Texas 78701, USA

摘要

摘要: 以圣卡洛斯（San Carlos）单晶橄榄石为研究对象，结合交流阻抗谱和金刚石对顶砧（DAC）技术，在300 K、0～19 GPa条件下对其电导率的各向异性进行系统研究。压力标定根据红宝石荧光谱线的漂移以及硅油的拉曼光谱。实验结果表明：在300 K、0～19 GPa条件下，橄榄石[100]方向上的电导率最大，从3.8×10^–8 S/m增加到9.0×10^–8 S/m，[010]与[001]方向上的电导率接近，约为[100]方向电导率的1/2～1/3；橄榄石电导率随着压力线性增加，其中[100]方向的电导率随着压力变化的斜率最大。在室温条件下，橄榄石主要的导电机制是小极化子导电，且具有负的活化体积。研究结果表明，在含水量较低的上地幔区域，随着深度增加，压力效应可能导致电导率横向和纵向的不均一性增强。
- 橄榄石 /
- 复阻抗谱 /
- 金刚石压腔 /
- 电导率 /
- 各向异性
Abstract: The electrical conductivity of single-crystal San Carlos olivine was measured up to 19 GPa at room temperature in a diamond-anvil cell, coupled with a complex impedance spectroscopy. The pressure was determined by in-situ ruby luminescence and Raman shift of silicone fluid. We found that the electrical conductivity along [100] is largest, increasing approximately from 3.8×10^–8 S/m at 0 GPa to 9.0×10^–8 S/m at 18 GPa at room temperature. The conductivity along [010] is comparable to that of [001], approximately as 1/2 to 1/3 as that of [100]. Furthermore, the conductivity linearly increases with the pressure, while it changes faster with the pressure along [100] than that of [010] and [001]. At room temperature, the charge transport mechanism of olivine is dominant from the Fe²⁺–Fe³⁺ (small polarons) with a negative activation volume. The present results suggest that the pressure effect could lead to larger lateral and vertical heterogeneity in electrical conduction for a dry upper mantle.
- olivine /
- complex impedance spectroscopy /
- diamond anvil cell /
- electrical conductivity /
- anisotropy

HTML全文

近年来随着液化天然气(LNG，Liquefied Natural Gas)产业在全球迅速发展，天然气的液化技术和设备也在不断发展完善、日渐成熟。我国LNG领域内的相关研究起步较晚，许多技术远远落后发达国家水平，在天然气液化工艺及装置的生产等方面缺乏自主产权。因此，开展天然气液化工艺及装置的研究，对于实现液化装置的国产化、高效化有十分重要的意义^[1-3]。

超声速旋流分离技术是一种新兴的天然气加工处理技术，被较为广泛地用于天然气脱水、脱重烃、脱酸等方面，近年来开始逐渐应用于天然气液化方面^[4-7]。天然气超声速液化的原理是：高压天然气混合物在Laval喷管内达到一定的温度、压力条件，开始凝结成核，最终凝结成液滴，在后续工艺中进行进一步气液分离。与传统的天然气液化技术相比，具有结构工艺简单、支持无人操作(适用于海底天然气处理)、对水合物抑制剂依赖性小、投资和运行成本低等优势^[8-9]。

为了探究入口复杂多变的压力条件对天然气超声速液化特性的影响，对甲烷-乙烷气体混合物的超声速凝结流动特性进行研究，在凝结成核与生长理论的基础上建立了适用于甲烷-乙烷双可凝气体混合物的凝结流动数学模型，重点研究了入口压力对天然气混合物在Laval喷管内主要流动与凝结参数的影响规律。

1. Laval喷管结构设计

Laval喷管结构主要包括入口段、渐缩段、喉部及扩张段4部分^[10-11]。各部分参数如表 1所示，L₀为入口长度，r₁为渐缩段入口半径，L₁为渐缩段长度，r_cr为喉部截面半径，L₂为渐扩段长度，r₂为渐扩段出口半径。为尽量减小流场涡流的影响，渐缩段采用双三次曲线设计，喉部采用一段平缓光滑的圆弧作为过渡曲线，渐扩段采用等膨胀率设计，膨胀率取为10 000 s^-1。考虑到实验加工方便，保证曲面的精度，且能够更加直观地观察Laval喷管内部的流场分布情况，所设计Laval喷管截面采用矩形截面，三维结构如图 1所示。

表 1 Laval喷管各部分参数

Table 1. Parameters of Laval nozzle

L₀/mm	r₁/mm	L₁/mm	r_cr/mm	L₂/mm	r₂/mm
50.00	20.00	56.01	2.50	71.28	6.15

下载: 导出CSV

| 显示表格

图 1 Laval喷管三维结构

Figure 1. Structure of Laval nozzle

下载: 全尺寸图片幻灯片

2. 超声速凝结流动数学模型及计算方法

2.1 数学模型

采用欧拉双流体模型开展数值计算，控制方程主要包括气相流动方程和液相流动方程。在无滑移假设及欧拉双流体模型的前提下分别建立气相及液相流动控制方程组，液滴数目守恒方程及液滴半径、数目、湿度关系式分别添加到对应源相方程中^[12-14]。

气相流动控制方程组

$\frac{{\partial {\rho _{\rm{v}}}}}{{\partial t}} + \frac{\partial }{{\partial {x_j}}}\left( {{\rho _{\rm{v}}}{u_j}} \right) = {S_{\rm{m}}}$

(1)

$\frac{{\partial {\rho _{{\rm{c2}}}}}}{{\partial t}} + \frac{\partial }{{\partial {x_j}}}\left( {{\rho _{{\rm{c2}}}}{u_j}} \right) = {S_{{\rm{m}}, {\rm{c}}2}}$

(2)

$\begin{array}{l} \frac{\partial }{{\partial t}}\left( {{\rho _{\rm{v}}}{u_i}} \right) + \frac{\partial }{{\partial {x_j}}}\left( {{\rho _{\rm{v}}}{u_j}{u_i}} \right) = \\ - \frac{{\partial {p_{\rm{v}}}}}{{\partial {x_i}}} + \frac{\partial }{{\partial {x_j}}}\left[{\mu \left( {\frac{{\partial {u_j}}}{{\partial {x_i}}} + \frac{{\partial {u_i}}}{{\partial {x_j}}}-\frac{2}{3}{\delta _{{\rm{ij}}}}\frac{{\partial {u_j}}}{{\partial {x_j}}}} \right)} \right] + \frac{\partial }{{\partial {x_j}}}\left( { -{\rho _v}\overline {{{u'}_i}{{u'}_j}} } \right) + {S_{\rm{u}}} \end{array}$

(3)

$\frac{\partial }{{\partial t}}\left( {{\rho _{\rm{v}}}E} \right) + \frac{\partial }{{\partial {x_j}}}\left( {{\rho _{\rm{v}}}{u_j}E + {u_j}{p_{\rm{v}}}} \right) = \frac{\partial }{{\partial {x_j}}}\left( {{k_{{\rm{eff}}}}\frac{{\partial T}}{{\partial {x_j}}} + {u_i}{\tau _{{\rm{eff}}}}} \right) + {S_{\rm{h}}}$

(4)

液相流动控制方程组

$\frac{\partial }{{\partial t}}\left( {\rho Y} \right) + \frac{\partial }{{\partial {x_j}}}\left( {\rho {u_j}Y} \right) = {S_{\rm{Y}}}$

(5)

$\frac{{\partial \rho N}}{{\partial t}} + \frac{\partial }{{\partial {x_j}}}\left( {\rho N{u_j}} \right) = J$

(6)

${r_{\rm{d}}} = \sqrt[3]{{3Y/\left( {4\pi {\rho _{\rm{L}}}N} \right)}}$

(7)

式中:u_i、u_j为时均速度分量，m/s；ρ_v为气相密度，kg/m³；ρ为气液混合相密度，kg/m³；p_v为时均压力，Pa；μ为黏度，kg/(m·s)；δ_ij为Kronecker delta数；E为总能量，J/kg；k_eff为有效导热系数，W/(m·K)；τ_eff为有效应力张量，无量纲；Y为液相质量分数，无量纲；r_d为液滴半径，m；dr_d/dt为液滴生长速率，m/s；N为液滴数目，kg^-1。

成核模型采用文献[15-16]中提出的双组分气体自发凝结成核模型修正方法。液滴生长过程采用Gyarmathy液滴生长模型，模型中液滴与气体间的传热系数^[17-18]为

${k_{\rm{r}}} = \frac{{{\lambda _{\rm{v}}}}}{{{r_{\rm{d}}}}}\frac{1}{{1 + \frac{{2\sqrt {8\pi } }}{{1.5{P_{\rm{v}}}}}\frac{\gamma }{{1 + \gamma }}{K_{\rm{n}}}}}$

(8)

依据传热、传质过程，可推导得到液滴生长速率计算模型

$\frac{{{\rm{d}}{r_{\rm{d}}}}}{{{\rm{d}}t}} = \frac{{{\lambda _{\rm{v}}}}}{{{\rho _L}{h_{{\rm{LV}}}}}}\frac{{\left( {{T_{\rm{s}}}-T} \right)\left( {1-\frac{{{r_{\rm{c}}}}}{{{r_{\rm{d}}}}}} \right)}}{{{r_{\rm{d}}}\left( {1 + \frac{{2\sqrt {8\pi } }}{{1.5{P_{\rm{v}}}}}\frac{\gamma }{{1 + \gamma }}{K_{\rm{n}}}} \right)}}$

(9)

式中:λ_v为气体导热系数，W/(m·K)；P_v为气体Prandtl数；γ为气体比热比；h_LV为凝结潜热，J/kg；T_s为气体压力对应的饱和温度，K；K_n表示Kundsen数。由于双组分气体不存在压力对应的饱和温度T_s这一概念，将双组分相图中露点线类比于单组分中饱和曲线。

针对气体状态方程的选择，由于低温气体已偏离理想气体假设，本研究采用了NIST真实气体模型进行计算。

2.2 湍流模型

湍流发生时会导致流体之间相互交换动量、能量，也会造成浓度的改变。本研究建模时忽略相间速度的滑移，即液滴产生不影响湍流，因此只考虑气相的湍流方程。FLUENT中提供了以下几种湍流模型：S-A模型、标准k-ε模型、RNG k-ε模型、Realizable k-ε模型、k-ω模型以及雷诺应力模型。S-A模型主要应用流动分离区附近模拟，标准k-ε模型、RNG k-ε模型一般用于各向同性的均匀湍流，k-ω模型可用于带压梯度的流动模拟和跨声速激波模拟，雷诺应力模型主要用于龙卷风、燃烧室等强烈旋转流动的模拟。对Laval喷管内跨声速流动，采用k-ω模型可以获得较为理想的计算精度和计算速度，故本研究采用该模型进行数值计算。

2.3 计算方法

气体在Laval喷管中的流动属于高速可压缩流动，采用密度基进行求解，流动控制方程组、湍流动能方程、湍流耗散率方程均采用二阶迎风格式进行离散。

根据双组分气体在Laval喷管内的高速可压缩的流动特性，入口和出口边界设置为压力入口边界和压力出口边界条件，对于气体在Laval喷管内的超声速流动，由于所有的流动参数都可从Laval喷管内部外推得到，故在出口处不进行相应设置，固体壁面边界设置为无滑移、无渗流、绝热边界条件。

在数学模型中，由于气相方程添加了源相方程，液相方程定义了标量以及引入的真实气体方程，这些仅靠在FLUENT自带的模型和材料物性无法满足要求，需要编写相应的用户自定义函数(UDF)。本研究编写UDF时，分别定义DEFINE AJUST、DEFINE SOURCE和DEFINE PROPERTY 3个宏函数。DEFINE AJUST宏用来定义过饱和度、过冷度、成核速率、液滴生长率、液滴半径、液滴质量以及液滴表面张力等参数，DEFINE SOURCE宏用来定义控制方程中的质量、动量和能量源相，DEFINE PROPERTY用来定义数值计算中用到的真实气体的热力学参数如黏度系数、导热系数等。

3. 实验验证

为验证所建立的双组分气体凝结数学模型及数值计算方法的准确性，采用本研究所设计的Laval喷管结构，在中国石油大学(华东)超声速气体凝结流动实验系统开展了水-乙醇双可凝组分气体凝结相变实验研究。实验条件为：Laval喷管入口压力0.586 MPa，入口温度288.05 K，气体湿度98.1%，水与乙醇摩尔体积比84:16，气体流量为323.78标方每小时，实验测得的Laval喷管沿程压力分布如图 2所示，可以看出，压力分布实验结果与数值计算结果吻合较好，说明本研究所建立的双组分气体超声速凝结流动特性数学模型及数值计算方法具有一定的准确性和可靠性。

图 2 Laval喷管沿程压力分布数据对比

Figure 2. Comparison of pressure distribution data in Laval nozzle

下载: 全尺寸图片幻灯片

4. 不同压力条件下天然气超声速液化特性

保持入口温度及组成(甲烷体积分数90%、乙烷体积分数10%)不变，研究不同的入口压力对Laval喷管内部甲烷-乙烷双组分气体凝结过程中压力、温度、成核率、液滴生长率、液滴半径、液相质量分数的影响。在数值计算中设定的入口温度为270 K，设定入口压力分别为5.5、6.0和6.5 MPa。Laval喷管内双组分气体凝结参数的变化趋势及对比如图 3~图 8所示。

图 3 Laval喷管内气体压力分布

Figure 3. Gas pressure distribution in Laval nozzle

下载: 全尺寸图片幻灯片

图 4 Laval喷管内气体温度分布

Figure 4. Gas temperature distribution in Laval nozzle

下载: 全尺寸图片幻灯片

图 5 Laval喷管内成核率分布

Figure 5. Nucleation rate distribution in Laval nozzle

下载: 全尺寸图片幻灯片

图 6 Laval喷管内液滴半径分布

Figure 6. Droplet radius distribution in Laval nozzle

下载: 全尺寸图片幻灯片

图 7 Laval喷管内液滴生长率分布

Figure 7. Droplet growth rate distribution in Laval nozzle

下载: 全尺寸图片幻灯片

图 8 Laval喷管内液相质量分数分布

Figure 8. Liquid mass fraction distribution in Laval nozzle

下载: 全尺寸图片幻灯片

从压力与温度分布可以看出，气体进入Laval喷管后压力、温度不断降低，当达到一定过冷度时，气体发生凝结并释放潜热，但凝结突跃现象对压力造成的影响并不显著，压力在Laval喷管渐扩段减小到了一个比较稳定状态，温度在减小到最小值后又略微上升，这主要是由于液滴凝结释放潜热引起的。随着入口压力的增大，出口压力略微升高，温升位置有所提前，出口温度也越高，这是因为，随着入口压力的增大，液滴成核与凝结量也随之增多，释放的潜热也就越多。

从成核率分布可以看出，保持其他条件一致，当压力发生变化时，成核速率的变化趋势几乎相同，在刚进入Laval喷管的一段距离内为零，在某一位置处开始，成核率从零开始突跃一直增大到峰值后迅速减小至零。随着入口压力从5.5 MPa增大到6.5 MPa，成核的发生位置(Wilson点)不断向前移动，逐渐向Laval喷管喉部靠拢，且成核率的最值逐渐增大。当压力为5.5 MPa时，成核发生位置为x=0.120 6 m，且在x=0.147 5 m处达到极限成核，为4.044×10²⁰ m^-3·s^-1；当压力为6 MPa时，成核发生位置较5.5 MPa时向前移动，为x=0.119 9 m，极限成核位置也随压力的增大而前移，在x=0.139 3 m处达到8.062×10²⁰ m^-3·s^-1；当压力继续增大为6.5 MPa时，成核发生位置较6 MPa时更加靠近喉部，为x=0.118 2 m，极限成核位置在x=0.132 8 m处，为9.015×10²⁰ m^-3·s^-1。

从液滴半径分布可以看出，随着入口压力的增大，Laval喷管内平均液滴半径越大，出口液滴半径也随之增大，当压力为5.5、6.0、6.5 MPa时，对应最大液滴半径尺寸分别为415.86、447.88和477.44 nm。由此可知，压力的升高有利于液滴的生长。

从液滴生长率分布可以看出，液滴生长率在气体刚进入Laval喷管时一直为零，当液滴开始发生成核凝结时液滴生长率开始突增，变化到最大值后又迅速减小，最终减小为零。综合图 6和图 7还可以看出，随着入口压力的升高，在成核开始时液滴生长率较大，液滴半径增长速度较快，但一段距离后液滴生长率下降更大，液滴半径增长速度也明显放缓。

从液相质量分数分布可以看出，伴随着混合气体的凝结成核，液相质量分数也不断增大，且随着入口压力的升高，Laval喷管出口处的湿度值随之增大，当压力为5.5 MPa时，湿度的最大值为3.989 2%，当压力增大到6.5 MPa时湿度最终增大到7.382 0%。

5. 结论

(1) 建立了三维双组分天然气混合物超声速凝结流动数学模型，对Laval喷管内双组分混合物凝结流动进行了数值模拟，得出沿Laval喷管轴向的参数分布，通过开展双可凝组分气体凝结相变实验，对比发现数值模拟与实验结果基本一致，说明了所建立数学模型及计算方法的正确性。

(2) 利用数值模型研究了入口参数对天然气混合物超声速液化特性的影响，结果表明，保持Laval喷管入口温度及组成不变，增大入口压力，混合气体成核位置前移，成核率、平均液滴半径、液相质量分数均随之增大，即入口压力越大，混合气体在Laval喷管内越易发生凝结。在实际生产中，可以通过调节入口压力来促进天然气的凝结，提高Laval喷管的液化效率。

图 1 （a）DAC中样品、电极俯视图；（b）[100]方向橄榄石样品显微图像；（c）DAC横截面图；（d）四电极法样品显微图像

Figure 1. (a) Top view of the sample and electrodes in a DAC; (b) Microscopic image of the olivine sample along [100] direction;(c) Cross section view of DAC; (d) Microscopic image of the sample and electrodes using the four-probe method

下载: 全尺寸图片幻灯片

图 2 PMX-200硅油v_{2 906}拉曼频移与压力关系

Figure 2. The relation between pressure and Raman shift of v_{2 906}

下载: 全尺寸图片幻灯片

图 3 橄榄石的代表性阻抗谱及其等效电路

Figure 3. Representative impedance spectra and equivalent circuit

下载: 全尺寸图片幻灯片

图 4 300 K下橄榄石单晶的电导率随压力的变化

Figure 4. The conductivity of the single-crystal olivine as a function of pressure at 300 K

下载: 全尺寸图片幻灯片

图 5 橄榄石的平均电导率随压力的变化

Figure 5. Electrical conductivity of olivine aggregates as a function of pressure

下载: 全尺寸图片幻灯片

图 6 橄榄石小极化子导电机制的活化体积随铁含量的变化

Figure 6. Total iron content versus activation volume for the small polaron conduction mechanism in olivine

下载: 全尺寸图片幻灯片

图 7 干的橄榄石不同轴向电导率随深度变化曲线

Figure 7. Conductivity profile of dry olivine in different orientations

下载: 全尺寸图片幻灯片

表 1 样品的EPMA分析结果（质量分数）

Table 1. EPMA analysis results of sample (Mass fraction) %

FeO	MgO	CaO	MnO	SiO₂	Cr₂O₃	NiO	Total
8.78	50.10	0.07	0.11	40.49	0.03	0.35	99.94

下载: 导出CSV

表 2 样品电阻拟合数据

Table 2. Sample resistance fitting data

[100] direction		[010] direction		[001] direction
p/GPa	R/GΩ	p/GPa	R/GΩ	p/GPa	R/GΩ
0	80.8	0	91.7	0	88.1
2.2	67.2	1.8	79.0	1.6	84.1
3.7	61.8	5.3	71.6	3.7	76.9
6.2	57.0	6.6	61.7	6.3	69.7
8.1	52.0	8.0	58.1	8.5	62.1
10.1	46.2	10.1	57.0	10.6	54.3
12.5	41.1	12.7	51.6	12.5	53.2
14.8	34.7	14.7	50.1	13.7	52.1
18.0	28.0	15.8	44.0	14.9	50.3
		19.1	36.0	16.8	49.0
				17.8	47.1

下载: 导出CSV

参考文献(39)

[1]	LARSEN J C. Low frequency (0.1-6.0 CPD) electromagnetic study of deep mantle electrical conductivity beneath the Hawaiian islands [J]. Geophysical Journal International, 1975, 43(1): 17–46. doi: 10.1111/j.1365-246X.1975.tb00626.x
[2]	FILLOUX J H. Ocean-floor magnetotelluric sounding over North Central Pacific [J]. Nature, 1977, 269(5626): 297–301. doi: 10.1038/269297a0
[3]	OLDENBURG D W. Conductivity structure of oceanic upper mantle beneath the Pacific plate [J]. Geophysical Journal International, 1981, 65(2): 359–394. doi: 10.1111/j.1365-246X.1981.tb02717.x
[4]	SHANKLAND T J, O’CONNELL R J, WAFF H S. Geophysical constraints on partial melt in the upper mantle [J]. Reviews of Geophysics, 1981, 19(3): 394–406. doi: 10.1029/RG019i003p00394
[5]	EVANS R L, HIRTH G, BABA K, et al. Geophysical evidence from the MELT area for compositional controls on oceanic plates [J]. Nature, 2005, 437(7056): 249–252. doi: 10.1038/nature04014
[6]	NAIF S, KEY K, CONSTABLE S, et al. Melt-rich channel observed at the lithosphere-asthenosphere boundary [J]. Nature, 2013, 495(7441): 356–359. doi: 10.1038/nature11939
[7]	GAILLARD F, MALKI M, IACONO-MARZIANO G, et al. Carbonatite melts and electrical conductivity in the asthenosphere [J]. Science, 2008, 322(5906): 1363–1365. doi: 10.1126/science.1164446
[8]	KARATO S I. The role of hydrogen in the electrical conductivity of the upper mantle [J]. Nature, 1990, 347(6290): 272–273. doi: 10.1038/347272a0
[9]	WANG D, MOOKHERJEE M, XU Y, et al. The effect of water on the electrical conductivity of olivine [J]. Nature, 2006, 443(7114): 977–980. doi: 10.1038/nature05256
[10]	YOSHINO T, MATSUZAKI T, YAMASHITA S, et al. Hydrous olivine unable to account for conductivity anomaly at the top of the asthenosphere [J]. Nature, 2006, 443(7114): 973–976. doi: 10.1038/nature05223
[11]	YOSHINO T, MATSUZAKI T, SHATSKIY A, et al. The effect of water on the electrical conductivity of olivine aggregates and its implications for the electrical structure of the upper mantle [J]. Earth and Planetary Science Letters, 2009, 288(1/2): 291–300.
[12]	POE B T, ROMANO C, NESTOLA F, et al. Electrical conductivity anisotropy of dry and hydrous olivine at 8 GPa [J]. Physics of the Earth and Planetary Interiors, 2010, 181(3/4): 103–111.
[13]	DUBA A G, SHANKLAND T J. Free carbon & electrical conductivity in the Earth’s mantle [J]. Geophysical Research Letters, 1982, 9(11): 1271–1274. doi: 10.1029/GL009i011p01271
[14]	LASTOVICKOVÁ M. A review of laboratory measurements of the electrical conductivity of rocks and minerals [J]. Physics of the Earth and Planetary Interiors, 1991, 66(1/2): 1–11.
[15]	DAI L, KARATO S. High and highly anisotropic electrical conductivity of the asthenosphere due to hydrogen diffusion in olivine [J]. Earth and Planetary Science Letters, 2014, 408: 79–86. doi: 10.1016/j.jpgl.2014.10.003
[16]	DAI L, KARATO S. The effect of pressure on the electrical conductivity of olivine under the hydrogen-rich conditions [J]. Physics of the Earth and Planetary Interiors, 2014, 232: 51–56. doi: 10.1016/j.pepi.2014.03.010
[17]	YANG X. Orientation-related electrical conductivity of hydrous olivine, clinopyroxene and plagioclase and implications for the structure of the lower continental crust and uppermost mantle [J]. Earth and Planetary Science Letters, 2012, 317: 241–250.
[18]	XU Y, SHANKLAND T J, DUBA A G. Pressure effect on electrical conductivity of mantle olivine [J]. Physics of the Earth and Planetary Interiors, 2000, 118(1/2): 149–161.
[19]	YOSHINO T, SHIMOJUKU A, SHAN S, et al. Effect of temperature, pressure and iron content on the electrical conductivity of olivine and its high-pressure polymorphs [J]. Journal of Geophysical Research: Solid Earth, 2012, 117(B8): 205–220.
[20]	YOSHINO T, ZHANG B, RHYMER B, et al. Pressure dependence of electrical conductivity in forsterite [J]. Journal of Geophysical Research: Solid Earth, 2017, 122(1): 158–171. doi: 10.1002/2016JB013555
[21]	BORUP K A, FISCHER K F, BROWN D R, et al. Measuring anisotropic resistivity of single crystals using the van der Pauw technique [J]. Physical Review B, 2015, 92(4): 045210. doi: 10.1103/PhysRevB.92.045210
[22]	SHEN Y, KUMAR R S, PRAVICA M, et al. Characteristics of silicone fluid as a pressure transmitting medium in diamond anvil cells [J]. Review of Scientific Instruments, 2004, 75(11): 4450–4454. doi: 10.1063/1.1786355
[23]	MAO H K, XU J A, BELL P M. Calibration of the ruby pressure gauge to 800 kbar under quasi-hydrostatic conditions [J]. Journal of Geophysical Research: Solid Earth, 1986, 91(B5): 4673–4676. doi: 10.1029/JB091iB05p04673
[24]	刘锦, 孙樯. 硅油作为压力计的拉曼光谱研究 [J]. 光谱学与光谱分析, 2010, 30(9): 2390–2392. doi: 10.3964/j.issn.1000-0593(2010)09-2390-03 LIU J, SUN Q. Raman spectroscopic study on silicone fluid as pressure gauge [J]. Spectroscopy and Spectral Analysis, 2010, 30(9): 2390–2392. doi: 10.3964/j.issn.1000-0593(2010)09-2390-03
[25]	王晓霞, 李志慧, 陈晨, 等. 硅油的高压拉曼散射 [J]. 高等学校化学学报, 2014, 35(11): 2384–2389. WANG X X, LI Z H, LI C, et al. High pressure Raman spectra of silicone oil [J]. Chemical Journal of Chinese Universities, 2014, 35(11): 2384–2389.
[26]	ROBERTS J J, TYBURCZY J A. Frequency dependent electrical properties of polycrystalline olivine compacts [J]. Journal of Geophysical Research: Solid Earth, 1991, 96(B10): 16205–16222. doi: 10.1029/91JB01574
[27]	SINCLAIR D C, WEST A R. Impedance and modulus spectroscopy of semiconducting BaTiO₃ showing positive temperature coefficient of resistance [J]. Journal of Applied Physics, 1989, 66(8): 3850–3856. doi: 10.1063/1.344049
[28]	JOHNSON D. ZView: a software program for IES analysis. Version 2.8 [CP/OL]. Southern Pines, NC: Scribner Associates [2019-03-05]. http://www.scribner.com.
[29]	ZHA C, DUFFY T S, DOWNS R T, et al. Brillouin scattering and X-ray diffraction of San Carlos olivine: direct pressure determination to 32 GPa [J]. Earth and Planetary Science Letters, 1998, 159(1/2): 25–33.
[30]	DU FRANE W L, TYBURCZY J A. Deuterium-hydrogen exchange in olivine: implications for point defects and electrical conductivity [J]. Geochemistry, Geophysics, Geosystems, 2012, 13(3): Q03004.
[31]	GODDAT A, PEYRONNEAU J, POIRIER J P. Dependence on pressure of conduction by hopping of small polarons in minerals of the Earth’s lower mantle [J]. Physics and Chemistry of Minerals, 1999, 27(2): 81–87. doi: 10.1007/s002690050243
[32]	KATSURA T, SATO K, ITO E. Electrical conductivity of silicate perovskite at lower-mantle conditions [J]. Nature, 1998, 395(6701): 493–495. doi: 10.1038/26736
[33]	LIN J F, WEIR S T, JACKSON D D, et al. Electrical conductivity of the lower-mantle ferropericlase across the electronic spin transition [J]. Geophysical Research Letters, 2007, 34(16): L16305.
[34]	OHTA K, HIROSE K, ONODA S, et al. The effect of iron spin transition on electrical conductivity of (Mg,Fe)O magnesiowüstite [J]. Proceedings of the Japan Academy Series B, 2007, 83(3): 97–100. doi: 10.2183/pjab.83.97
[35]	YOSHINO T, ITO E, KATSURA T, et al. Effect of iron content on the spin transition pressure of ferropericlase [J]. Journal of Geophysical Research: Solid Earth, 2011, 116(B4): B04202.
[36]	DOBSON D P, RICHMOND N C, BRODHOLT J P. A high-temperature electrical conduction mechanism in the lower mantle phase (Mg, Fe)_{1- x}O [J]. Science, 1997, 275(5307): 1779–1781. doi: 10.1126/science.275.5307.1779
[37]	NEAL S L, MACKIE R L, LARSEN J C, et al. Variations in the electrical conductivity of the upper mantle beneath North America and the Pacific Ocean [J]. Journal of Geophysical Research: Solid Earth, 2000, 105(B4): 8229–8242. doi: 10.1029/1999JB900447
[38]	TARITS P, HAUTOT S, PERRIER F. Water in the mantle: results from electrical conductivity beneath the French Alps [J]. Geophysical Research Letters, 2004, 31(6): 265–282.
[39]	JUNG H, KARATO S. Water-induced fabric transitions in olivine [J]. Science, 2001, 293: 1460–1463. doi: 10.1126/science.1062235