地球物理学报  2020, Vol. 63 Issue (8): 3037-3049   PDF    
钻孔岩心重定向的古地磁学新方法及其在华南铀矿床岩心定向中的初步应用
葛坤朋1,2, 谢基海2, 应阳根3, 仇登登4, 王宇钦2, 邓居智2     
1. 东华理工大学核资源与环境国家重点实验室, 南昌 330013;
2. 东华理工大学地球物理与测控技术学院, 南昌 330013;
3. 江西省核工业地质局261大队, 江西鹰潭 344400;
4. 广东省核工业地质局293大队, 广州 510800
摘要:基于岩石剩余磁化强度的矢量特性,使用古地磁学方法对钻孔岩心进行重新定向,对于解译地下构造、分析应力特征、判定成矿方向等问题具有重要意义.针对华南地区复杂构造下的倾斜钻孔岩心,本文构建了古地磁学岩心定向新方法理论,结合多次坐标转换和磁场一致性扫描,实现了基于黏滞剩磁的岩心重定向;设计了定向所需的成套古地磁学岩心定向、测量、取样装置,并将该方法和装置应用于华南铀矿钻孔的重定向,恢复后的岩心特征剩磁计算获得的早白垩世古地磁极位置与同期华南地块的极位置一致.该方法对地球物理深部探矿技术的改进有一定意义.此外,所设计的研究装置在古地磁学研究对象的拓展中也具有应用前景.
关键词: 古地磁学      岩心      重定向      坐标转换      剩磁     
A new paleomagnetic method of borehole core reorientation and its preliminary applications in reorientation of cores of uranium deposits in southern China
GE KunPeng1,2, XIE JiHai2, YING YangGen3, QIU DengDeng4, WANG YuQin2, DENG JuZhi2     
1. State Key Laboratory of Nuclear Resources and Environment, East China University of Technology, Nanchang 330013, China;
2. School of Geophysics and Measurement-control Technology, East China University of Technology, Nanchang 330013, China;
3. Jiangxi Province Nuclear Industry Geology Bureau 261 Dadui, Yingtan Jiangxi 344400, China;
4. Guangdong Province Nuclear Industry Geology Bureau 293 Dadui, Guangzhou 510800, China
Abstract: Based on the vector characteristics of the rock remanent magnetization, the reorientation of the borehole core by paleomagnetism method is of great significance for interpreting the underground structure, analyzing the stress characteristics and determining the ore-forming direction. Aiming at the inclined borehole cores under complex structures in south China, a new orientation method for paleomagnetic cores is proposed in this paper, which combines multiple coordinate transformations and consistent magnetic field scanning to realize core reorientation based on viscous remanent magnetization. A complete set of paleomagnetic core orientation, measurement and sampling devices for orientation is designed, and the method is applied to the reorientation of boreholes in south China uranium deposits. The calculated position of paleomagnetic pole in early Cretaceous based on recovered core characteristic remanent magnetization is consistent with that of the South China Plate in the same period. This research has an important significance in improving deep geophysical prospecting technology. In addition, the designed research devices also have application prospect in the development of objects of study in paleomagnetism.
Keywords: Paleomagnetism    Core    Reorientation    Coordinates transformation    Remanence    
0 引言

岩心作为地下地质信息的直接载体,包含裂缝成矿、古应力、古水流方向等重要地质信息(Deng et al., 2013; Ge et al., 2017).然而大部分岩心在从钻孔提取时已经发生旋转,导致其方向性信息无法直接被利用,因此其重定向意义重大.作为恢复板块运动历史、重建全球气候变化和揭示地球内部过程等重大地质问题的主要手段之一(Dunlop, 2002; 刘青松等, 2007; Tauxe, 2010; Liu et al., 2012),以矢量方向特性的地磁场为研究基础的古地磁学也是一种重要的岩心重定向方法(谈晓东等,1994孟小红和周海民, 1997).由于古地磁方法不仅在成本上较摄影工具、钻孔成像技术、倾角仪、全岩心扫描技术等方法低,而且其数据精度也能够得到保证(谢基海等,2019).近年来,随着古地磁学的不断发展,应用古地磁技术对岩心进行重定向的这一方法越来越受到学者的重视(杨振宇和Moreau, 1998; 李学森等, 2006; 韦乐乐等, 2013).

古地磁学岩心定向的基本原理是岩心的黏滞剩磁分量方向与现代地磁场方向在统计意义上一致,而岩心的高温剩磁分量方向与岩石形成或重磁化时的地磁场方向一致(朱岗崑, 2005).古地磁学岩心定向的关键在于:(1)是否能够有效提取剩磁矢量;(2)能否将旋转后的天然剩磁矢量方向恢复至岩心原始状态.通过对岩心样品进行逐步热退磁或交变退磁(Butler, 1992; Tauxe, 2010),可分离提取岩心的低温黏滞剩磁分量和高温特征剩磁分量.通过坐标转换可将岩心样品由样品坐标系校正到地理坐标系,将地理坐标系下的低温黏滞剩磁分量方向与现代地磁场方向对比,高温特征剩磁分量方向与当时地磁场方向进行对比,即可实现岩心定向(杨斌谊等, 2003).因此古地磁学岩心定向最为常用的是黏滞剩磁和特征剩磁校正(侯守信和田国荣,1999),此外也可采用其他的次生剩磁,但对该次生剩磁的准确年龄要求较高(Hailwood and Ding, 1995).此外,韩非等(2017)指出在特定情况下可使用特征剩磁进行岩心重定向,但在构造复杂的地区其应用存在困难.

