高压合成双钙钛矿Y<sub>2</sub>NiIrO<sub>6</sub>中场冷诱导的巨磁电阻效应

彭毅; 邓正; 李文敏; 史鲁川; 赵建发; 张俊; 望贤成; 靳常青

doi:10.11858/gywlxb.20230781

高压合成双钙钛矿Y₂NiIrO₆中场冷诱导的巨磁电阻效应

doi: 10.11858/gywlxb.20230781

彭毅^{1, 2,},
邓正^{1, 2,*, ,},
李文敏^{1, 3},
史鲁川^{1, 2},
赵建发¹,
张俊¹,
望贤成¹,
靳常青^{1, 2,*, ,}

1.
中国科学院物理研究所, 北京　100190
2.
中国科学院大学物理学院, 北京　101408
3.
河南省科学院量子材料与物理研究所, 河南郑州　450046

基金项目: 北京市自然科学基金（2212049）；国家自然科学基金（11974407）

详细信息

作者简介:
彭　毅（1993－），男，博士研究生，主要从事磁性功能材料研究. E-mail：ypeng@iphy.ac.cn

通讯作者:
邓　正（1985－），男，博士，副研究员，主要从事高压功能材料研究. E-mail：dengzheng@iphy.ac.cn

靳常青（1965－），男，博士，研究员，主要从事高压功能材料研究. E-mail：Jin@iphy.ac.cn

中图分类号: O521.2
计量
- 文章访问数: 216
- HTML全文浏览量: 72
- PDF下载量: 35
出版历程
- 收稿日期: 2023-11-08
- 修回日期: 2023-12-14
- 录用日期: 2024-01-17
- 网络出版日期: 2024-01-29
- 刊出日期: 2024-02-05

Cooling Fields Induced Giant Magnetoresistance in High-Pressure Synthesized Double Perovskite Y₂NiIrO₆

PENG Yi^{1, 2
,},
DENG Zheng^{1, 2,*
, ,},
LI Wenmin^{1, 3},
SHI Luchuan^{1, 2},
ZHAO Jianfa¹,
ZHANG Jun¹,
WANG Xiancheng¹,
JIN Changqing^{1, 2,*
, ,}

1.
Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
2.
School of Physics, University of Chinese Academy of Sciences, Beijing 101408, China
3.
Institute of Quantum Materials and Physics, Henan Academy of Sciences, Zhengzhou 450046, Henan, China

摘要

摘要: 双钙钛矿材料Y₂NiIrO₆的亚铁磁转变温度为192 K，因其奇异的交换偏置效应而受到广泛关注。系统研究了Y₂NiIrO₆的低温晶体结构、电导行为及磁电阻性能，发现该材料在130 K时保持了290 K时的晶体结构，并在130～300 K的温区内表现出半导体电导行为。在居里温度以上的顺磁状态，其导电行为可以用Efros-Shklovskii变程跃迁模型拟合，在居里温度以下，亚铁磁有序使电阻行为偏离该模型。更为有趣的是，亚铁磁序诱导了材料的负磁电阻效应，并且7.0 T的场冷诱导了–10%的巨磁电阻效应。这一新机制为探索新型巨磁电阻材料提供了全新的研究思路。
- 高压合成 /
- 双钙钛矿 /
- 交换偏置 /
- 巨磁电阻效应
Abstract: Double perovskite Y₂NiIrO₆ is a ferrimagnetic material with Curie temperature of 192 K. It has drawn wide attention owing to its remarkable exchange-bias effect. Here, we studied the low-temperature crystal structure, electron transport properties and magnetoresistance of Y₂NiIrO₆. The crystal structure under 130 K is almost identical with that of room-temperature. The material shows semiconducting behavior in the temperature range of 130 to 300 K. Above Curie temperature it can be well describe as Efros-Shklovskii variable-range hopping model. Below Curie temperature, a departure occurs due to the forming of long-range ferrimagnetic ordering. It is interesting to find that the magnetic ordering results into negative magneto-resistance. Moreover, giant magnetoresistance up to –10% is induced by cooling field of 7.0 T. This mechanism of this remarkable effect provides a new boulevard to discover new type of giant magneto-resistance materials.
- high-pressure synthesis /
- double perovskite /
- exchange-bias effect /
- giant magnetoresistance

HTML全文

岩体可爆性是指岩石或岩体受到炸药爆炸作用时对炸药释放的能量进行吸收或阻挡的能力。岩体可爆性可以用来衡量岩体对爆破作业的适应性和阻碍自身破碎的难易程度。基于岩体可爆性评级结果，可以进行爆破参数设计和优化，以提高爆破效率、降低爆破风险，达到预期的工程目标。这包括确定合适的爆破药量、装药方式、装药布置以及起爆方式等。总之，岩体可爆性评级和基于评级结果的爆破参数设计是提高爆破效率和安全性的基础。

国内外学者对爆破作用下岩体力学性质和岩体可爆性进行了大量研究^[1–4]。张紫晗等^[5]将岩体脆性指数引入可爆性分级指标中，使用熵权理论对岩体可爆性进行定性和定量评级，该方法可以提供系统的评级结果，帮助爆破参数设计与优化；史涵虚等^[6]基于加权聚类分析原理，对魏家峁煤矿区的岩体可爆性进行了研究，主要探讨了不同岩体的可爆性特征，并通过聚类分析的方法进行评级；王文军等^[7]运用爆破指数法对镜铁山矿业的黑沟山矿区进行了可爆性分级，爆破指数法是一种常用的爆破参数设计方法，通过评估岩体的可爆性特征来确定适当的爆破参数；Wu等^[8]提出了一种基于概念格和粗糙集理论的可爆性分级方法，该方法通过引入概念格和粗糙集理论，提高了分级效率和准确性；Zhou等^[9]建立了一种改进的RES-多维云岩团爆破分类模型，该模型结合了不同的评估指标，现场实际应用效果良好；Alipour等^[10]采用Mamdani模糊算法，用模糊集来表示爆破稳定性指数，并分析了岩体对爆破碎裂的抵抗能力。此外，还有其他研究方法应用于可爆性分级，如神经网络法^[11]、集对分析理论^[12]、CRITIC法与Vague集理论^[13]、综合赋权云模型^[14]等。这些方法通过不同的理论和模型，对岩体可爆性进行多方位、多角度的评估与分级，取得了良好的效果。

