一种基于磁梯度张量反演磁化率张量的方法与流程

[0001]
本发明属于磁化率张量测量技术领域，具体涉及一种基于岩石样本的磁梯度张量测试方法，特别涉及一种基于磁梯度张量反演磁化率张量的方法。


背景技术：

[0002]
岩石内部的矿物质颗粒或者矿物质晶体在空间分布、相对比例配置上所表现出来的某种方向性，也应该在其相关的物理性质上有所表现，岩矿石的磁化率张量就是岩石磁学性质之一。岩矿石的磁化率张量包含着丰富的地质信息，对于解决地质构造、岩石生成等地质问题发挥着重要的作用。此外，只测量岩石的磁化率张量要比常规的岩组构测量速度快很多。所以，在地质工作、地球物理、测井技术中应该重视并进一步研究测定岩石的磁化率张量。
[0003]
目前，最常用的磁化率张量测定方法为磁感应法和静态法。在磁感应法中，需将岩石标本放入通电线圈中，线圈的磁感应通量就会变化。保持其他条件不变的情况下，线圈通量的改变量与磁化率大小成正比，比例系数为定值，该定值可用已知的磁化率标准样品来测定。在静态法中，岩石标本只受到地磁场的磁化，磁钢的影响可以不用考虑，这就使得把岩石标本的磁矩和磁化强度联系起来，磁化强度中既包含感应磁化强度(m
i
)又包含剩余磁化强度(m
r
)，感应磁化强度与磁化率大小成正比，通过测量感应磁化强度与外磁场，便能得到磁化率值。
[0004]
但上述两种测量方式都需通过专门的仪器来测量，磁感应法通常使用卡帕桥磁化率仪，静态法通常使用无定向磁力仪来完成。而通过反演的方法能多种途径来获得磁化率张量数据，并且，能与上述仪器相互配合，进行验证。
[0005]
在传统的磁测方法中，通常使用地磁总场强度数据来反演求磁化率值，但在实际应用当中存在两个缺点：一是地磁总场容易受到地理位置的影响；二是容易受到时变磁场的影响，需单独在配备日变校正仪，并且求得的是磁化率标量，而不是磁化率张量数据，这对于研究古流向造成的磁性矿物的定向排列，以及构造应力作用引起的岩石内磁性矿物的定向重结晶，定向排列及韧性变形不具有任何价值。


技术实现要素：