过去应用古地磁学实现岩心定向的前提是(谈晓冬等,1994朱日祥等,2003周亚楠等,2015韩非等,2017):(1)钻孔近于垂直;(2)地层近于水平.然而实际钻井中存在大量斜井,对于斜井岩心定向,除了李学森等(2006)水平取样下依据等面积图内D、I扫描推导了斜井岩心定向算法和刘敬寿等(2015)通过岩层产状法推导了斜井岩心坐标转换公式,其余理论研究较少.并且前人主要集中于对华北东北地区的钻孔开展研究,而对岩浆岩变质强烈、热液蚀变机理复杂的华南地区,古地磁学定向方法鲜有研究.使得该地区岩心裂隙方位和地下成矿构造属性的解译存在困难,一定程度降低了华南地区固体矿产的勘探效率(Min et al., 2005).

此外,由于古地磁学设备主要针对地表露头采样而设计,对于普遍较细(直径约30 mm)的华南地区岩心古地磁样品的精确采集,目前古地磁设备所使用的钻机、太阳罗盘、定向仪及相关处理软件均无法适用.现在处理岩心磁学样品时,多使用非专业古地磁设备,如台钻、线切割机等,其加工效率和误差无法估计,影响了古地磁数据的准确采集.再者,因微古地磁系统设备的缺乏,对岩心古地磁样品的采集被限制于对垂直钻孔的定向采集(李学森等,2006韩非等,2017),而对华南地区普遍的多裂隙、强变质岩心样品仍然难以处理.因此亟须设计一套专用微古地磁设备系统,能够满足华南地区岩心古地磁学样品处理需求.

综上,本研究即在古地磁学基本理论框架下,在前人无定向岩心的剩磁矢量恢复基础上,基于相关古地磁学剩磁理论,设计制作相关设备,形成一种无定向岩心的剩磁恢复进而重新定向的新方法;并以华南矿田岩心样品为研究对象,对研究方法进行检验;最终对这种古地磁学定向新方法的进一步提升和拓展应用提出展望.

1 理论基础 1.1 数学基础

磁力仪对岩石样品的测量是基于偶极子场形态测量的,在远离坐标原点一定距离时,圆柱体古地磁样品磁场接近磁偶极子场的程度,仅取决于圆柱体的底面直径与柱高之比(葛坤朋等,2011; Collinson, 2013).因此等比例缩小柱体体积,将不会影响磁偶极子场测量,唯一问题在于缩小后的柱体体积总磁矩是否可以被准确测量.随着古地磁仪器的进一步发展,超导磁力仪的灵敏度已经达到10-12 Am2(朱日祥等,2003),足以满足钻孔古地磁小样品的测量需求.

考虑到黏滞剩磁一般为阻挡温度300 ℃以下剩磁组分,相当于近十万年来地磁场平均作用方向(Butler, 1992),故黏滞剩磁矢量方向在统计上应与现今平均地磁场方向(地理北极)一致,而特征剩磁在统计上应与该地区同时代古地磁极对应磁场方位一致,因此可通过对比黏滞剩磁或特征剩磁矢量及其参考方向,将磁参数转换到地理坐标系上(图 1).具体包括:(ⅰ)微古地磁样品坐标系到固定测量仪坐标系的转换;(ⅱ)固定测量仪坐标系到岩心坐标系的转换;(ⅲ)岩心坐标系的旋转变换(搜索旋转角度ω);(ⅳ)旋转岩心坐标系到地理坐标系的转换.

图 1 古地磁小样品和岩心样品坐标系,包括古地磁小样品(a)和小样品在岩心样品中坐标系(b)以及岩心样品在地层下的坐标系(c).其中(a)到(b)为小样品坐标系到固定测量仪的转换; (b)到(c)为固定测量仪到岩心坐标系的转换; (c)的旋转扫描为岩心坐标系的旋转变换; 旋转后(c)的地理改正为岩心坐标系到地理坐标系的转换 Fig. 1 Coordinate system of small paleomagnetic samples and core samples, including coordinate system (a) and (b) of small paleomagnetic samples in core samples and coordinate system (c) of core samples in stratum. Where, (a) to (b) is the conversion of small sample coordinate system to fixed measuring instrument; (b) to (c) is the transformation from the fixed measuring instrument to the core coordinate system; The rotation scan is the rotation transformation of the core coordinate system; The geographical correction of (c) after rotation is the transformation from the core coordinate system to the geographical coordinate system

设(X0, Y0, Z0)是微古地磁样品坐标系,样品方位角为φ,样品倾角为θ,(X1, Y1, Z1)是固定测量仪坐标系,基于两者的关系有

(1)

设(X2, Y2, Z2)是岩心坐标系,其与(X1, Y1, Z1)的关系有

(2)

设(X3, Y3, Z3)是旋转后的岩心坐标系,旋转角度为ω,其与(X2, Y2, Z2)的关系有

(3)