以上学者的研究为岩体可爆性分级提供了多种方法和工具，有助于爆破参数的设计与优化，从而提高爆破效率和安全性。然而，在可爆性分级研究中，大多数学者只是对某种岩体进行可爆性评级，评级结果在工程实际中的应用却鲜见报道。基于此，综合考虑指标获取的难易程度，本研究利用相关性分析提高评级指标的独立性，采用正交试验设计确定评级指标的权重，运用模糊决策理论进行岩体可爆性评级，并基于评级结果预测岩体爆破所需炸药单耗的范围。

1. 岩体可爆性分级方法及指标选取

1.1 模糊决策原理

模糊决策方法可以处理不确定和模糊性的问题，适用于岩体可爆性评级。在模糊决策中，可以建立模糊集和隶属函数来表示评级指标与可爆性等级之间的关系。通过模糊推理和模糊逻辑运算，可以实现对岩体可爆性的评级。具体分为5个步骤。

1.1.1 确定研究对象

对于模糊决策而言，首先确定其研究对象：(x₁, x₂, ··· , x_n)为待分类对象的全体。每个样本x_n=(x_1n, x_2n, ··· , x_mn)^T可用m个指标的特征值表示，则样本集可用m×n阶指标特征值矩阵X表示

$\boldsymbol {X} = \left[ {\begin{array}{*{20}{c}} {{x_{11}}}&{{x_{12}}}&{...}&{{x_{1n}}} \\ {{x_{21}}}&{{x_{22}}}&{...}&{{x_{2n}}} \\ {...}&{...}&{...}&{...} \\ {{x_{m1}}}&{{x_{m2}}}&{...}&{{x_{mn}}} \end{array}} \right] = ({x_{ij}})$

(1)

式中：x_ij为第j个样本对应的第i个指标的特征值。

1.1.2 确定评级指标与评级标准

假设有m个评级指标，共分为c个等级，那么其评级标准表对应的特征值矩阵Y为

$\boldsymbol Y = \left[ {\begin{array}{*{20}{c}} {{y_{11}}}&{{y_{12}}}&{...}&{{y_{1c}}} \\ {{y_{21}}}&{{y_{22}}}&{...}&{{y_{2c}}} \\ {...}&{...}&{...}&{...} \\ {{y_{m1}}}&{{y_{m2}}}&{...}&{{y_{mc}}} \end{array}} \right] = ({y_{ih}})$

(2)

式中：y_ih为第i个指标对应的第h个等级所对应的特征值。

1.1.3 评级标准的规格化处理

由于评级指标存在量纲上的差异，为消除其对结果的影响，需要对指标进行规格化处理，使其数据大小处于0～1之间，且量纲均为1。规格化处理可采用

${s_{ih}} = \frac{{{y_{ih}} - {y_{i1}}}}{{{y_{ic}} - {y_{i1}}}}$

(3)

式中：s_ih为第i个指标对应的第h个等级进行规格化处理后的值，且0≤s_ih≤1；y_ic、y_i1分别为第i项指标对应的第c个等级和第1个等级的特征值。

利用式(3)对评级标准进行规格化处理，得到

$\boldsymbol S = \left[ {\begin{array}{*{20}{c}} {{s_{11}}}&{{s_{12}}}&{...}&{{s_{1c}}} \\ {{s_{21}}}&{{s_{22}}}&{...}&{{s_{2c}}} \\ {...}&{...}&{...}&{...} \\ {{s_{m1}}}&{{s_{m2}}}&{...}&{{s_{mc}}} \end{array}} \right] = ({s_{ih}})$

(4)

1.1.4 确定样本的相对隶属度矩阵

设模糊概念β分为c级，在进行模糊决策的时候应先消除指标量纲的影响，基本思路如下。对模糊概念β而言，特征值越大，级别越高，即极大型指标，可以认为：当x_ij大于或等于c级特征值y_ic时，x_ij对于β的相对隶属度为1；反之，当x_ij小于或等于1级特征值y_i1时，其相对隶属度为零；若x_ij介于y_ic与y_i1之间时，则用线性插值法确定其隶属度，隶属度介于0～1之间。相反，若对模糊概念β而言，特征值越小，级别越高，即极小型指标，则认为：当x_ij小于或等于c级特征值y_ic时，x_ij对于β的相对隶属度为1；当x_ij大于或等于1级特征值y_i1时，其相对隶属度为零；x_ij介于两者之间时，则用线性插值确定。

综上，样本对模糊概念β的相对隶属度r_ij的公式为