[0006]
本发明的目的在于提供一种利用磁梯度张量数据来反演求解磁化率张量的方法，以解决现有磁测方法需专门的测量仪器，且易受地理位置及时变磁场影响，无法得到完整的磁化率张量信息的问题。本方法对静态法进行了改进，从静态法的理论出发，用磁梯度张量值替换了静态方程组中的磁矩，在实验中直接测量岩石样本的磁梯度张量值。
[0007]
本发明的目的是通过以下技术方案实现的：
[0008]
一种基于磁梯度张量反演磁化率张量的方法，包括以下步骤：
[0009]
a、将岩石样本置于三轴无磁转台3上，在岩石样本所处坐标系的xy平面沿x轴距离岩石样本中心r处放置磁梯度张量仪1；
[0010]
b、磁梯度张量仪1输出调零；
[0011]
c、将三轴无磁转台3绕z轴旋转45
°
、90
°
、135
°
，再使用磁梯度张量仪1分别测量g
ij
数据；
[0012]
d、将岩石样本绕岩石样本所处的坐标系的x轴旋转90
°
，使其z轴处于水平面内；
[0013]
e、使用磁梯度张量仪1测量g
ij
数据；
[0014]
f、可重复步骤c-e，测量多组数据求取平均值；
[0015]
g、获取磁化率张量数据信息。
[0016]
进一步地，所述岩石样本为正方体型。
[0017]
更进一步地，步骤g，具体是根据公式(7)获取磁化率张量数据信息，公式(7)为磁偶极子型磁梯度的张量表达式，其中，正方体型样本当做磁偶极子型；
[0018][0019]
进一步地，所述岩石样本为圆柱体型。
[0020]
更进一步地，步骤g，具体是根据公式(8)获取磁化率张量数据信息，公式(8)为线偶极子型磁梯度的张量表达式，其中，圆柱体型样本当做线偶极子型；
[0021][0022]
与现有技术相比，本发明的有益效果在于：
[0023]
本发明基于磁梯度张量反演磁化率张量的方法，利用磁梯度张量数据能很好的消除地理位置以及时变磁场的影响，受地理位置影响小以及时变磁场影响小，且能得到完整的磁化率张量信息；无需专门的磁化率张量测量仪器，在实际的测量中更加方便，且测量结果能与专门的磁化率张量仪进行相互验证。
附图说明
[0024]
图1为岩石样本及磁梯度张量仪放置示意图；
[0025]
图2为基于磁梯度张量反演磁化率张量的方法的流程图。
[0026]
图中：1.磁梯度张量仪 2.岩石样本放置处 3.三轴无磁转台 4.x轴旋转电机 5.y轴旋转电机 6.z轴旋转电机。
具体实施方式
[0027]
下面结合实施例对本发明作进一步说明：
[0028]
实施例
[0029]
岩石样本及磁梯度张量仪放置如图1所示，其中，r为磁梯度张量仪与岩石样本中心距离，x、y、z为岩石样本坐标系。
[0030]
如图2所示，本发明基于磁梯度张量反演磁化率张量的方法，包括以下步骤：
[0031]
a、将岩石样本(正方体型或圆柱体型)置于三轴无磁转台3上，在岩石样本所处的x-y-z坐标系下的x轴上距离岩石样本中心r处放置磁梯度张量仪1，其中，r放置的位置满足张量仪的信噪比大于10db(以保证计算结果的相对误差小于5％)，如图1所示；
[0032]
b、校正磁梯度张量仪1，使其输出为零；
[0033]
c、启动三轴无磁转台3绕z轴旋转45
°
，磁梯度张量仪1测量此时岩石样本的磁梯度张量数据；继续使三轴无磁转台3绕z轴旋转45
°
，磁梯度张量仪1测量此时岩石样本的磁梯度张量；再继续使三轴无磁转台3绕z轴旋转45
°
，磁梯度张量仪1测量此时岩石样本的磁梯度张量；
[0034]
d、恢复三轴无磁转台3到初始状态，启动三轴无磁转台3绕y轴旋转90
°
，磁梯度张量仪1测量此时岩石样本的磁梯度张量；
[0035]
e、分别取对应角度测量到的磁梯度张量数据的前三组数据，即g
xx
、g
xy
、g
x
；
[0036]
f、根据公式(17)、(18)计算岩石样本的磁化率张量以及剩余磁化矢量。圆柱体型样本当做线偶极子型，正方体型样本当做磁偶极子型。
[0037]
其中，公式(17)推导过程如下：
[0038]
m＝m
×
v＝(m
i
+m
r
)
×
v
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0039][0040]
其中：m
i
和m
r
分别是感应磁化强度和剩余磁化强度；t为当地地磁总场值；v为岩石样本体积；k为岩石样本磁化率。
[0041]
其磁化率张量如等式(3)所示
[0042][0043]
[0044]
将等式(3)代入等式(4)中，得到：
[0045][0046]
k
xx
、k
xy
、k
xz
、k
yy
、k
yz
、k
zz
是磁化率张量k的六个独立分量；t
x
、t
y
、t
z
分别是当地地磁总场三分量；m
rx
、m
ry
、m
rz
分别是岩石剩余磁化强度三分量；v是岩石样本的体积；μ0是真空磁导率，是一个常量。
[0047]
整理等式(5)可得到如下等式(6)
[0048][0049][0050][0051]
其中，等式(7)为磁偶极子型磁梯度的张量表达式；等式(8)为线偶极子型磁梯度的张量表达式。
[0052]
将等式(6)分别代入等式(7)、(8)中可得到等式(9)、(10)
[0053][0054][0055]
其中，等式(9)为磁偶极子型表达式；等式(10)为线偶极子型表达式。
[0056]
其中：a1＝3x(3r
2-5x2) a2＝3y(r
2-5x2) a3＝3z(r
2-5x2) a4＝3x(r
2-5x2) a5＝15xyz a6＝3x(r
2-5z2) a7＝x(6r
2-8x2) a8＝y(2r
2-8x2) a9＝z(2r
2-8x2) a
10
＝8xyz
[0057]
b1＝3x(r
2-5y2) b2＝3y(3r
2-5y2) b3＝3y(r
2-5y2) b4＝3y(r
2-5z2) b5＝3z(r
2-5y2) b7＝x(2r
2-8y2) b8＝y(6r
2-8y2) b9＝z(2r
2-8y2) b
10
＝x(2r
2-8z2)
[0058][0059]
c1＝3x(4r
2-5x
2-5y2) c2＝3y(4r
2-5x
2-5y2) c3＝3z(2r
2-5x
2-5y2)
[0060]
例如绕z轴旋转45
°
后磁化率张量的表达式可写成如下形式：
[0061][0062]
其旋转后的磁化率张量六个独立分量可写成公式(11.b)
[0063][0064]
如果被测岩石样本是立方体型就将公式(11.b)代入公式(9)；如果被测岩石样本是圆柱体型就将公式(11.b)代入公式(10)，得到磁偶极子型表达式(12.a)、线偶极子型表达式(12.b)。
[0065]
g'＝abk
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(12.a)
[0066]
g'＝a'bk
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(12.b)
[0067]
其中：
[0068][0069]
[0070][0071]
绕z轴旋转一定角度或者绕x轴旋转一定角度只需要把公式(12)中的旋转矩阵b替换掉即可。
[0072]
岩石样品绕z轴旋转90
°
所得到的旋转矩阵b
[0073][0074][0075]
岩石样品绕z轴旋转135
°
所得到的旋转矩阵b
[0076][0077][0078]
岩石样品绕x轴旋转90
°
所得到的旋转矩阵b
[0079]
[0080][0081]
由于，受到每个角度测量到的磁梯度数据组成的矩阵表达式的秩为3的限制，想利用一次测量数据求解9个未知量显然是不能实现的，所以需分别取对应角度测量到的磁梯度张量数据的前三组数据，即gxx、gxy、gxz；组成新的矩阵表达式g”＝a”b'k，新矩阵表达式的秩为9来求解9个未知量。
[0082][0083]
只需对公式(16)进行逆运算即可求得磁化率张量的六个独立分量以及剩余磁化强度的三个分量，如公式(17)。
[0084]
[0085]
因为系数矩阵a”中包含的参数较多，无法一一表达出来，故选择元素进行代替。等式(17)中系数矩阵a”中的元素表达式具体为：
[0086]
a
12
＝d(-a1t
x
+a2t
y
)、)、a
16
＝da3t
z
、a
19
＝fa3、a
22
＝d(-a2t
x
+b1t
y
)、)、a
26
＝-da5t
z
、、a
29
＝fa3、a
32
＝-d(a3t
x
+a5t
y
)、)、a
36
＝da6t
z
、、a
39
＝fa6、a
41
＝da2t
y
、a
42
＝-d(a1t
y
+a2t
x
)、a
43
＝d(a2t
z
+a3t
y
)、a
44
＝da1t
x
、a
45
＝-d(a1t
z
+a3t
x
)、a
46
＝da3t
z
、a
47
＝fa2、a
48
＝-fa1、a
49
＝fa3、a
51
＝db1t
y
、a
52
＝-d(a2t
y
+b1t
x
)、a
53
＝d(b1t
z-a5t
y
)、a
54
＝da2t
x
、a
55
＝-d(a2t
z-a5t
x
)、a
56
＝-da5t
z
、a
57
＝fb1、a
58
＝-fa2、a
59
＝-fa5、a
61
＝-da5t
y
、a
62
＝-d(a3t
y-a5t
x
)、a
63
＝d(a6t
y-a5t
z
)、a
64
＝da3t
x
、a
65
＝-d(a3t
z
+a6t
x
)、a
66
＝da6t
z
、a
67
＝-fa5、a
68
＝-fa3、a
69
＝fa6、a
72
＝d(a1t
x-a2t
y
)、)、a
76
＝da
3 t
z
、a
79
＝fa3、a
82
＝d(a2t
x-b
1 t
y
)、
a
86
＝-da5t
z
、、a
89
＝-fa5、a
92
＝d(a3t
x
+a5t
y
)、)、a
96
＝da6t
z
、a
99
＝fa6、a
101
＝da1t
x
、a
102
＝d(a1t
z
+a3t
x
)、a
103
＝-d(a
1 t
y
+a2t
x
)、a
104
＝da3t
z
、a
105
＝-d(a2t
z
+a3t
y
)、a
106
＝da2t
y
、a
107
＝fa1、a
108
＝fa3、a
109
＝fa2、a
111
＝da2t
x
、a
112
＝d(a2t
z-a5t
x
)、a
113
＝-d(a2t
y
+b1t
x
)、a
114
＝-da5t
z
、a
116
＝db1t
y
、a
117
＝fa2、a
118
＝-fa5、a
119
＝fb1、a
121
＝da3t
x
、a
122
＝d(a
3 t
z
+a6t
x
)、a
123
＝-d(a3t
y-a5t
x
)、a
124
＝da6t
z
、a
125
＝-d(a6t
y-a5t
z
)、a
126
＝-da
5 t
y
、a
127
＝fa3、a
128
＝fa6、a
129
＝fa5。
[0087]
上述是表示岩石样本为立方体(磁偶极子)的反演表达式，如果岩石样本为圆柱体型(线偶极子)，还得用系数矩阵a”'来替换系数矩阵a”。
[0088][0089]
替换后可得到如下等式(18)
[0090][0091]
等式(18)中系数矩阵a”'中的元素表达式具体为：
[0092]
a
12
＝-d1(a7t
x-a8t
y
)、)、a
16
＝da9t
z
、a
19
＝fa9、a
22
＝-d1(a8t
x-b7t
y
)、)、a
26
＝-d1a
10
t
z
、、a
29
＝-f1a
10
、a
32
＝-d(a9t
x
+a
10
t
y
)、)、a
36
＝d1b
10
t
z
、、a
39
＝f1b
10
、a
41
＝d1a8t
y
、a
42
＝-d1(a7t
y
+a8t
x
)、a
43
＝d1(a8t
z
+a9t
y
)、a
44
＝d1a7t
x
、a
45
＝-d1(a7t
z
+a9t
x
)、a
46
＝d1a9t
z
、a
47
＝-f1a7、a
48
＝-f1a7、a
49
＝f1a9、a
51
＝d1b7t
y
、a
52
＝-d(a7t
y
+b7t
x
)、a
53
＝d(b7t
z-a
10
t
y
)、a
54
＝d1a8t
x
、a
55
＝-d(a8t
z-a
10
t
x
)、a
56
＝-d1a
10
t
z
、a
57
＝f1b7、a
58
＝-f1a8、a
59
＝-f1a
10
、a
61
＝-d1a
10
t
y
、a
62
＝-d1(a9t
y-a
10
t
x
)、a
63
＝d1(b
10
t
y-a
10
t
z
)、a
64
＝d1a9t
x
、a
65
＝-d1(a9t
z
+b
10
t
x
)、a
66
＝d1b
10
t
z
、a
67
＝-f1a
10
、a
68
＝-f1a9、a
69
＝f1b
10
、
a
72
＝d1(a7t
x-a8t
y
)、)、a
76
＝d1a
9 t
z
、a
79
＝f1a9、a
82
＝d1(a8t
x-b7t
y
)、)、a
86
＝-d1a
10
t
z
、、a
89
＝-f1a
10
、a
92
＝d1(a9t
x
+a
10
t
y
)、)、a
96
＝d1b
10
t
z
、、a
99
＝f1b
10
、a
101
＝d1a7t
x
、a
102
＝d1(a7t
z
+a9t
x
)、a
103
＝-d1(a7t
y
+a8t
x
)、a
104
＝d1a9t
z
、a
105
＝-d1(a8t
z
+a9t
y
)、a
106
＝d1a8t
y
、a
107
＝f1a7、a
108
＝f1a9、a
109
＝-f1a8、a
111
＝d1a8t
x
、a
112
＝d1(a8t
z-a
10
t
x
)、a
113
＝-d1(a8t
y
+b7t
x
)、a
114
＝-d1a
10
t
z
、a
116
＝d1b7t
y
、a
117
＝f1a8、a
118
＝-f1a
10
、a
119
＝-f1b7、a
121
＝d1a9t
x
、a
122
＝d1(a9t
z
+b
10
t
x
)、a
123
＝-d1(a9t
y-a
10
t
x
)、a
124
＝d1b
10
t
z
、a
125
＝-d1(b6t
y-a
10
t
z
)、a
126
＝-d1a
10
t
y
、a
127
＝f1a9、a
128
＝f1b
10
、a
129
＝f1a
10
。