设(X4, Y4, Z4)是地理坐标系,钻孔测斜方位角为Φ,测斜倾角为Θ,其与(X3, Y3, Z3)的关系有

(4)

通过以上四次坐标变换,得到地理坐标系下磁化强度矢量M4,单位矢量为EM4,设参考地磁场磁化强度单位矢量为Eref.则通过检索旋转角度ω,求取min[(EM4-Eref)2]时的最佳旋转角度ωp即可获得最接近参考地磁场方向的磁化强度矢量(特征剩磁或黏滞剩磁)的旋转角度.此时,将原始岩心样品按照标志线顺时针旋转ωp,即将岩心恢复到了原位.

需要注意的是,由于早期钻孔岩心经过搬运震动或者实验时错误标记可能导致岩心反向放置,影响岩心的复位和对剩磁方向的判断.判断岩心是否反向放置,可考虑黏滞剩磁的磁化方向最小原理,即磁化强度方向倾向于沿着磁场方向磁化,原位黏滞剩磁与地磁场方向不大于90°的可能性最大.此时可以考虑第(ⅲ)次坐标变换即旋转扫描时,360°扫描情况下参考矢量M3的均值M3与参考场Mref的夹角是否为钝角,如果〈M3, Mref〉大于90°,则判定为反向放置.此时有两种复位方式:其一求取max[(EM4-Eref)2]的角度ωp.实际操作中需要将岩心按原心对称调转180°,并将原始岩心样品按照标志线逆时针旋转ωp,即将岩心恢复到了原位;其二是在第二次坐标变换向量M2变更为-M2,依次进行第(ⅲ)、(ⅳ)次坐标变换,搜寻min[(EM4-Eref)2]情况下的ωp,实际操作中将岩心按原心对称调转180°,并将原始岩心样品按照标志线顺时针旋转ωp,同样可实现岩心复位.

对于原始特征分量ChRM的获取,可将其三分量按照上述进行第(ⅰ)和(ⅱ)次坐标变换,随即旋转ωp角度后,进行地理改正后的结果即为原位特征分量.

1.2 装置设计

在使用水冷充电钻机搭配无/弱磁微钻头对野外露头进钻完毕后,需要对微古地磁样品进行定向和取样.由于古地磁专用太阳罗盘设备仅适用于标准古地磁取样,因此需要设计太阳罗盘连接器(图 2),搭配太阳罗盘进行样品定向.在使用铜丝刻画标志线后,使用设计制作的无磁取样器(如铜制)进行取样.取出来的样品,可直接用作一些古地磁学实验(如地磁古强度)测量;或定向放入自制标准石英盒(方形盒和柱状盒)中,形成可用于古地磁学标准实验的样品.

图 2 太阳罗盘连接器(a)、取样器(b)及岩心固定测量仪(c)示意图,其中各单元部件分别为: 1.太阳罗盘对接线; 2.微钻孔套管; 3.标志线刻画孔; 4.固定套; 5.太阳罗盘套管; 6.标准样品模型; 7.太阳罗盘插槽; 8.取样器壁; 9.取样器套管; 10.支撑板; 11.旋转式水平调节件; 12.下卡槽; 13.下底座; 14.上卡槽; 15.方位标志线; 16.上底座; 17.岩心基准线; 18.螺纹孔 Fig. 2 Schematic diagram of sun compass connector (a), sampler (b) and core fixed measuring instrument (c), in which each unit parts. 1. Respectively the sun compass to terminal; 2. Micro drilling casing; 3. Marker depicting hole; 4. Fixed set; 5. Sun compass sleeve; 6. Standard sample model; 7. The sun compass slot; 8. Wall of the sampler; 9. Sampler casing; 10. Support plate; 11. Rotating horizontal adjustment; 12. Under slot; 13. Under base; 14. Upper slot; 15. Mark line of direction; 16. Upper base; 17. Core reference line; 18. Threaded hole

在本项目中,因样品采集于无定向的岩心上,因此需要将岩心样品按一定方向放置,并标记用于恢复岩心方位的预旋转标志线,才能逐步恢复其原始方位.通过设计制作岩心微古地磁固定测量仪(图 2),将岩心样品固定于测量仪上卡槽5和下卡槽9中,保证方向标志线10指向北方向,岩心基准线8与岩石预旋转标志线对齐,随后进行微古地磁定向和采集.

2 实际应用

因研究获得的方法对处理复杂、破碎地区的岩心钻孔具有重要应用价值,因此本文在华南铀矿床地区进行了应用测试.

2.1 华南铀矿床地区地质简介及研究对象

华南铀成矿省是我国脉型铀矿分布最密集的铀矿省.在该成矿省江西、广东境内,以及赣粤、赣湘粤、赣闽、湘粤交界处已发现多处火山岩型、花岗岩型铀矿田,如相山矿田、下庄矿田等(Hu et al., 2008).其中下庄铀矿位处于贵东岩体东部,华夏古陆西缘,岩性以燕山期侵入的中粒黑云母花岗岩为主体,其次为边缘相的白云母花岗岩(图 3).区域内构造亦十分发育,断裂带呈NNE向展布,下庄矿田铀矿床的产出即受到该断裂带控制.该地区中基性岩脉极为发育,广泛出露后期侵入的辉绿岩,呈NWW方向展布.辉绿岩脉、硅化断裂带、花岗岩围岩在热液作用下,最终形成了以辉绿岩脉和硅化断裂带交点成矿的“交点型”为主要类型的铀矿床(Ge et al., 2017).