${r}_{ij}=\left\{\begin{array}{lc}1& \quad{x}_{ij}\geqslant {y}_{ic}\\ \dfrac{{x}_{ij}-{y}_{i1}}{{y}_{ic}-{y}_{i1}}& \quad{y}_{i1} < {x}_{ij} < {y}_{ic}\\ 0&\quad {x}_{ij}\geqslant {y}_{i1}\end{array}\right.\;\;\;\;{\text{或}}\;\;\;\; {r}_{ij}=\left\{\begin{array}{lc}1&\quad {x}_{ij}\leqslant {y}_{ic}\\ \dfrac{{x}_{ij}-{y}_{i1}}{{y}_{ic}-{y}_{i1}}&\quad{y}_{ic} < {x}_{ij} <{y}_{i1} \\ 0& \quad{x}_{ij}\leqslant {y}_{i1}\end{array}\right.$

(5)

式中：r_ij为样本第i个指标的第j个特征值对β的相对隶属度；左边公式对应极大型指标，右侧公式对应极小型指标。

利用式(5)对样本进行处理，得到其相对隶属度矩阵R

$\boldsymbol R = \left[ {\begin{array}{*{20}{c}} {{r_{11}}}&{{r_{12}}}&{...}&{{r_{1c}}} \\ {{r_{21}}}&{{r_{22}}}&{...}&{{r_{2c}}} \\ {...}&{...}&{...}&{...} \\ {{r_{m1}}}&{{r_{m2}}}&{...}&{{r_{mc}}} \end{array}} \right] = ({r_{ij}})$

(6)

1.1.5 确定模糊决策矩阵

样本与评级标准之间的贴近程度用加权距离来衡量，在已知模糊标准矩阵以及各指标的权重矩阵后，模糊决策矩阵中的元素可用下式得出

${\mu _{hj}} = \dfrac{1}{{\displaystyle\sum\limits_{k = 1}^c {{{\left[ {\dfrac{{{{\displaystyle\sum\limits_{i = 1}^m {\left[ {{\omega _i}({r_{ij}} - {s_{ih}})} \right]^p} }}}}{{{{\displaystyle\sum\limits_{i = 1}^m {\left[ {{\omega _i}({r_{ij}} - {s_{ik}})} \right]^p} }}}}} \right]}^{{2}/{p}}}} }}$

(7)

式中：μ_hj为样本j从属于类别h的隶属度；ω_i为第i项指标的权重，$ \sum\limits_{i = 1}^m {{\omega _i}} = 1 $；p为系数。

1.2 岩体可爆性分级指标的选取

岩体可爆性的影响因素复杂多样，针对不同的研究内容，选择的评价指标也有所不同。岩石的单轴抗压强度、单轴抗拉强度、动载冲击强度、抗剪强度、密度、脆性指数、岩体完整性系数等是目前众多学者采用相对较多的评级指标。本研究拟通过相关性分析找出相关性较低的评级指标用于可爆性评级，以避免冗余信息，提高评级指标的独立性。

1.2.1 评级指标的相关性分析

根据岩体爆破理论，岩体破坏是爆破冲击波与爆生气体共同作用的结果。当爆破冲击波和爆生气体对岩体的作用力超过岩体本身的抗剪强度时，岩体就会发生破坏。因此，岩石的抗剪强度可以作为评估岩体可爆性的指标。

脆性^[15]既是一种变形特性，也是一种材料特性。脆性可以用脆性指数度量，脆性指数等于抗压强度与抗拉强度的比值。爆破是一种对岩体施加动态冲击的过程，脆性程度代表岩体对外部冲击的抵抗能力，所以将脆性指数作为评估岩体可爆性的指标是合理的。当岩石的脆性较高时，岩石受到爆破冲击时更容易发生破碎。

璩世杰等^[16]通过对14种岩石的多种物理力学参数进行相关性分析，得到了岩石的静载抗拉强度、密度、完整性系数可以较好地评价岩体的可爆性的结论。在选取岩石单轴抗拉强度、密度、完整性系数作为可爆性评级指标的基础上，为确定能否将岩石的抗剪强度和脆性指数也引入评级指标，选取文献[17–18]中岩石的抗拉强度和抗剪强度以及文献[5, 15, 19–20]中岩石的物理力学参数进行相关性分析。岩石的抗拉强度（σ_t）和抗剪强度（τ）见表1。岩石的物理力学参数见表2，其中：ρ为岩石密度，η为岩体的完整性系数，B为岩石脆性指数。共分析测试39块岩石样品的物理力学参数，表2仅列出12块岩石样品的数据。

表 1 岩石的抗拉强度和抗剪强度

Table 1. Tensile and shear strength of rocks

Rock category	σ_t/MPa	τ/MPa	Rock category	σ_t/MPa	τ/MPa
Magnetite ore	12	43	Migmatite	8.5	38
Chlorite schist	2.6	12	Quartz schist	7	16
Magnetite rich ore	11.5	40	Mixed metamorphic rock	4.81	21.66
Amphibolite	7.07	22	Mixed granite	12.93	47.85
k shell	8	32	Fault breccia	8.53	33.02

下载: 导出CSV

| 显示表格

表 2 岩石的物理力学参数

Table 2. Physical and mechanical parameters of rock

Rock No.	ρ/(g·cm⁻³)	σ_t/MPa	η	B
1	2.618	9.30	0.767	12.74
2	2.674	14.25	0.760	5.46
3	3.810	5.85	0.510	16.20
4	2.710	7.05	0.410	22.70
5	2.646	14.62	0.780	4.40
6	2.588	7.61	0.830	5.27
7	2.919	15.33	0.793	3.94
8	2.631	9.71	0.763	16.55
9	3.418	16.30	0.230	11.87
10	3.029	21.04	0.197	7.27
11	2.725	14.98	0.480	8.01
12	3.479	18.26	0.779	10.55

下载: 导出CSV

| 显示表格

分别对表1与表2中的数据进行相关性分析，结果见表3，其中：β₀、β₁为回归方程系数，R为相关度。从结果看：岩石的抗拉强度与抗剪强度的相关度为0.8138，相关度较高，不可将两者同时用作可爆性评级指标；岩石脆性指数与岩石密度、单轴抗拉强度、完整性系数的相关度分别为0.0084、0.2386、0.0273，基本线性无关，即可以在岩石密度、单轴抗拉强度、完整性系数作为评级指标的前提下，将脆性指数作为评级指标的一部分。

表 3 相关性分析结果

Table 3. Correlation analysis results

Results	β₀	β₁	R
τ=0.3658+3.5369σ_t	0.3658	3.5369	0.8138
B=5.2149+2.0844ρ	5.2149	2.0844	0.0084
B=19.9728−0.7139σ_t	19.9728	−0.7139	0.2386
B=14.3477−5.4339η	14.3477	−5.4339	0.0273

下载: 导出CSV

| 显示表格

相对于抗拉强度，抗剪强度（内聚力与内摩擦角）的获取途径复杂，为此，本研究最终确定以岩石密度、单轴抗拉强度、完整性系数、岩石脆性指数作为可爆性的评级指标。

1.3 权重系数的确定

在岩体可爆性评级过程中，各分级指标权重的确定是至关重要的一步，它体现岩体的物理力学参数对最终评级结果的贡献度。本研究选择应用每个因子3个水平的正交试验设计确定4个评级指标的敏感性，进而确定其权重。通过正交设计，可以高效地确定权重，减少主观性和随机性对权重的影响。正交设计中各因素水平见表4，L₉(3⁴)型正交试验设计见表5，正交试验设计结果见表6。利用极差分析方法对结果进行分析，将各列水平数相同的结果相加，记为K，3个水平的计算结果分别为K₁、K₂、K₃。表6中的等级基于文献[12]中提到的方法得出。

表 4 正交设计各因素水平

Table 4. Levels of various factors in orthogonal design

Level	ρ/(g·cm⁻³)	σ_t/MPa	B	η
1	2.618	8.61	10.77	0.360
2	2.975	13.38	7.79	0.513
3	3.270	20.14	5.46	0.783

下载: 导出CSV

| 显示表格

表 5 正交设计试验L₉(3⁴)

Table 5. Orthogonal design test L₉(3⁴)

Test No.	Level
Test No.	ρ	σ_t	B	η
1	1	1	1	1
2	1	2	2	2
3	1	3	3	3
4	2	1	2	3
5	2	2	3	1
6	2	3	1	2
7	3	1	3	2
8	3	2	1	3
9	3	3	2	1

下载: 导出CSV

| 显示表格

表 6 正交试验设计结果

Table 6. Orthogonal design test results

Test No.	ρ/(g·cm⁻³)	σ_t/MPa	B	η	Level^[12]
1	2.618	8.61	10.77	0.360	2
2	2.618	13.38	7.79	0.513	4
3	2.618	20.14	5.46	0.783	5
4	2.975	8.61	7.79	0.783	7
5	2.975	13.38	5.46	0.360	3
6	2.975	20.14	10.77	0.513	4
7	3.270	8.61	5.46	0.513	6
8	3.270	13.38	10.77	0.783	6
9	3.270	20.14	7.79	0.360	5
K₁	11	15	12	10
K₂	14	13	16	14
K₃	17	14	14	18
Range	2	0.667	1.333	2.667
Sensitiveness	2	4	3	1

下载: 导出CSV

| 显示表格

通过表6可以看出，对于评价岩体可爆性的4个指标，按照敏感性由大到小依次为：完整性系数、岩体密度、岩石脆性指数、岩石单轴抗拉强度，因此在确定评级指标的权重时应予以考虑。根据极差分析结果，完整性系数、岩体密度、岩石脆性指数、岩石单轴抗拉强度的极差比值为4∶3∶2∶1，由此确定其权重为(0.4, 0.3, 0.2, 0.1)。文献[6,18]指出，岩石的单轴抗拉强度对可爆性评级的影响较大，最终将权重调整为(0.35, 0.29, 0.21, 0.15)。

2. 岩体可爆性等级评定

选取文献[5]中的12种岩体试样作为样本，样本集指标特征值矩阵X为

${ \boldsymbol {X}} ={\left[ {\begin{array}{*{20}{c}} {2.618}&{2.674}&{2.604}&{3.070}&{2.875}&{2.813}&{2.962}&{3.368}&{3.388}&{3.445}&{3.325}&{3.479} \\ {9.30}&{14.25}&{11.12}&{17.38}&{9.41}&{6.61}&{9.00}&{15.83}&{25.66}&{9.34}&{15.27}&{18.26} \\ {12.74}&{5.46}&{11.64}&{9.02}&{7.79}&{3.13}&{8.69}&{9.11}&{6.41}&{13.90}&{10.77}&{10.55} \\ {0.767}&{0.760}&{0.796}&{0.783}&{0.783}&{0.767}&{0.783}&{0.783}&{0.774}&{0.783}&{0.763}&{0.779} \end{array}} \right]} = {({x_{ij}})_{4 \times 12}}$

(8)

基于以上分析，以岩石单轴抗拉强度、岩石密度、岩石脆性指数、岩体完整性系数为评级指标，以文献[5, 21]为参考，建立可爆性分级标准，如表7所示，共分为7个等级，即最易爆Ⅰ、易爆Ⅱ、较易爆Ⅲ、中等Ⅳ、较难爆Ⅴ、难爆Ⅵ、极难爆Ⅶ。

表 7 可爆性分级标准

Table 7. Explosivity classification standards

Rock mass explosivity level	ρ/(g·cm⁻³)	σ_t/MPa	B	η	Level
Ⅰ	2.50	6.6	17.5	0.0494	Easiest
Ⅱ	2.60	10.0	15.5	0.2555	Easy
Ⅲ	2.75	13.0	13.3	0.3654	Easier
Ⅳ	2.90	17.0	11.8	0.5122	Medium
Ⅴ	3.16	20.0	9.5	0.6021	More difficult
Ⅵ	3.30	23.0	6.0	0.7122	Difficult
Ⅶ	3.45	26.0	3.5	0.8232	Most difficult

下载: 导出CSV

| 显示表格

取表7中的数据建立可爆性评级标准对应的特征值矩阵Y

$\boldsymbol Y = \left[ {\begin{array}{*{20}{c}} {2.50}&{2.60}&{2.75}&{2.90}&{3.16}&{3.30}&{3.45} \\ {6.6}&{10.0}&{13.0}&{17.0}&{20.0}&{23.0}&{26.0} \\ {17.5}&{15.5}&{13.3}&{11.8}&{9.5}&{6.0}&{3.5} \\ {0.049\,4}&{0.255\,5}&{0.365\,4}&{0.512\,2}&{0.602\,1}&{0.712\,2}&{0.823\,2} \end{array}} \right] = {({y_{ij}})_{4 \times 7}}$

(9)

采用式(3)对Y进行规格化处理，得到消除量纲影响后的标准特征值矩阵S

$\boldsymbol S = \left[ {\begin{array}{*{20}{c}} 0&{0.105\,3}&{0.263\,2}&{0.421\,1}&{0.694\,7}&{0.842\,1}&1 \\ 0&{0.175\,3}&{0.329\,9}&{0.536\,1}&{0.690\,7}&{0.845\,4}&1 \\ 0&{0.142\,9}&{0.300\,0}&{0.407\,1}&{0.571\,4}&{0.821\,4}&1 \\ 0&{0.266\,3}&{0.408\,4}&{0.598\,1}&{0.714\,3}&{0.856\,6}&1 \end{array}} \right] = {({s_{ih}})_{4 \times 7}}$

(10)

利用式(5)对X进行规格化处理，确定样本的相对隶属度矩阵R

$\boldsymbol R = {\left[ {\begin{array}{*{20}{c}} {0.124\,2}& {0.183\,2}& {0.109\,5}& {0.600\,0}& {0.394\,7}& {0.329\,5}& {0.486\,3}& {0.913\,7}& {0.934\,7}& {0.994\,7}& {0.868\,4}& 1.000\,0 \\ {0.139\,2}& {0.394\,3}& {0.233\,0}& {0.555\,7}& {0.144\,8}& {0.000\,5}& {0.123\,7}& {0.475\,8}& {0.9825}& {0.141\,2}& {0.446\,9}& {0.601\,0} \\ {0.340\,0}& {0.860\,0}& {0.418\,6}& {0.605\,7}& {0.693\,6}& {1.000\,0}& {0.629\,3}& {0.599\,3}& {0.792\,1}& {0.257\,1}& {0.480\,7}& {0.496\,4} \\ {0.927\,4}& {0.918\,3}& {0.964\,8}& {0.948\,0}& {0.948\,0}& {0.927\,4}& {0.948\,0}& {0.948\,0}& {0.936\,4}& {0.948\,0}& {0.922\,2}& {0.942\,9} \end{array}} \right] } = {({r_{ij}})_{4 \times 12}}$

(11)

通过式(7)进行计算，取p=2，计算得到样本的模糊决策矩阵U

$\boldsymbol U = {\left[ {\begin{array}{*{20}{c}} {0.061\,5}& {0.047\,3}& {0.059\,2}& {0.017\,2}& {0.037\,0}& {0.047\,8}& {0.032\,8}& {0.013\,9}& {0.004\,4}& {0.034\,9}& {0.016\,9}& {0.016\,1} \\ {0.124\,7}& {0.088\,1}& {0.116\,2}& {0.031\,2}& {0.069\,2}& {0.083\,4}& {0.060\,9}& {0.022\,7}& {0.007\,0}& {0.055\,3}& {0.028\,5}& {0.025\,7} \\ {0.194\,2}& {0.145\,2}& {0.180\,2}& {0.055\,7}& {0.119\,7}& {0.132\,0}& {0.106\,3}& {0.037\,0}& {0.011\,0}& {0.084\,2}& {0.047\,6}& {0.040\,7} \\ {0.282\,2}& {0.250\,8}& {0.274\,9}& {0.144\,6}& {0.242\,1}& {0.212\,9}& {0.227\,0}& {0.075\,4}& {0.021\,6}& {0.145\,5}& {0.102\,6}& {0.079\,2} \\ {0.163\,6}& {0.205\,3}& {0.173\,9}& {0.355\,3}& {0.247\,7}& {0.210\,9}& {0.273\,5}& {0.237\,1}& {0.038\,2}& {0.270\,1}& {0.324\,1}& {0.233\,2} \\ {0.105\,5}& {0.161\,6}& {0.118\,1}& {0.284\,1}& {0.180\,8}& {0.188\,6}& {0.193\,6}& {0.417\,2}& {0.519\,2}& {0.243\,6}& {0.333\,0}& {0.379\,6} \\ {0.068\,3}& {0.101\,8}& {0.077\,3}& {0.112\,1}& {0.103\,5}& {0.124\,4}& {0.105\,7}& {0.196\,7}& {0.368\,6}& {0.166\,4}& {0.147\,3}& {0.225\,5} \end{array}} \right]} = {({u_{hj}})_{7 \times 12}}$

(12)

根据最大隶属度原则得到可爆性等级结果，见表8。参考文献[5, 21]，采用熵权理论法和属性识别法也得到可爆性等级，见表8。

表 8 可爆性分级结果

Table 8. Explosion classification results

Rock No.	Proposed method	Attribute recognition method	Entropy weight theory method
1	Ⅳ	Ⅲ	Ⅲ
2	Ⅳ	Ⅲ	Ⅲ
3	Ⅳ	Ⅴ	Ⅴ
4	Ⅴ	Ⅴ	Ⅵ
5	Ⅴ	Ⅴ	Ⅴ
6	Ⅳ	Ⅳ	Ⅳ
7	Ⅴ	Ⅳ	Ⅴ
8	Ⅵ	Ⅵ	Ⅵ
9	Ⅵ	Ⅵ	Ⅵ
10	Ⅴ	Ⅵ	Ⅶ
11	Ⅵ	Ⅵ	Ⅵ
12	Ⅵ	Ⅵ	Ⅵ

下载: 导出CSV

| 显示表格

从结果来看，本研究得到的评级结果与其他两种方法得到的结果具有较好的一致性，说明模糊决策方法具有一定的可行性，同时也证明采用正交试验设计确定各评级指标的权重具有合理性。然而，从实际应用的角度来看，模糊决策方法更适合处理数据，它基于模糊集理论，可以处理模糊的、有限的、不完全的和不精确的数据，使决策过程更加灵活，可以适应各种实际情况。此外，模糊决策方法还可以方便地植入数学软件中进行计算和分析，有许多基于模糊决策方法的软件工具和数据库可供使用，它们提供了专门的算法和函数，能够快速进行模糊推理和评级计算。这些软件工具使实施模糊决策方法变得简单，大大减少了计算工作量和人工操作的复杂性。

岩体的可爆性评级是多个因素共同作用的结果，任何一种因素发生变化都会对评级结果产生影响。在本研究的评级方法中，岩石1和岩石2的评级结果比其他两种方法高一级，而岩石3和岩石10的评级结果则比其他两种方法低，这是因为3种评级方法的权重分配方案不同。对于岩石10，其密度、抗拉强度、完整性系数、脆性指数分别为3.445 g/cm³、9.34 MPa、0.783、13.90，与其他岩石相比，其脆性指数较大，且脆性指数是极小型指标，其值越大代表岩体可爆性等级越低，岩体越容易爆破。相对其余两种方法，本研究赋予脆性指数的权重略大，因此岩石10的评级结果比其他两种方法低。

3. 可爆性等级对应炸药单耗区间的预测

基于可爆性评级结果预测岩体爆破所需炸药单耗的范围是评价岩体可爆性的一种补充应用，可为工程实际应用提供更详细的指导信息。本研究利用多元回归分析方法，对可爆性等级与炸药单耗之间的关系建立数学模型，从而实现对炸药单耗范围预测。

3.1 炸药单耗预测公式推导

对于非煤矿山地下采矿浅孔爆破，爆破每立方米岩体所需的炸药量q通常要根据不同岩石的坚固性系数f进行调整，具体如表9^[22]所示。

表 9 炸药单耗与岩石坚固性系数^[22]

Table 9. Unit consumption and rock soundness coefficient of explosive^[22]

f	q/(kg·m⁻³)	f	q/(kg·m⁻³)
0–4	0.25–0.65	10–15	1.60–2.60
5–8	0.65–1.00	>15	>2.80
8–10	1.00–1.60

下载: 导出CSV

| 显示表格

参考《工程岩体分级标准》（GB/T 50218—2014）^[23]，基于岩体基本质量指标BQ，将岩体分为5级，将其与表9结合得到表10。

表 10 $BQ $与炸药单耗

Table 10. $BQ $ and explosive unit consumption

BQ	q/(kg·m⁻³)	BQ	q/(kg·m⁻³)
≤250	0.25–0.65	451–550	1.60–2.60
251–350	0.65–1.00	>551	>2.80
351–450	1.00–1.60

下载: 导出CSV

| 显示表格

利用MATLAB将表10中的数据进行拟合，得到炸药单耗q与BQ之间的关系，见表11，q_ll为炸药单耗的下限，q_ul为炸药单耗的上限。

表 11 炸药单耗与$BQ $拟合结果

Table 11. Fitting results of explosive unit consumption and $BQ $

Result	β₀	β₁	R
q_ll=–0.6180+0.0055BQ	–0.6180	0.0055	0.9629
q_ul=–0.7425+0.0075BQ	−0.7425	0.0075	0.9425

下载: 导出CSV

| 显示表格

岩体基本质量指标BQ计算公式为

$BQ=100+3R_{\mathrm{c}}+250\eta$

(13)

式中：R_c为岩石的单轴饱和抗压强度。当R_c>90η+30时，将η和R_c=90η+30代入式(13)求BQ；当η>0.04R_c+0.4时，将η=0.04R_c+0.4和R_c代入式(13)求BQ的值。

基于此，得到39种岩石的BQ，如表12所示，与表2一致，这里仅列其中12种岩石的BQ。

表 12 岩石的$BQ $

Table 12. Rock $BQ $

Rock No.	BQ	Rock No.	BQ
1	588	7	479
2	523	8	588
3	455	9	309
4	403	10	597
5	489	11	440
6	428	12	595

下载: 导出CSV

| 显示表格

根据工程实际，对岩体的BQ进行修正，利用表2中线性相关度较弱的数据，将岩体密度、抗拉强度、脆性指数、完整性系数与BQ进行拟合，得到的拟合公式为：BQ=57.4－5.5ρ+7.41σ_t+460.2η+5.22B，相关度为0.913。将该式与表11中炸药单耗拟合公式结合，得到基于评级指标的炸药单耗预测公式。

完整性系数可以由岩体与岩石波速计算得到，即知道岩体与岩石波速、岩体的密度、抗拉强度、抗压强度，就可通过式(14)和式(15)计算出炸药单耗的上限和下限。

$q_{\mathrm{ul}}=-0.311\,9-0.041\,2\rho+0.055\,6\sigma\mathrm{_t}+3.451\,4\eta+0.039\,1B$

(14)

$q\mathrm{_{ll}}=-0.302\,2-0.030\,2\rho+0.040\,7\sigma\mathrm{_t}+2.531\eta+0.028\,7B$

(15)

3.2 可爆性等级对应的炸药单耗

以39种岩石为样本，采用本研究的方法对其进行评级，利用式(14)和式(15)，计算炸药单耗的范围，结果见表13。

表 13 评级与炸药单耗计算结果

Table 13. Rating and explosive unit consumption calculation results

Level	Number	q/(kg·m⁻³)	Level	Number	q/(kg·m⁻³)
Ⅰ	0	<0.77	Ⅴ	12	1.91–2.71
Ⅱ	3	0.77–1.22	Ⅵ	5	2.48–3.49
Ⅲ	8	1.17–1.69	Ⅶ	0	>3.49
Ⅳ	11	1.58–2.26

下载: 导出CSV

| 显示表格

由表13可以看出，39种岩石中可爆性等级为Ⅱ、Ⅲ、Ⅳ、Ⅴ、Ⅵ级的分别有3、8、11、12、5种。将评级结果与式(14)和式(15)计算得到的炸药单耗区间进行对应，即可得到各可爆性等级对应的炸药单耗预测范围。

3.3 炸药单耗预测的有效性验证

为验证炸药单耗预测范围的有效性，在云南某锡矿开展岩石的物理力学实验，并进行现场爆破试验。试验地点位于三坑1480中段，爆破试验的对象岩体为大理岩和氧化矿。

3.3.1 岩石的室内实验

采用微机电液伺服试验机HYE-2000对矿岩试样进行单轴抗压和单轴抗拉实验，采用ZT801岩体参数测定仪测量矿岩试样的波速，通过称量法测量矿岩试样的密度，得到的岩石物理力学数据见表14。

表 14 岩矿试样的物理参数

Table 14. Physical parameters of rock specimen

Rock type	ρ/(g·cm⁻³)	σ_t/MPa	B	η
Marble	2.72	4.66	12.9	0.566
Oxidized ore	2.47	4.86	15.1	0.156

下载: 导出CSV

| 显示表格

根据表14所列数据，采用本研究的方法对岩矿进行可爆性评级：大理岩的可爆性等级为Ⅳ级，氧化矿的可爆性等级为Ⅰ级。由表13即可得到其对应的炸药单耗范围，大理岩的炸药单耗范围为1.58～2.26 kg/m³，氧化矿的炸药单耗小于0.77 kg/m³。

3.3.2 现场爆破试验

分别进行1次大断面大理岩巷道掘进的爆破试验以及1次氧化矿采场爆破试验。试验采用大直径空孔直眼掏槽，装药孔与空孔的距离为200 mm，为空孔直径D的2.0倍，起爆网路采用微差起爆技术。大理岩掏槽孔布置见图1，爆破效果见图2。氧化矿装药作业以及连线作业分别见图3和图4。

图 1 大理岩掏槽孔布置

Figure 1. Layout of marble cutting holes

下载: 全尺寸图片幻灯片

图 2 大理岩爆破效果

Figure 2. Effect of marble blasting

下载: 全尺寸图片幻灯片

图 3 氧化矿装药作业

Figure 3. Operation of oxide ore charging

下载: 全尺寸图片幻灯片

图 4 氧化矿连线作业

Figure 4. Operation of oxidation ore connection

下载: 全尺寸图片幻灯片

爆破后，对出矿、进尺、孔深等相关数据进行统计，发现无特殊大块产生，大理岩爆破的平均炸药单耗为2.04 kg/m³，氧化矿爆破的平均炸药单耗为0.71 kg/m³，处于预测的炸药单耗范围内，这在一定程度上证明了预测范围的合理性，为相似工况的炸药单耗设计提供了一定的参考。

4. 结　论

(1) 通过相关性分析得到了抗拉强度、岩石密度、脆性指数和岩体完整性系数基本不相关的结论，4个量能提供独立的信息，可同时作为可爆性评级指标；

(2) 通过设计一系列正交试验，确定了评级指标对结果的影响，在此基础上确定了各指标权重；

(3) 模糊决策理论是一套科学合理的决策理论，为岩体可爆性评级提供了一种新思路，通过评级结果的对比，证明了该方法的可靠性；

(4) 结合可爆性评级指标和岩体基本质量指标，推导出了炸药单耗的预测公式，进一步给出了各可爆性等级对应的炸药单耗预测区间，现场爆破试验证明了其合理性。

图 1 YNIO在290和130 K的XRD谱及Rietveld精修结果

Figure 1. XRD spectra and the corresponding Rietveld refinements of YNIO at 290 and 130 K

下载: 全尺寸图片幻灯片

图 2 0.5 T外场下磁矩随温度的变化曲线

Figure 2. Temperature dependence of magnetization under external field of 0.5 T

下载: 全尺寸图片幻灯片

图 3 130 K下经历零场冷、0.5 T场冷和7.0 T场冷后的磁滞回线

Figure 3. Field dependence of magnetization after zero field cooling, field cooling of 0.5 and 7.0 T at 130 K

下载: 全尺寸图片幻灯片

图 4 不同磁场下YNIO的电阻率-温度变化曲线

Figure 4. Temperature dependence of resistivity of YNIO under different fields

下载: 全尺寸图片幻灯片

图 5 零场下ρ(T)的变程跃迁拟合

Figure 5. Fitting of Efros-Shklovskii variable-rangehopping at zero field cooling

下载: 全尺寸图片幻灯片

图 6 零场冷和7.0 T场冷下磁电阻随温度的变化

Figure 6. Temperature dependence of magnetoresistance under ZFC and FC at 7.0 T

下载: 全尺寸图片幻灯片

图 7 130 K时零场冷、0.5 T场冷和7.0 T场冷后磁电阻与温度的变化曲线

Figure 7. Field dependence of magnetoresistance after zero field cooling, field cooling of 0.5 and 7.0 T at 130 K

下载: 全尺寸图片幻灯片

表 1 YNIO在290和130 K下的晶格参数

Table 1. Structural parameters of YNIO at 290 and 130 K

T/K	Space group	a/Å		b/Å	c/Å	β/(^o)	B'-B'' antisite/%
290	P2₁/n	5.265 2		5.684 7	7.584 6	90.14	7.5
130	P2₁/n	5.263 1		5.683 0	7.580 9	90.16	7.5

T/K	Position
T/K	Y	Ni	Ir	O₁	O₂		O₃
290	0.021 8, 0.077 5 0.249 0	1/2, 0, 1/2	1/2, 0, 0	0.184 3, –0.191 4, 0.055 3	0.619 5, –0.042 6, 0.254 7		0.323 3, 0.307 3, 0.059 8
130	0.022 1, 0.078 8, 0.249 1	1/2, 0, 1/2	1/2, 0, 0	0.183 9, –0.189 3, 0.051 7	0.618 4, –0.047 1, 0.253 3		0.322 2, 0.308 0, 0.061 2

下载: 导出CSV

参考文献(25)

[1]	WU S M, CYBART S A, YU P, et al. Reversible electric control of exchange bias in a multiferroic field-effect device [J]. Nature Materials, 2010, 9(9): 756–761.
[2]	PARKIN S S P, ROCHE K P, SAMANT M G, et al. Exchange-biased magnetic tunnel junctions and application to nonvolatile magnetic random access memory (invited) [J]. Journal of Applied Physics, 1999, 85(8): 5828–5833.
[3]	LEIGHTON C, SONG M, NOGUÉS J, et al. Using magnetoresistance to probe reversal asymmetry in exchange biased bilayers [J]. Journal of Applied Physics, 2000, 88(1): 344–347.
[4]	BOBO J F, GABILLET L, BIBES M. Recent advances in nanomagnetism and spin electronics [J]. Journal of Physics: Condensed Matter, 2004, 16(5): S471–S496.
[5]	OHLDAG H, SCHOLL A, NOLTING F, et al. Correlation between exchange bias and pinned interfacial spins [J]. Physical Review Letters, 2003, 91(1): 017203.
[6]	SHIRATSUCHI Y, NOUTOMI H, OIKAWA H, et al. Detection and in situ switching of unreversed interfacial antiferromagnetic spins in a perpendicular-exchange-biased system [J]. Physical Review Letters, 2012, 109(7): 077202.
[7]	LIU Z H, ZHANG Y J, ZHANG H G, et al. Giant exchange bias in Mn₂FeGa with hexagonal structure [J]. Applied Physics Letters, 2016, 109(3): 032408.
[8]	DING L, CHU L H, MANUEL P, et al. Giant spontaneous exchange bias in an antiperovskite structure driven by a canted triangular magnetic structure [J]. Materials Horizons, 2019, 6(2): 318–325.
[9]	MANIV E, MURPHY R A, HALEY S C, et al. Exchange bias due to coupling between coexisting antiferromagnetic and spin-glass orders [J]. Nature Physics, 2021, 17(4): 525–530.
[10]	PATRA M, MAJUMDAR S, GIRI S. Exchange bias effect and intragranular magnetoresistance in Nd_0.84Sr_0.16CoO₃ [J]. Journal of Physics: Condensed Matter, 2009, 21(48): 486003. doi: 10.1088/0953-8984/21/48/486003
[11]	GIRI S, PATRA M, MAJUMDAR S. Exchange bias effect in alloys and compounds [J]. Journal of Physics: Condensed Matter, 2011, 23(7): 073201. doi: 10.1088/0953-8984/23/7/073201
[12]	JANA S, MIDDEY S, RAY S. Spin-valve-type magnetoresistance: a generic feature of ferromagnetic double perovskites [J]. Journal of Physics: Condensed Matter, 2010, 22(34): 346004. doi: 10.1088/0953-8984/22/34/346004
[13]	DENG Z, WANG X, WANG M Q, et al. Giant exchange-bias-like effect at low cooling fields induced by pinned magnetic domains in Y₂NiIrO₆ double perovskite [J]. Advanced Materials, 2023, 35(17): 2209759.
[14]	CHENG J G, ZHOU J S, GOODENOUGH J B, et al. High-pressure synthesis and physical properties of perovskite and post-perovskite Ca_{1− x}Sr_xIrO₃ [J]. Physical Review B, 2011, 83(6): 064401. doi: 10.1103/PhysRevB.83.064401
[15]	PENG S Z, ZHU D Q, LI W X, et al. Exchange bias switching in an antiferromagnet/ferromagnet bilayer driven by spin-orbit torque [J]. Nature Electronics, 2020, 3(12): 757–764.
[16]	WEI Q Q, WANG H L, ZHAO X P, et al. Electron mobility anisotropy in (Al, Ga)Sb/InAs two-dimensional electron gases epitaxied on GaAs (001) substrates [J]. Journal of Semiconductors, 2022, 43(7): 072101. doi: 10.1088/1674-4926/43/7/072101
[17]	DENG Z, RETUERTO M, LIU S Z, et al. Dynamic ferrimagnetic order in a highly distorted double perovskite Y₂CoRuO₆ [J]. Chemistry of Materials, 2018, 30(20): 7047–7054.
[18]	DENG Z, KANG C J, CROFT M, et al. A pressure-induced inverse order-disorder transition in double perovskites [J]. Angewandte Chemie International Edition, 2020, 59(21): 8240–8246. doi: 10.1002/anie.202001922
[19]	TOBY B H. EXPGUI, a graphical user interface for GSAS [J]. Journal of Applied Crystallography, 2001, 34(2): 210–213. doi: 10.1107/S0021889801002242
[20]	BUCHNER M, HENNE B, NEY V, et al. Transition from a hysteresis-like to an exchange-bias-like response of an uncompensated antiferromagnet [J]. Physical Review B, 2019, 99(6): 064409. doi: 10.1103/PhysRevB.99.064409
[21]	NAYAK A K, NICKLAS M, CHADOV S, et al. Design of compensated ferrimagnetic Heusler alloys for giant tunable exchange bias [J]. Nature Materials, 2015, 14(7): 679–684. doi: 10.1038/nmat4248
[22]	NAYAK A K, NICKLAS M, CHADOV S, et al. Large zero-field cooled exchange-bias in bulk Mn₂PtGa [J]. Physical Review Letters, 2013, 110(12): 127204. doi: 10.1103/PhysRevLett.110.127204
[23]	NAKATSUJI S, DOBROSAVLJEVIĆ V, TANASKOVIĆ D, et al. Mechanism of hopping transport in disordered mott insulators [J]. Physical Review Letters, 2004, 93(14): 146401. doi: 10.1103/PhysRevLett.93.146401
[24]	FISHER B, GENOSSAR J, CHASHKA K B, et al. Variable range hopping in A₂MnReO₆ (A=Ca, Sr, Ba) [J]. Journal of Applied Physics, 2008, 104(3): 033716. doi: 10.1063/1.2967820
[25]	LIU X Y, RINEY L, GUERRA J, et al. Colossal negative magnetoresistance from hopping in insulating ferromagnetic semiconductors [J]. Journal of Semiconductors, 2022, 43(11): 112502. doi: 10.1088/1674-4926/43/11/112502