图 3 下庄地区地质构造示意图及取样钻孔位置(修改自Ge et al., 2017). 图中简称NCB:华北地块;YB:扬子地块;CB:华夏地块 Fig. 3 Geological sketch map of Xiazhuang uranium orefield and the locations of sampling cores (modified after Ge et al. 2017). Abbreviated symbols: NCB: North China block; YB: Yangtze block, CB: Cathaysian block

本研究陆续采集了下庄地区湖子北、对山、小水等地区8个“交点型”铀矿钻孔的辉绿岩脉岩心样品64块(图 4),编号为J1-8,随即使用本文研究装置定向采集加工为80块圆柱状古地磁学小样品.根据偶极子效应(Collinson,2013),保证高径比等同于标准样品(~0.8).经线切割机处理后,小样品尺寸为高~0.8 cm,底面直径~1 cm,随即装入标准样品陶瓷盒子(高2 cm,底面直径2.54 cm),形成标准样品.同时准备一组平行样品研磨成粉末,进行相关岩石磁学测试.

图 4 预采集下庄岩心样品及部分在研器件的使用图.即使用实验室台钻(a)在水中对岩心样品进行进钻(b),对钻取的古地磁小样品(c)进行原位岩心定向(d),随后放入托盘进行系统标记(e),最后置于标准样品的陶瓷盒子等待测量(f) Fig. 4 Pre-collected and in-process core samples of Xiazhuang using developed devices. Namely, the laboratory bench drill (a) was used to drill into the core samples in water (b), the drilled paleomagnetic small samples (c) were oriented in situ (d), and then put into the tray for system marking (e), and finally placed in the ceramic box of the standard samples for measurement (f)
2.2 岩石磁学方法

室温岩石磁学测量使用了中科院地质与地球物理研究所古地磁与年代学实验室VSM3900振动样品磁力仪进行.磁滞回线和一阶反转曲线(FORC)的测量设置最大外加场为1~1.5 T,可同时获得Mrs, Ms, Bc等磁滞参数.随即测量反向场退磁曲线,以获得剩磁矫顽力Bcr.

样品的磁化率(χ-T,质量归一化)随温度变化曲线应用MFK-FA磁化率仪及其CS-4温度控制系统测量.所有样品的最高加热温度700 ℃,温度梯度为2 ℃·min-1,为避免样品在加热过程中氧化,测量均在氩气环境中进行,氩气流量为50 mL·min-1.

三轴热退磁实验使用2G660脉冲磁化仪沿样品Z,YX轴分别施加2.5 T,0.5 T和0.15 T磁场,然后使用MMTDSC热退磁炉进行逐步热退磁,并在每步退磁后使用2G760超导磁力仪进行测量.

利用MPMS 5XL低温磁性测量系统对样品进行低温测量.首先在零场中降温到10 K,然后使样品在2.5 T的外场下获得SIRM.逐渐升温到300 K,升温梯度为5 K·min-1,测量低温下饱和剩磁随温度的变化.

2.3 古地磁学方法

古地磁学实验使用古地磁与年代学实验室TD-48热退磁炉进行逐步热退磁.TD-48炉内残余磁场在10 nT以下,温度控制误差小于±2 ℃.退磁采用热退磁方法,退磁温度段设计为NRM,60,80,100,120,150,200,250,300,350,400,450,500,525,560,570,575,585 ℃,随即使用2G755超导磁力仪进行磁化强度矢量的实验测量.获得数据使用主成分分析法进行运算,获取初始的MDI等参数.

3 实验结果 3.1 岩石磁学结果

样品常温磁学特性如图 5所示,样品J3-1磁滞回线具有一定顺磁特征,顺磁校正后显示出假单畴(PSD)或多畴(MD)行为,其FORC曲线显示相互作用场展布较宽,呈现近MD特征.而样品J6-5含有更多的顺磁组分,其铁磁组分的饱和磁化强度较低,但显示出更高的矫顽力Bc和剩磁比Mrs/Ms.其FORC曲线能见到更多同心轮廓,相互作用场仍有展布,整体呈现出细颗粒PSD特征.

图 5 典型样品磁滞回线和FORC图 Fig. 5 Hysteresis loops and FORC diagrams of typical samples

典型样品的高低温磁学实验如图 6所示,三轴退磁实验显示样品J3-1具有明显的三种磁性矿物组分,其中软磁性组分载磁最高,在约600 ℃磁性退尽.磁化率加热曲线,在升温初期曲线平直,样品在580 ℃时磁化率出现明显下降,具有MD特征.低温SIRM获得曲线显示120 K左右的Verwey转化,证实磁铁矿为主要软磁性载磁组分.对于样品J6-5,仍以软磁性组分载磁最高,但该组分在300~400 ℃快速降低,显示主要载磁矿物可能为磁黄铁矿等铁硫化物.热磁曲线和SIRM获得曲线显示样品仍含有少量磁铁矿.

图 6 典型样品的三轴热退磁、χ-T和低温SIRM获得曲线图 Fig. 6 Thermal demagnetization of three-component IRM, χ-T and low temperature SIRM acquirement curves of typical samples

综上岩石磁学结果,可以认为所采集辉绿岩样品主要分为两类,第一,以MD磁铁矿颗粒为主要载磁矿物的新鲜辉绿岩样品,其可能包含赤铁矿等硬磁组分,样品具有较强的铁磁性特征;第二,以磁黄铁矿为主要载磁矿物的弱蚀变辉绿岩样品,并含有少量PSD磁铁矿颗粒,样品的超顺磁性可能来源于黄铁矿化的作用.

3.2 古地磁学结果

经过系统的热退磁实验,可获得岩心样品天然剩磁的热退磁曲线.在筛除了主成分分析法MAD大于20°的磁学数据后,共计获得了56个可以提取黏滞剩磁和特征剩磁的样品.典型热退磁曲线如图 7所示,Z氏图主要显示低温和中高温两个分量.样品在300 ℃以下可以分离出较稳定的黏滞剩磁,其强度可达到天然剩磁总强度的20%~40%.300 ℃以上至500~585 ℃可以分离出中高温分量.在基于黏滞剩磁的坐标转换校正后,古地磁偏角明显变小,接近现代场平均方向,但古地磁视倾角变化并不显著.

图 7 典型岩心样品坐标校正前(a-f)后(a′-f′)的热退磁正交投影图,空心圆形(正方形)分别代表磁化强度矢量在垂直XZ面(水平XY面)的投影,虚线箭头代表VRM的方向 Fig. 7 The orthogonal projections of thermal demagnetization before (a-f) and after (a′-f′) coordinate correction of typical samples. The hollow circle (square) represents the projections of magnetization vector on the vertical XZ plane (horizontal XY plane), and the dotted arrow represents the directions of VRM

黏滞剩磁在进行坐标转换校正前,在等面积投影图呈现较分散的状态,但并非完全分散(图 8).这是由于当测斜倾角为90°,即岩心直立时,倾角应统一接近于现代场角度.而实际测斜倾角较高,为~75°,因此岩心取出时已经发生无定向旋转,但仍有一定集中性.在基于现代场的坐标转换后,黏滞剩磁实现了与现代场最为接近的旋转,改正后的磁偏角接近0°,而倾角则在现代场倾角附近均匀分散.特征剩磁在坐标转换校正前,同样呈现较分散的状态.在进行基于黏滞剩磁方向的坐标转换后,置信圆半径明显变化,显示数据一定程度的聚集,并且统计平均值偏离了原来的方位.

图 8 黏滞剩磁和特征剩磁校正前(a, c)和校正后(b, d)的等面积投影图 Fig. 8 Equal area projections of VRM and ChRM before (a, c) and after (b, d) coordinate correction

基于坐标旋转,将获得的校正黏滞剩磁进行Fisher统计平均后,转换为VGP经纬度为83.6°N,322.2°E,A95 (, Khramov, 1987)为3.72°,其与参考地磁极90°N,0°E接近(图 9).基于李献华等(1997)邹东风等(2011)对下庄-寨下辉绿岩进行Ar-Ar年代学定年研究,获得该区域辉绿岩具体年龄为109.9±2.3 Ma,为早白垩世末期侵入岩体.因此我们将Fisher统计后的校正特征剩磁DI值换算为VGP,与早白垩世扬子地块已有文献的VGP进行了比较(Besse and Courtillot, 1991; Enkin et al., 1992吴汉宁等, 1998黄宝春等,2008).实际统计结果与文献统计VGP结果接近,并落在全部文献VGP置信区域内部,进一步证实了岩心定向新方法的准确性和在华南铀矿床岩心定向的可行性.

图 9 特征剩磁ChRM经坐标变换后求取的VGP(五星)以及文献报道早白垩世VGP及其置信圆(A95) Fig. 9 Equal area projections of VGP (five-star) obtained by the ChRM after the coordinate transformation, and the referred VGPs in early cretaceous and the corresponding confidence circles (A95)
4 讨论与展望 4.1 岩心重定向的古地磁学新方法及意义

本重定向方法对任意岩心的定向须结合所设计装置、古地磁学实验和基于1.1节的Matlab数据处理程序协作完成.具体的步骤为:(ⅰ)沿钻孔岩心柱面标记岩心标志线,并固定于岩心固定测量仪上(图 2c)进行台钻或电钻进钻,使用太阳罗盘结合连接器(图 2a)进行定向并用铜丝刻线,随即用岩心小样品取样器(图 2b)取出;(ⅱ)根据岩石磁学鉴定样品载磁矿物的主要特征,将小样品套入标准样品盒,使用超导磁力仪进行标准古地磁学实验,获取热退磁曲线,使用主成分分析法提取黏滞剩磁(Dv, Iv)和特征剩磁方位(Dc, Ic);(ⅲ)编制包括钻孔测斜方位角和倾角、岩心小样品坐标系方位角和倾角、黏滞剩磁和特征剩磁DI方位的基本数据表;(ⅳ)使用本文编制的Matlab数据处理程序调用数据文件,进行基于现代场方位的四次坐标变换,可获得钻孔岩心标志线的应旋转角度ωp,实现对倾斜岩心的重定向.

本定向方法的优势在于:(ⅰ)能够处理垂直岩心和倾斜岩心的重定向问题;(ⅱ)能够实现岩心任意位置的小样品钻取定向和取样,可以有效规避裂隙,对复杂构造下的破裂岩心重定向具有显著的优势;(ⅲ)能够处理特殊岩心样品,如丢失顶底方向的老岩心,可以通过古地磁学方法判断并恢复方位.特别是对于早期缺乏精确钻井成像技术的钻井,并且钻孔已经填埋或塌陷的情况,能够对老岩心样品进行重定向.而本方法局限性仅在于古地磁学本身的精确性,特别是需要进行精细岩石磁学对岩心样品的载磁性质进行判定,对黏滞剩磁和特征剩磁方向进行准确提取.

岩心定向后的古地磁学继续研究对于华南地区铀矿田勘探研究具有重要意义.将特征剩磁恢复至岩心原位方位,可对地下岩心样品进行传统古地磁学构造、应力展布、受热历史研究(Butler, 1998; Parés et al., 2007),能够拓展华南地区古地磁学研究对象,加深对华南地区的构造成矿认识;通过精确的岩石剩磁测试,可得到华南铀矿田天然剩磁特征,厘定富矿构造的磁学属性,进而能够对磁法勘探异常进行正确的解释和反演研究;同时,结合岩心裂隙相对位置,可进一步确定岩心裂缝方向和破碎带走向,锁定铀矿成矿方位.

因此,在将来使用钻孔岩心重定向的古地磁学新方法及其后续传统古地磁学方法,通过对华南地区丰富钻孔岩心样品和资料的再研究,不仅能够圈定华南铀矿田深部成矿构造,而且在提高勘探效率的同时降低了铀矿深部找矿的成本,对地球物理深部探矿技术的改进有一定意义.

4.2 其他应用展望

由于古地磁学取样具有一定的通用性,开发出的微古地磁设备系统将同时适用于一般露头微古地磁学取样(图 2a-b).将来的研究将进一步完善装置,匹配Matlab处理程序,形成与标准古地磁学样品及数据处理相对应的微古地磁样品数据采集和处理系统.这对于环境保护地区的样品(考古遗址、地质公园等)和难获取的样品(岩心样品、缺水地区样品、破碎样品、贵重样品等)的古地磁学、古强度学、岩石磁学、环境磁学、生物磁学定向取样都具有适用性(图 10).因此本研究在古地磁学设备和技术手段的改善上具有重要应用前景.

图 10 实验研制的太阳罗盘连接器的应用及一般野外微古地磁取样.(a-b)使用手持钻机对野外露头进行原位进钻. 因为钻样尺寸很小,可在不平整露头(c)进行成功取样,最后使用太阳罗盘连接器对接太阳罗盘进行定向取样(d) Fig. 10 Application of the solar compass connector in general field micro-paleomagnetic sampling. (a, b) In-situ drilling of the field outcrop with a hand drill, because the size of the drill sample is small, the sample can be successfully sampled on the uneven outcrop (c) and finally the solar compass connector is used to conduct directional sampling with the solar compass (d)
5 结论

本文通过提出一种基于古地磁学的斜井岩心自由定向的数据转换方法,编写成Matlab数据处理程序,配合设计相应的钻孔、定向和取样装置,实现了基于黏滞剩磁扫描的古地磁学岩心重定向.通过早白垩世古地磁极的一致性验证,将该新方法在华南铀矿地区得到了成功的初步应用,为进一步的钻孔磁学研究和华南铀矿田的深部找矿技术开发提供了新思路.结合装置在古地磁学中的适用性,本方法对古地磁学设备和技术手段的改善及其研究对象的拓展上具有重要应用前景.

致谢  感谢广东核工业地质局293大队汪洋、陈渝罗和江西核工业地质局261大队冷建华、林福江对地质勘查取样的帮助.作者与中科院地质与地球物理研究所秦华峰、刘双迟在装置设计和磁学实验进行了详细探讨,同中国地质大学徐慧茹博士在古地磁学计算上进行了有益讨论,在此一并表示感谢.并特别感谢潘永信老师和两位匿名审稿人对本文的指正和建议.
References
Besse J, Courtillot V. Revised and synthetic apparent polar wander paths of the African, Eurasian, North American and Indian Plates, and true polar wander since 200 Ma. Journal of Geophysical Research: Solid Earth, 96: 4029-4050.
Butler R F. 1992. Paleomagnetism:Magnetic Domains to Geologic Terranes. Boston: Blackwell Scientific Publications.
Collinson D W. 2013. Methods in Rock Magnetism and Palaeomagnetism:Techniques and Instrumentation. London: Springer Science & Business Media.
Deng C L, He H, Pan Y, et al. 2013. Chronology of the terrestrial Upper Cretaceous in the Songliao Basin, northeast Asia. Palaeogeography, Palaeoclimatology, Palaeoecology, 385: 44-54. DOI:10.1016/j.palaeo.2012.07.028
Dunlop D J. 2002. Theory and application of the Day plot (Mrs/Ms versus Hcr/Hc) 1. Theoretical curves and tests using titanomagnetite data. TheoJournal of Geophysical Research:Solid Earth, 107(B3): EPM 4-1-EPM 4-22. DOI:10.1029/2001JB000487
Ge K P, Liu Q S, Zhu R X. 2011. The secondary magnetic field generated by magnetic samples in the magnetic shielding room. Progress in Geophysics (in Chinese), 26(3): 843-849. DOI:10.3969/j.issn.1004-2903.2011.03.009
Ge K P, Liu Q S, Deng J Z, et al. 2017. Rock magnetic investigation and its geological significance for vein-type uranium deposits in southern China. Geochemistry, Geophysics, Geosystems, 18(4): 1333-1349. DOI:10.1002/2016GC006756
Enkin R J, Yang Z Y, Chen Y, et al. 1992. Paleomagnetic constrains on the geodynamic history of the major blocks of China from the Permian to the present. Journal of Geophysical Research, 97: 13953-13989. DOI:10.1029/92JB00648
Hailwood E A, Ding F. 1995. Palaeomagnetic reorientation of cores and the magnetic fabric of hydrocarbon reservoir sands. Geological Society, London, Special Publications, 98(1): 245-258. DOI:10.1144/GSL.SP.1995.098.01.15
Han F, Qin H F, Ji Q, et al. 2017. Rock magnetism and paleomagnetism of Late Cretaceous sedimentary rocks in borehole CCSD-LK-I (Shandong province) and implications for the feasibility of core orientation recovery using natural remanent magnetization. Chinese Journal of Geophysics (in Chinese), 60(12): 4730-4740. DOI:10.6038/cjg20171215
Hou S X, Tian G R. 2000. The application of viscous remanent magnetization in orientation determination by using drill cores. Chinese Journal of Rock Mechanics and Engineering (in Chinese), 19(S1): 1128-1131.
Hu R Z, Bi X W, Zhou M F, et al. 2008. Uranium metallogenesis in South China and its relationship to crustal extension during the Cretaceous to Tertiary. Economic Geology, 103(3): 583-598. DOI:10.2113/gsecongeo.103.3.583
Huang B C, Zhou Y X, Zhu R X. 2008. Discussions on Phanerozoic evolution and formation of continental China, based on paleomagnetic studies. Earth Science Fontiers, 15(3): 348-359.
Khramov A N. 1987. Paleomagnetology. Grudina V P, Trans. Berlin Heidelberg: Springer-Verlag.
Li X H, Hu R Z, Rao B. 1997. Geochronology and geochemistry of cretaceous mafic dikes from northern Guangdong, SE China. Geochimica (in Chinese), 26(2): 14-31.
Li X S, Wu H N, Zhang X H, et al. 2006. The paleomagnetic method applied to reorient core and its reliability. Petroleum Geophysics (in Chinese), 4(2): 27-32.
Liu J S, Dai J S, Wang K, et al. 2015. An approach to correct the core fracture attitude in deviated boreholes and its application. Acta Petrolei Sinica (in Chinese), 36(1): 67-73.
Liu Q S, Deng C L, Pan Y X. 2007. Temperature-dependency and frequency-dependency of magnetic susceptibility of magnetite and maghemite and their significance for environmental magnetism. Quaternary Sciences (in Chinese), 27(6): 955-962.
Liu Q S, Roberts A P, Larrasoaña J C, et al. 2012. Environmental magnetism:principles and applications. Reviews of Geophysics, 50(4): RG4002. DOI:10.1029/2012RG000393
Meng X H, Zhou H M. 1997. Magneticstudy of drill-cores in Jidong oil field. Acta Geophysica Sinica (in Chinese), 40(6): 809-820.
Min M Z, Xu H F, Chen J, et al. 2005. Evidence of uranium biomineralization in sandstone-hosted roll-front uranium deposits, northwestern China. Ore Geology Reviews, 26(3-4): 198-206. DOI:10.1016/j.oregeorev.2004.10.003
Parés J M, Hassold N J C, Rea D K, et al. 2007. Paleocurrent directions from paleomagnetic reorientation of magnetic fabrics in deep-sea sediments at the Antarctic Peninsula Pacific margin (ODP Sites 1095, 1101). Marine Geology, 242(4): 261-269.
Tan X D, Wang P Y, Jiang L P. 1994. Progress in techniques for borehole paleomagnetic research. Progress in Geophysics (in Chinese), 9(3): 98-103.
Tauxe L. 2010. Essentials of Paleomagnetism. California: University of California Press.
Wei L L, Cheng X, Liu X T, et al. 2013. A paleomagnetical method to redirect cracks in cores. Chinese Journal of Geology (in Chinese), 48(4): 1286-1294. DOI:10.3969/j.issn.0563-5020.2013.04.023
Wu H N, Zhu Rixiang, Bai L X, et al. 1998. Revised apparent polar wander path of the Yangtze Block and its tectonic implications. Science in China Series D: Earth Sciences, 1998, 41: 78-90.
Xie J H, Ge K P, Xu H R et al. 2019. Review of paleomagnetic core orientation method. Progress in Geophysics (in Chinese), http://kns.cnki.net/kcms/detail/11.2982.P.20190926.1350.010.html.
Yang B Y, Gong J J, Chen J J. 2003. Paleomagnetism research on orientation of drilling cores in oil field. Northwestern Geology (in Chinese), 36(4): 79-83.
Yang Z Y, Moreau M G. 1998. Paleomagnetic study of the Montcornet drilled core. Acta Geophysica Sinica (in Chinese), 41(5): 652-657.
Zhou Y N, Cheng X, Ma L, et al. 2015. The principles and methods on determining spatial orientation of sandstones and fractures in subsurface reservoir. Progress in Geophysics (in Chinese), 30(3): 1243-1250. DOI:10.6038/pg20150333
Zhu G K. 2005. Paleomagnetism-Foundation, Theory, Method, Achievement and Application (in Chinese). Beijing: Science Press.
Zhu R X, Huang B C, Pan Y X, et al. 2003. A brief guide to the laboratory of rock magnetism and paleomagnetism at the institute of geology and geophysics, Chinese academy of sciences. Progress in Geophysics (in Chinese), 18(2): 177-181.
Zou D F, Li F L, Zhang S, et al. 2011. Timing of No. 335 ore deposit in Xiazhuang uranium orefield, northern Guangdong province:evidence from LA-ICP-MS U-PB dating of pitchblende.Mineral Deposits (in Chinese), 30(5): 912-922. DOI:10.1007/s12583-011-0163-z
葛坤朋, 刘青松, 朱日祥. 2011. 磁屏蔽室内磁性样品所产生干扰磁场的探讨. 地球物理学进展, 26(3): 843-849. DOI:10.3969/j.issn.1004-2903.2011.03.009
韩非, 秦华峰, 季强, 等. 2017. 鲁科一井上白垩统沉积岩岩石磁学和天然剩磁分析:兼论利用剩磁方向恢复钻孔岩芯原始方位的可行性. 地球物理学报, 60(12): 4730-4740. DOI:10.6038/cjg20171215
侯守信, 田国荣. 2000. 粘滞剩磁(VRM)岩芯定向的应用. 岩石力学与工程学报, 19(S1): 1128-1131.
黄宝春, 周烑秀, 朱日祥. 2008. 从古地磁研究看中国大陆形成与演化过程. 地学前缘, 15(3): 348-359. DOI:10.3321/j.issn:1005-2321.2008.03.031
李献华, 胡瑞忠, 饶冰. 1997. 粤北白垩纪基性岩脉的年代学和地球化学. 地球化学, 26(2): 14-31. DOI:10.3321/j.issn:0379-1726.1997.02.004
李学森, 吴汉宁, 张小浩, 等. 2006. 钻井岩心重定向的古地磁方法及其可靠性研究. 油气地球物理, 4(2): 27-32.
刘敬寿, 戴俊生, 王珂, 等. 2015. 斜井岩心裂缝产状校正方法及其应用. 石油学报, 36(1): 67-73. DOI:10.3969/j.issn.1001-8719.2015.01.011
刘青松, 邓成龙, 潘永信. 2007. 磁铁矿和磁赤铁矿磁化率的温度和频率特性及其环境磁学意义. 第四纪研究, 27(6): 955-962. DOI:10.3321/j.issn:1001-7410.2007.06.010
孟小红, 周海民. 1997. 冀东油田钻井岩芯的磁学研究. 地球物理学报, 40(6): 809-820. DOI:10.3321/j.issn:0001-5733.1997.06.010
谈晓冬, 王朋岩, 姜莉萍. 1994. 钻孔岩芯古地磁研究方法进展. 地球物理学进展, 9(3): 98-103.
韦乐乐, 程鑫, 刘秀婷, 等. 2013. 岩心裂缝重定向的古地磁方法研究. 地质科学, 48(4): 1286-1294. DOI:10.3969/j.issn.0563-5020.2013.04.023
吴汉宁, 朱日祥, 白立新, 等. 1998. 扬子地块显生宙古地磁视极移曲线及地块运动特征. 中国科学:地球科学, 28(S1): 69-78.
谢基海, 葛坤朋, 徐慧茹, 邓居智.古地磁学岩心定向方法回顾.地球物理学进展. http://kns.cnki.net/kcms/detail/11.2982.P.20190926.1350.010.html.
杨斌谊, 龚建军, 陈建军. 2003. 确定油田钻井岩心原始方位的古地磁学方法探讨. 西北地质, 36(4): 79-83. DOI:10.3969/j.issn.1009-6248.2003.04.014
杨振宇, Moreau M G. 1998. Montcornet钻孔岩芯的古地磁研究. 地球物理学报, 41(5): 652-657. DOI:10.3321/j.issn:0001-5733.1998.05.008
周亚楠, 程鑫, 马轮, 等. 2015. 确定地下埋藏砂体与裂缝空间方位的原理和方法. 地球物理学进展, 30(3): 1243-1250. DOI:10.6038/pg20150333
朱岗崑. 2005. 古地磁学:基础, 原理, 方法, 成果与应用. 北京: 科学出版社.
朱日祥, 黄宝春, 潘永信, 等. 2003. 岩石磁学与古地磁实验室简介. 地球物理学进展, 18(2): 177-181. DOI:10.3969/j.issn.1004-2903.2003.02.001
邹东风, 李方林, 张爽, 等. 2011. 粤北下庄335矿床成矿时代的厘定--来自LA-ICP-MS沥青铀矿U-Pb年龄的制约. 矿床地质, 30(5): 912-922. DOI:10.3969/j.issn.0258-7106.2011.05.012