The processing effect comparison and optimization of small subdomain filtering method in potential field data
-
摘要:
重磁位场数据处理解释的主要任务之一是研究推断地质体的边界位置，而通过对重磁位场数据进行一定方法滤波处理后的相关图件可直观显示地质体的边界。小子域滤波法是众多滤波方法中的一种，前人利用小子域滤波法识别地质体边界位置有较多的研究，目前已有10种类型的小子域滤波法，但是鲜有对10种方法滤波效果的系统研究，给使用者带来诸多不便。因此，设计2组理论模型对比了10种方法的处理效果，并根据试验分析优选了最佳的滤波参数。理论模型与实际资料处理结果表明，滑动形子域滤波法对地质体边界增强的效果最好，且滤波窗口选择5×5数据点、迭代1次即可；对称形子域滤波法在滤除随机噪声的同时能够较好的保持原始异常场的梯度变化特征，滤波窗口选择5×5数据点、迭代2次较好。
Abstract:The study of the edge location of the geological body is the most important task in the processing and interpretation of gravity and magnetic potential field data. The results of relevant map can show the geological body's edge location directly by means of filtering processing with the gravity and magnetic potential field data. Small subdomain filtering is one of the filtering methods, and it has been used by many researchers to recognize the edge location of the geological body. There are 10 kinds of small subdomain filtering methods at present; however, rare studies have analyzed the filtering effect about those methods systematically, and this brings much inconvenience to users. Therefore, the authors designed 2 theoretical models to compare their filtering results and then selected the optimal filtering parameters. The results of theoretical model and field data show that the moving shape small subdomain filtering method has the best effect on enhancing the edge location of the geological body, and the filtering window chooses 5×5 points and 1 iteration. The symmetrical shape small subdomain filtering method can preserve the main gradient direction when filtering the random noise, and the filtering window chooses 5×5 points and 2 iterations.
-
新生代印度-欧亚大陆的碰撞汇聚作用，导致晚始新世—中新世喜马拉雅地壳物质发生大规模的熔融作用，并形成巨量的淡色花岗岩、花岗岩体[1]。这些岩体沿碰撞造山作用有关的伸展构造分布[2-11]，并与造山带演化过程密切相关，是探索构造-岩浆相互作用关系和机制的重要“岩石探针”[12]。
在喜马拉雅山脉的脊部，自西向东出露一系列不连续的淡色花岗岩侵入体，岩体的分布受区域性伸展构造藏南拆离系（STDS）控制，并构成STDS韧性剪切带的一部分[13-14]，因此，又称其为STDS淡色花岗岩。前人对造山带西部的淡色花岗岩进行了较多的年代学和成因机制方面的研究，取得的主要认识和存在的争议包括：①淡色花岗岩的时代，基本厘定STDS淡色花岗岩的年代格架，时代集中于36~13Ma[2, 4-5, 7, 9, 15-22]，暗示长时间的地壳深熔作用；但多数年龄来自喜马拉雅中西部地区的岩体，缺少东部岩体年龄的补充和限定。②淡色花岗岩的源区和源岩，大多数学者认为，STDS淡色花岗岩源区和源岩为大喜马拉雅变质沉积岩[23-26]；而杨晓松等[27]根据部分熔融实验研究，认为源岩为大喜马拉雅中的黑云斜长片麻岩；Guo等[28]通过系统的地球化学和实验岩石学模拟研究，认为源区为大喜马拉雅变质沉积岩和来自小喜马拉雅流体的两端元混合区。③淡色花岗岩的形成机制，淡色花岗岩源自白云母脱水熔融这一结论已基本达成共识[29-30]，但对触发部分熔融的因素仍有争议；Le Fort等[31]认为，流体（如水）是导致部分熔融的因素，然而，实验岩石学研究表明，在无流体参与时，部分熔融行为仍可进行[24]；Har⁃ris等[25]、Davidson等[32]认为，构造减压是部分熔融的直接触发因素。后期，又有学者提出断裂的剪切摩擦生热和放射性同位素生热可以导致部分熔融[17, 30-33]，但理论计算表明，单纯的剪切摩擦热或地壳放射性元素生热难以触发大规模的岩浆作用[34]。由此可见，STDS淡色花岗岩的研究仍存在不足，已取得的部分认识多基于对造山带中西部岩体的分析，而缺少来自东部岩体的年代学、地球化学等的约束。
本次研究采集喜马拉雅造山带东部错那地区淡色花岗岩样品，进行锆石U-Pb年龄、全岩地球化学和Sr-Nd同位素研究，为STDS淡色花岗岩研究提供了新的数据，以期厘定深熔作用时代，查明其源区和源岩，并探讨成因机制和构造意义。
1. 区域地质背景
狭义的喜马拉雅造山带指雅鲁藏布江缝合带与主前锋逆冲断裂（MFT）之间的、由新生代印度-欧亚大陆碰撞形成的强烈变形、变质带。造山带自北向南发育的构造-岩石单元依次为特提斯喜马拉雅沉积系（THS）、STDS、大喜马拉雅结晶杂岩（GHC）、主中央逆冲断裂（MCT）、小喜马拉雅沉积岩系（LHS）、主边界逆冲断裂（MBT）、西瓦里克前陆盆地沉积（SS）和MFT（图 1）[35]。THS位于最北部，南北分别以STDS和雅鲁藏布江缝合带为界，主要发育一套古生代—始新世、经历极低级变质的硅质碎屑岩和碳酸盐岩沉积，在中部自西向东分布一系列片麻岩穹窿——北喜马拉雅片麻岩穹窿（NH⁃GD），是北喜马拉雅伸展构造的重要组成部分[36]。GHC位于MCT和STDS之间，由一套源岩时代为古元古代—奥陶纪的中高级变质结晶杂岩系组成[37]，上部紧邻STDS处出露大量淡色花岗岩。LHS位于MBT和MCT之间，主要由碎屑沉积岩和低级变质岩组成。最南部为SS，由古近纪—中新世的海相和陆相沉积组成。
![]()
图 1 喜马拉雅造山带中东段地质简图（据参考文献[35]修改）GCT—大反冲断层；STDS—藏南拆离系；MCT—主中央逆冲断裂；MBT—主边界逆冲断裂；MFT—主前锋逆冲断裂；NHGD—北喜马拉雅片麻岩穹窿Figure 1. Sketch geological map of the middle and east Himalayan orogen研究区位于喜马拉雅造山带东部的错那地区（图 1）。研究区以STDS为界，北部为THS，由三叠纪—白垩纪沉积岩组成，主要岩性为砂岩、粉砂岩、泥页岩，以及低级变质的板岩和千枚岩（图 2）。南部为GHC中高级变质的肉切村群和聂拉木群，其中，肉切村群主要岩石包括钙质石英片岩、泥质片岩、大理岩等（图 2），内部发育少量切层的淡色花岗岩脉（图 3-a、3-b）。聂拉木群主要由黑云斜长片麻岩、透辉大理岩、含石榴子石矽线黑云斜长变粒岩、含矽线黑云石英片岩、石英岩等组成。在STDS附近出露淡色花岗岩体（图 2、图 3-c），岩石整体呈灰白色，显微镜下观察主要矿物成分为石英（40%）、斜长石（25%）、钾长石（20%）、白云母（10%），以及少量的黑云母和石榴子石（图 3-d）。
![]()
图 3 藏南错那地区主要岩石野外及显微照片a—泥质片岩及侵入其中的淡色花岗岩脉；b—大理岩及侵入其中的淡色花岗岩脉；c—淡色花岗岩野外露头；d—淡色花岗岩矿物组合，正交偏光。Qtz—石英；Pl—斜长石；Kfs—钾长石；Mus—白云母；Bi—黑云母；Grt—石榴子石Figure 3. Photographs and photomicrograph of the rocks in Cuona area, southern Tibet2. 实验方法
2.1 LA-MC-ICP-MS锆石U-Pb定年
采用常规方法分选出锆石单矿物并制作成环氧树脂靶，打磨靶面使锆石中心暴露，拍摄锆石的可见光和阴极发光显微照片，然后进行LA-MCICP-MS锆石U-Pb同位素测试。锆石阴极发光图像（CL）在北京大学造山带与地壳演化教育部重点实验室的扫描电镜上完成。LA-MC-ICP-MS锆石U-Pb同位素测试在天津地质矿产研究所同位素实验室进行。所用仪器为Thermo Fisher公司生产的Neptune型多接受器等离子体质谱仪（MC-ICPMS），并结合美国ESL公司生产的UP193-FX ArF准分子激光器，激光剥蚀束斑直径为35μm，激光能量密度为2.5J/cm2，频率为10Hz，激光剥蚀物质以氦为载气送入Neptune的电感耦合等离子体，GJ-1作为外部锆石年龄标准进行U-Pb同位素分馏校正。在测试过程中，每测定7个样品点之前、之后各测定2次锆石标样GJ-1。数据的离线处理采用ICPMS⁃DataCal程序[38]完成，锆石U-Pb谐和图用Isoplot/ Ex3.00程序绘制[39]。
2.2 主量和微量元素
野外采集的样品挑选较新鲜、无蚀变的部分碎至200目，在加拿大ACME实验室完成全岩主量和微量元素测定。主量元素采用电感耦合等离子体发射光谱仪（ICP-ES）测定，测试精度优于5%；微量元素采用Agilent 7500a型电感耦合等离子体质谱仪（ICP-MS）测定，测试精度优于5%~10%。
2.3 Sr-Nd同位素
Rb-Sr、Sm-Nd同位素测试样品的化学分离在北京大学造山带与地壳演化教育部重点实验室完成，Sr、Nd同位素比值测定在天津地质矿产研究所热电离质谱仪TRITON上完成。同位素比值测试过程中，扇形磁分析器的有效半径为81cm，加速电压为10kV时分析质量数范围为3~320amu，分辨率大于450（10%峰谷定义），灵敏度大于3ion/100μmol，丰度灵敏度≤10×10-9。在样品测试过程中，测定的JN-DI Nd同位素标样和NBS-987 Sr同位素标样的同位素值分别为143Nd/144Nd=0.512107±0.000004 （2σ）和87Sr/86Sr=0.710265±0.000004 （2σ）。
3. 测试结果
3.1 LA-MC-ICP-MS锆石U-Pb年龄
错那淡色花岗岩样品CN-06的锆石阴极发光图像显示，锆石颗粒自形，长柱状，棱角清晰，颗粒长度为200~400μm，长宽比介于2∶1~4∶1之间。大多数锆石发育核-幔-边结构，核部和幔部为继承锆石，核部色调较暗，幔部形状不规则，色调灰白；边部宽度为40~60μm，色调黑暗并显示密集的韵律生长环带，表明为岩浆成因。少数锆石仅发育核-边结构，核部色调较亮，围绕核部发育具有韵律环带的暗色生长边（图 4-a）。
![]()
图 4 藏南错那淡色花岗岩代表性锆石CL图像(a)和U-Pb谐和图(b)Figure 4. CL images of representative zircons (a) and zircon U-Pb concordia diagram (b) of Cuona leucogranite, southern Tibet对样品CN-06共选取23个点位进行U-Pb同位素测试，全部分析点均位于锆石边部的暗色韵律环带上，测试数据列于表 1。23个分析点的Th和U含量分别为187×10-6~509×10-6和8186×10-6~8970×10-6，Th/U值极低，为0.02~0.06，均不大于0.10，与深熔成因的锆石特征吻合[40]。23个分析点的206Pb/ 238U年龄介于18.9~16.4Ma之间，给出的年龄加权平均值为17.7±0.3Ma（MSWD=2.0）（图 4-b），代表错那淡色花岗岩的最终结晶年龄。
表 1 藏南错那淡色花岗岩样品（CN-06）LA-MC-ICP-MS锆石U-Th-Pb同位素测试数据Table 1. LA-MC-ICP-MS zircon U-Th-Pb isotopic data of the leucogranite sample CN-06 from Cuona, southern Tibet分析点号 Pb Th U Th/U 同位素比值 表观年龄/Ma 10-6 207Pb/235U ±1σ 206Pb/238U ±1σ 208Pb/232Th ±1σ 207Pb/235U ±1σ 206Pb/238U ±1σ 208Pb/232Th ±1σ 01 2598 339 8357 0.04 0.0170 0.0009 0.00255 0.00013 0.00108 0.00006 17.1 0.9 16.4 0.8 21.9 1.3 02 2687 208 8481 0.02 0.0177 0.0009 0.00275 0.00014 0.00087 0.00010 17.8 0.9 17.7 0.9 17.6 2.0 03 2716 213 8475 0.03 0.0181 0.0010 0.00279 0.00014 0.00086 0.00009 18.2 0.9 17.9 0.9 17.5 1.9 04 2889 226 8970 0.03 0.0175 0.0009 0.00267 0.00013 0.00085 0.00007 17.7 0.9 17.2 0.9 17.3 1.3 05 2612 275 8418 0.03 0.0174 0.0009 0.00270 0.00014 0.00086 0.00007 17.5 0.9 17.4 0.9 17.5 1.4 06 2944 309 8375 0.04 0.0201 0.0010 0.00284 0.00014 0.00177 0.00011 20.2 1.0 18.3 0.9 35.7 2.3 07 2826 187 8499 0.02 0.0185 0.0011 0.00287 0.00015 0.00166 0.00027 18.6 1.1 18.5 1.0 33.6 5.4 08 2630 298 8395 0.04 0.0177 0.0009 0.00275 0.00014 0.00083 0.00006 17.9 0.9 17.7 0.9 16.7 1.1 09 2634 350 8344 0.04 0.0180 0.0009 0.00267 0.00014 0.00103 0.00007 18.1 0.9 17.2 0.9 20.9 1.3 10 2747 207 8481 0.02 0.0182 0.0010 0.00283 0.00014 0.00111 0.00010 18.3 1.0 18.2 0.9 22.4 2.0 11 3001 250 8433 0.03 0.0202 0.0011 0.00277 0.00014 0.00284 0.00022 20.3 1.1 17.9 0.9 57.3 4.5 12 2707 439 8254 0.05 0.0181 0.0010 0.00264 0.00014 0.00136 0.00010 18.2 1.0 17.0 0.9 27.4 2.1 13 2834 371 8321 0.04 0.0195 0.0010 0.00257 0.00013 0.00232 0.00015 19.6 1.0 16.5 0.9 46.8 2.9 14 2721 354 8336 0.04 0.0181 0.0009 0.00274 0.00014 0.00102 0.00007 18.2 0.9 17.6 0.9 20.7 1.4 15 2719 193 8496 0.02 0.0182 0.0009 0.00277 0.00014 0.00117 0.00012 18.3 0.9 17.8 0.9 23.7 2.4 16 3128 225 8456 0.03 0.0192 0.0019 0.00275 0.00014 0.00430 0.00048 19.3 1.9 17.7 0.9 86.7 9.6 17 2612 192 8501 0.02 0.0174 0.0010 0.00271 0.00014 0.00090 0.00025 17.5 1.0 17.4 0.9 18.3 5.0 18 2699 337 8354 0.04 0.0185 0.0010 0.00276 0.00014 0.00106 0.00008 18.6 1.0 17.8 0.9 21.4 1.6 19 2704 278 8412 0.03 0.0187 0.0010 0.00282 0.00014 0.00112 0.00009 18.8 1.0 18.1 0.9 22.7 1.7 20 2832 305 8379 0.04 0.0188 0.0010 0.00294 0.00015 0.00100 0.00009 18.9 1.0 18.9 1.0 20.2 1.8 21 3020 198 8483 0.02 0.0181 0.0010 0.00292 0.00015 0.00300 0.00024 18.2 1.0 18.8 1.0 60.6 4.8 22 2760 400 8290 0.05 0.0178 0.0009 0.00279 0.00014 0.00126 0.00011 18.0 0.9 18.0 0.9 25.4 2.3 23 2677 509 8186 0.06 0.0182 0.0011 0.00265 0.00014 0.00158 0.00022 18.3 1.1 17.1 0.9 32.0 4.4 3.2 主量和微量元素
样品全岩主量和微量元素数据见表 2。错那淡色花岗岩具有较高的SiO2（74.46%~75.57%）、Al2O3 （14.07%~14.64%）和K2O（4.19%~4.85%），以及较低的CaO（0.51%~0.67%）、MgO（0.07%~0.10%）和TiO2（0.02% ~0.04%）含量。K2O/Na2O值较高，介于1.09~1.31之间。在SiO2-K2O图解中，样品点位于高钾钙碱性系列（图 5-a）。A/CNK值为1.15~1.25，均大于1.1，在A/CNK-A/NK图解中，样品点全部落入过铝质区域（图 5-b）中。总体来看，错那淡色花岗岩属于高钾钙碱性过铝质S型花岗岩。
表 2 藏南错那淡色花岗岩全岩主量、微量和稀土元素测试数据Table 2. Whole-rock major, trace and rare earth element analyses of the Cuona leucogranites in southern Tibet样品号 CN-01 CN-02 CN-03 CN-04 CN-05 SiO2 74.59 74.58 74.46 74.87 75.57 Al2O3 14.48 14.64 14.60 14.31 14.07 Fe2O3 0.70 0.83 0.90 0.77 0.69 MgO 0.08 0.09 0.10 0.07 0.08 CaO 0.51 0.53 0.59 0.64 0.67 Na2O 3.69 3.64 3.51 3.90 3.84 K2O 4.85 4.68 4.41 4.48 4.19 TiO2 0.02 0.03 0.04 0.03 0.03 P2O5 0.20 0.19 0.18 0.18 0.15 MnOz 0.04 0.04 0.06 0.05 0.04 烧失量 0.8 0.7 1.1 0.7 0.7 合计 99.16 99.25 98.85 99.30 99.33 A/NK 1.28 1.32 1.38 1.27 1.30 A/CNK 1.18 1.22 1.25 1.15 1.16 Na2O+K2O 8.54 8.32 7.92 8.38 8.03 K2O/Na2O 1.31 1.29 1.26 1.15 1.09 CaO/Na2O 0.14 0.15 0.17 0.16 0.17 Be 7 12 10 16 14 Sc 3 3 4 3 3 Co 0.4 0.2 0.4 0.5 1.1 Ni 0.1 0.1 0.3 0.2 0.2 Zn 5 9 13 11 9 Ga 16.9 18.5 20.8 18.9 19.5 Rb 377.6 413.5 381.7 353.8 353.0 Sr 12.8 23.3 18.2 15.3 15.2 Y 9.7 12.8 20.1 13.3 11.2 Zr 32.6 39.9 52.9 41.0 25.2 Nb 14.2 17.7 20.5 18.3 18.0 Cs 25.3 30.9 31.2 25.7 23.7 Ba 17 47 37 24 23 La 3.4 5.5 6.8 5.2 4.0 Ce 7.2 11.2 14.9 10.0 8.9 Pr 0.86 1.37 1.79 1.17 0.99 Nd 3.6 3.5 6.4 3.3 3.3 Sm 1.03 1.43 2.15 1.55 1.28 Eu 0.09 0.14 0.14 0.13 0.10 Gd 1.29 1.83 2.57 1.49 1.42 Tb 0.27 0.39 0.56 0.35 0.32 Dy 1.39 2.07 3.24 1.95 1.93 Ho 0.26 0.36 0.63 0.41 0.36 Er 0.78 1.11 1.78 1.09 0.92 Tm 0.12 0.19 0.30 0.20 0.15 Yb 0.86 1.28 1.94 1.18 1.11 Lu 0.13 0.19 0.28 0.19 0.16 Hf 1.0 1.3 2.3 1.7 1.3 Ta 4.0 4.4 5.4 6.1 4.8 Pb 6.3 6.7 7.2 6.3 6.7 Th 2.4 3.4 5.0 2.9 2.4 U 5.4 5.7 8.1 4.1 4.2 Sn 14 17 20 14 12 Rb/Sr 29.50 17.75 20.97 23.12 23.22 LREE 16.18 23.14 32.18 21.35 18.57 HREE 5.10 7.42 11.30 6.86 6.37 TREE 21.28 30.56 43.48 28.21 24.94 (La/Sm)w 2.13 2.49 2.04 2.17 2.02 (Gd/Yb)N 1.24 1.18 1.10 1.04 1.06 (La/Yb)N 2.84 3.08 2.52 3.16 2.59 δEu 0.24 0.26 0.18 0.26 0.23 注：主量元素含量单位为%，微量和稀土元素含量单位为10-6；A/NK=摩尔Al2O3/ （Na2O+K2O），A/ CNK=摩尔Al2O3/ （CaO+Na2O+K2O）；δEu=2EuN/ （SmN+GdN），其中N为球粒陨石标准化值[41] ![]()
图 5 藏南错那淡色花岗岩的SiO2-K2O (a)、A/CNK-A/NK分类图解(b)及原始地幔标准化蛛网图(c)和球粒陨石标准化稀土元素配分模式图(d) (原始地幔和球粒陨石数值据参考文献[41])Figure 5. TAS diagram (a), A/CNK-A/NK diagram (b), primitive mantle-normalized trace element spider diagram (c) and chondrite-normalized REE patterns (d) of Cuona leucogranites, southern Tibet在原始地幔标准化微量元素蛛网图（图 5-c）中，显示Rb、Th、U和Pb的正异常，以及Ba、Nb、Sr和Zr的负异常。稀土元素总量较低，为21.28×10-6~ 43.48×10-6，轻、重稀土元素分馏较弱，（La/Yb）N为2.59~3.16。在球粒陨石标准化的稀土元素配分模式图（图 5-d）中，具有强烈的负Eu异常，δEu值为0.18~0.26。
3.3 Rb-Sr、Sm-Nd同位素
样品全岩Rb-Sr、Sm-Nd同位素测试结果列于表 3。经过t=18Ma校正后，4个样品的Isr值较高，为0.78982~0.79276，εNd(t)值较低，为-19.5~-18.2，SrNd同位素的初始比值变化不大。Nd同位素二阶段亏损地幔模式年龄为2469~2372Ma，表明淡色花岗岩可能源自古老地壳的重熔。
表 3 藏南错那淡色花岗岩全岩Rb-Sr和Sm-Nd同位素测试数据Table 3. Whole-rock Rb-Sr and Sm-Nd isotopic compositions of the Cuona leucogranites in southern Tibet样品号 t/Ma Rb/10-6 Sr/10-6 87Rb/86Sr 87Sr/86Sr 2σ Isr Sm/10-6 Nd/10-6 147Sm/144Nd 143Nd/144Nd 2σm TDM2/Ma εNd(t) CN-01 18.0 377.6 12.8 83.3240 0.81406 0.00020 0.79276 1.03 3.60 0.181524 0.511636 0.000004 2372 -19.5 CN-02 18.0 413.5 23.3 50.1266 0.81431 0.00006 0.80150 1.43 3.50 0.259220 0.511652 0.000004 2469 -19.4 CN-04 18.0 353.8 15.3 65.3153 0.80910 0.00001 0.79241 1.55 3.30 0.298001 0.511654 0.000004 2452 -19.4 CN-05 18.0 353.0 15.2 65.5963 0.80659 0.00001 0.78982 1.28 3.30 0.246091 0.511709 0.000006 2390 -18.2 注：87Rb/86Sr和147Sm/144Nd通过ICP-MS测试的微量元素Rb、Sr、Sm和Nd计算所得，计算公式为87Rb/86Sr=Rb/Sr×2.981，147Sm/144Nd=Sm/Nd×[0.531497+0.142521×（143Nd/144Nd）s]。ISr=（87Sr/86Sr）s+87Rb/86Sr（eλt–1），143Nd/144Nd（t）=（143Nd/144Nd)s+147Sm/144Nd （eλt-1）；εNd（t）=[ （143Nd/144Nd）s/ （143Nd/144Nd)CHUR-1]×104。（143Nd/144Nd）CHUR=0.512638，（147Sm/144Nd）CHUR=0.1967，（143Nd/144Nd)DM=0.51315，（147Sm/144Nd)DM=0.2137；λRb=1.42×10-12/年[42]，λSm=6.54×10-12/年[43]；二阶段模式年龄TDM2的计算见参考文献[41] 4. 讨论
4.1 STDS淡色花岗岩的时代
喜马拉雅淡色花岗岩的形成时代一直是地学界研究的重点，准确获得各岩体的年龄，不仅有助于限定相关构造的活动时代，而且有利于重建造山过程中岩浆作用的时空格局，进而反演造山带的岩浆演化历史。近期，吴福元等[1]根据目前已发表的年龄数据，重新勾画了喜马拉雅淡色花岗岩的年代格架，大致划分为始喜马拉雅阶段（Eo-Himalayan：44~26Ma）、新喜马拉雅阶段（Neo-Himalayan：26~ 13Ma）和后喜马拉雅阶段（Post-Himalayan：13~ 7Ma）。以上3个时期的花岗岩最显著的区别在于受造山带内不同构造的控制：始喜马拉雅阶段的花岗岩多分布于THS中部的NHGD核部，新喜马拉雅阶段的淡色花岗岩主要沿STDS分布，而后喜马拉雅阶段的花岗岩受控于南北向裂谷（NSTR）。错那淡色花岗岩锆石U-Pb年龄为17.7±0.3Ma （MSWD=2.0），且岩体沿STDS分布，应属于上述划分方案中的新喜马拉雅阶段淡色花岗岩。在造山带的中西部，相同峰期时代的STDS淡色花岗岩也广泛分布（图 1）。①在西部，玛纳斯鲁淡色花岗岩年龄为19.3~18.1Ma[30, 33, 45]。②在中部，吉隆STDS淡色花岗岩脉年龄为18.6~16.6Ma[14]，聂拉木淡色花岗岩年龄为17.1~16.8Ma[9]，希夏邦马淡色花岗岩年龄为17.3Ma[46]，达扎卡淡色花岗岩年龄为16.7Ma[4]，绒布淡色花岗岩年龄为16.8~16.2Ma[18, 47]，珠峰淡色花岗岩年龄为16.9~16.4Ma[2, 30]，定结淡色花岗岩年龄为16.9~15.1Ma[7, 48-49]。③在东部，东不丹地区淡色花岗岩年龄为17.6~16.8Ma[50]。这些淡色花岗岩均沿区域性伸展构造STDS分布，且具有近乎一致的峰期结晶年龄，暗示其形成可能与STDS的活动密切相关。
4.2 错那淡色花岗岩的源区和源岩
错那淡色花岗岩具有高SiO2和Na2O含量，铝过饱和，K2O＞Na2O，为高钾钙碱性过铝质花岗岩。全岩的Isr值较高（0.78982~0.79276），εNd(t)值较低（-19.5~-18.2），Rb/Sr值（17.75~29.50）远高于大陆地壳的平均值（0.32），与典型壳源成因花岗岩特征吻合。而淡色花岗岩体内部及附近均未发现幔源包体和相关的岩石，因此，笔者认为错那淡色花岗岩应为壳源成因。
目前，关于STDS淡色花岗岩的源岩问题仍存在分歧。杨晓松等[51]的高温高压岩石学模拟研究发现，在温压分别为770~980℃和1.0~1.4GPa，且无外来流体的条件下，GHC中黑云斜长片麻岩发生部分熔融，形成的熔体与STDS淡色花岗岩成分类似，且熔融后的残留物质与GHC下地壳麻粒岩成分相当，据此认为，淡色花岗岩源岩为GHC黑云斜长片麻岩[27]。Knesel等[52]通过实验岩石学研究发现，GHC变泥质岩早期发生白云母脱水熔融，形成的花岗质熔体成分与前人获得的STDS淡色花岗岩的地球化学特征一致，因此，认为源岩应该为GHC变质沉积岩。考虑到淡色花岗岩的壳源成因特性，将错那淡色花岗岩样品的Sr、Nd同位素成分与潜在的可能源岩，包括GHC变泥质岩、LHS变质沉积岩和北喜马拉雅片麻岩[53-56]，统一以t=18Ma重新进行计算并综合对比（图 6）。结果显示，4件样品均位于GHC变泥质岩区域。此外，错那淡色花岗岩的二阶段亏损地幔模式年龄TDM2值（2469~2372Ma）与GHC变泥质岩的相同，且与GHC变泥质岩中碎屑锆石获得的约2500Ma的峰值年龄重叠[57]，暗示其可能来自GHC变泥质岩的部分熔融。综上所述，错那淡色花岗岩的源岩为GHC变泥质岩。
![]()
图 6 藏南错那淡色花岗岩的Isr-εNd(t)图Figure 6. Isr-εNd(t diagram of the Cuona leucogranites, southern Tibet4.3 错那淡色花岗岩的形成机制
STDS淡色花岗岩是GHC变泥质岩部分熔融的产物，但部分熔融的机制仍有争议，争议点包括熔融的矿物是白云母还是黑云母？熔融过程是否有外来流体参与？
Knesel等[52]对GHC变泥质岩进行部分熔融实验研究发现，单一源区的变泥质岩在熔融的不同阶段具有不同的熔融机制，并产生化学成分迥异的熔体。早期在无水条件下，白云母发生脱水熔融，形成低Sr、高Isr的花岗质熔体。而后期伴随早期脱水熔融产生的含水流体的加入，主要发生水致部分熔融，形成高Sr、低Isr的熔体，分布于不同地区的花岗岩的Sr和Isr数据也证实了上述实验结论。本文搜集已发表的喜马拉雅花岗岩的Sr和Isr数据，统一以t=18Ma重新计算并归纳于图 7-a中。从该图可以看出，花岗岩主要集中分布于a1和a2两个截然不同的区域，其中a1区域内花岗岩表现出低Sr（10×10-6~ 80×10-6）、高Isr（0.75~0.80）的特征，与部分熔融实验中无水条件下白云母脱水熔融形成的熔体化学特征吻合。而a2区域内花岗岩显示高Sr（50×10-6~ 200×10-6）、低Isr（0.72~0.75）的特征，与水致部分熔融形成的熔体化学特征相符（图 7-a）。错那淡色花岗岩具有低Sr（12.8×10-6~23.3×10-6）和高Isr （0.78982~0.79276）的特征，在图 7-a中，与其他地区的STDS淡色花岗岩一样位于a1区域中，暗示其为无水条件下白云母脱水熔融的产物。
![]()
图 7 藏南错那淡色花岗岩的Sr-Isr (a)和Ba-Rb/Sr(b)图Mus(VP)—饱和蒸汽相的白云母熔融蚀变；Mus(VA)—缺乏蒸汽相的白云母熔融蚀变；Bi(VA)—缺乏蒸汽相的黑云母熔融蚀变Figure 7. Sr-Isr(a) and Ba-Rb/Sr(b) diagrams of the Cuona leucogranites, southern Tibet实验岩石学[52, 60]和理论计算结果[61-62]表明，随温压条件和含水量的变化，变泥质岩发生从白云母到黑云母的递进部分熔融，形成具有复杂RbSr组分的花岗质熔体[24, 52, 61]。在变泥质岩中，Rb主要赋存于云母类矿物中，而Sr主要赋存于长石矿物中，在有水参与的地壳熔融过程中，白云母和斜长石都发生部分熔融，且长石熔融的量远大于云母，导致熔体的Rb/Sr值偏低。在无水条件下，白云母大量熔融，形成的熔体Rb/Sr值较高[63]。因此，可以通过Rb/Sr值及其与Ba浓度的关系反推部分熔融机制。错那淡色花岗岩具有较高的Rb/Sr值（17.75~29.50），在Ba-Rb/Sr图解（图 7-b）中，Rb/ Sr值随着Ba浓度的增加而减小，二者呈负相关，与缺乏蒸汽相的白云母脱水熔融的趋势一致[63]，也与上述Sr-Isr系统分析结果一致。同时，由于无水条件下参与熔融的矿物主要为白云母，斜长石未进入熔体而残留在源区，导致错那淡色花岗岩具有较强的负Eu异常。综上所述，错那淡色花岗岩是GHC变泥质岩在无水条件下白云母脱水熔融的产物。
4.4 构造动力学意义
喜马拉雅造山带内发育多种类型的花岗岩，其形成过程对应不同的造山阶段，反映不同的构造背景和动力学过程。如晚始新世—渐新世（45~28Ma）的花岗岩主要出露于NHGD，如然巴穹窿（44Ma）[64]、麻布迦穹窿（35Ma）[65]、萨迦穹窿（28.1Ma）[66]、佩枯错穹窿（28.2Ma）[67]等。这些花岗岩的形成年龄与喜马拉雅变质岩的进变质时代相同[68-69]，是早期地壳逆冲增厚条件下部分熔融的产物[37]，熔体的产生使地壳弱化，并直接影响与STDS相关的中下地壳的伸展拆离，是触发STDS启动的关键因素。另一类为晚渐新世—早中新世（26~13Ma）[1]的淡色花岗岩，以STDS淡色花岗岩为代表，其地球化学成分与始新世—渐新世花岗岩不同，以高Rb/Sr值和高Isr为特征，结晶时代与喜马拉雅峰期变质作用时代相同，其形成可能与地壳伸展背景下GHC的快速折返/隆升导致的减压熔融有关[25]。最后一类为晚中新世（13~7Ma）花岗岩，这类花岗岩的矿物和化学成分与STDS淡色花岗岩类似，主要的区别在于前者主要沿NSTR分布，其形成与造山后期的东西向伸展有关。
在GHC变泥质岩部分熔融的过程中，温度和压力是主控因素。错那淡色花岗岩的锆石饱和温度为657~715℃，这一温度低于白云母脱水熔融的温度（720~770℃）。在缺少深部热源证据和剪切摩擦生热，地壳中放射性同位素衰变生热无法产生大规模岩浆作用的情况下[34]，减压成为部分熔融的必要条件，而STDS的启动不仅可以起到减压的作用，还可为岩浆上侵提供空间和通道，可能是触发熔融的重要因素。已有研究表明，STDS强烈活动的时间集中在28~13Ma[14]，与STDS淡色花岗岩的时代基本一致，且淡色花岗岩主要沿STDS分布，佐证了以上结论。错那淡色花岗岩的年龄（17.7±0.3Ma）与STDS活动时代吻合，其地球化学特征也与晚渐新世—早中新世淡色花岗岩相同。从构造-岩浆作用角度综合考虑，笔者认为，错那淡色花岗岩可能形成于STDS启动之后的伸展减薄背景下，反映了构造减压导致的中下地壳变泥质岩中含水矿物（白云母）脱水熔融的动力学过程。
5. 结论
（1）藏南错那淡色花岗岩的结晶年龄为17.7±0.3Ma，岩石属于高钾钙碱性过铝质花岗岩，富集Rb、Th和U，亏损Ba、Nb、Sr和Zr，具有强烈的负Eu异常；Rb/Sr值较高，Isr和εNd(t)值分别为0.78982~0.79276和-19.5~-18.2。
（2）错那淡色花岗岩为壳源成因，源岩是GHC变泥质岩，是无水条件下白云母脱水熔融的产物。熔体的形成与STDS的启动有关，反映了地壳伸展减薄背景下，构造减压导致中下地壳变泥质岩中含水矿物脱水熔融并沿STDS上升侵位的动力学过程。
致谢: 成文过程中得到中国地质调查局西安地质调查中心卢进才教授级高级工程师的指导，以及西安思坦仪器股份有限公司杨忠华工程师的帮助，在此表示衷心的感谢。 -
![]()
图 1 模型AB及其理论重力异常（参考文献[9]修改）
a—模型A及其分布；b—模型A理论重力异常；c—模型B及其分布；d—模型B理论重力异常
Figure 1. The model of AB and its theoretical gravity anomalies
![]()
图 2 模型A重力异常的小子域滤波处理结果
Figure 2. Results of model A by use of small subdomain filtering
![]()
图 3 模型B重力异常的小子域滤波处理结果
Figure 3. Results of model B by use of small subdomain filtering
![]()
图 4 不同滤波窗口、迭代次数对称形子域滤波法对模型B重力异常的处理结果
a—9×9数据点滤波窗口、迭代1次；b—13×13数据点滤波窗口、迭代1次；c—9×9数据点滤波窗口、迭代2次；d—13×13数据点滤波窗口、迭代2次
Figure 4. Results of model B by use of symmetrical shape small subdomain filtering
![]()
![]()
图 6 不同滤波窗口、迭代次数对称形子域滤波法对模型B含噪重力异常的处理结果
a—9×9数据点滤波窗口、迭代1次；b—13×13数据点滤波窗口、迭代1次；c—9×9数据点滤波窗口、迭代2次；d—13×13数据点滤波窗口、迭代2次
Figure 6. Results of noised model B by use of symmetrical shape small subdomain filtering with different filter windows and different iterations
![]()
图 7 不同滤波窗口3种类型小子域滤波法对模型B重力异常的处理结果
a—Tomita-Tsuji形子域9×9数据点滤波窗口、迭代1次；b—Tomita-Tsuji形子域13×13数据点滤波窗口、迭代1次；c— “田”字形子域9×9数据点滤波窗口、迭代1次；d— “田”字形子域13×13数据点滤波窗口、迭代1次；e—滑动形子域9×9数据点滤波窗口、迭代1次；
Figure 7. Results of model B by use of 3 different shapes of small subdomain filtering with different filter windows
![]()
图 8 不同迭代次数3种类型子域滤波法对模型B重力异常的处理结果
a—Tomita-Tsuji形子域5×5数据点滤波窗口、迭代2次；b—Tomita-Tsuji形子域5×5数据点滤波窗口、迭代3次；c— “田”字形子域5×5数据点滤波窗口、迭代2次；d— “田”字形子域5×5数据点滤波窗口、迭代3次；e—滑动形子域5×5数据点滤波窗口、迭代2次；f—滑动形子域5×5数据点滤波窗口、迭代3次
Figure 8. Results of model B by use of 3 different shapes of small subdomain filtering with different iterations
![]()
![]()
图 10 SH区块断裂体系及利用地震、电法剖面确定的断裂
a—断裂平面分布；b—A-A′地震剖面；c—B-B′地震剖面；d—C′-C电法视电阻率剖面
Figure 10. Fault structures of block SH and the evidence by seismic and electrical sections
表 1 理论模型参数
Table 1 Parameters of the theoretical model
模型编号 地质体编号 x，y角点坐标/m 顶面埋深/m 底面埋深/m 密度差/(g· cm-3) A1 （3000，9000）; （5000，9000）; （5000，11000）; （3000，11000） 100 2100 0.7 A A2 （9000，9000）; （11000，9000）; （11000，11000）; （9000，11000） 500 2500 0.7 A3 （15000，9000）; （17000，9000）; （17000，11000）; （15000，11000） 1100 3100 0.7 B1 （7125，7485）; （8050，7485）; （8050，8415）; （7125，8415） 100 200 0.7 B2 （3855，3640）; （1960，2640）; （5160，855）; （6250，1320） 400 1000 0.7 B B3 （1425，7292）; （1425，5278）; （3554，5278）; （2986，6920） 600 1000 0.7 B4 （0，8750）; （3100，8150）; （5800，4150）; （10000，2733）;
（10000，3200）; （6150，4500）; （3450，8500）; （0，9170）200 400 0.7 
下载: 导出CSV
-
杨高印.位场数据处理的一项新技术——小子域滤波法[J].石油地球物理勘探, 1995, 30(2):240-244. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=sydq199502011&dbname=CJFD&dbcode=CJFQ 王万银, 邱之云, 杨永, 等.位场边缘识别方法研究进展[J].地球物理学进展, 2010, 25(1):196-210. http://d.old.wanfangdata.com.cn/Periodical/dqwlxjz201001027 严良俊, 胡文宝, 姚长利.重磁资料面积处理中的滤波增强技术与应用[J].勘探地球物理进展, 2006, 29(2):102-103. http://d.old.wanfangdata.com.cn/Periodical/ktdqwljz200602006 马涛, 王铁成, 王雨.一种改进的网格数据保持梯度滤波方法[J].石油地球物理勘探, 2007, 42(2):198-201. http://kns.cnki.net/KCMS/detail/detail.aspx?filename=sydq200702014&dbname=CJFD&dbcode=CJFQ 张凤旭, 张凤琴, 刘财, 等.断裂构造精细解释技术——三方向小子域滤波[J].地球物理学报, 2007, 50(5):1543-1550. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dqwlxb200705031 张凤旭, 刘万崧, 张兴洲, 等.用重力三方向小子域滤波解释伊通盆地断裂[J].地球科学——中国地质大学学报, 2009, 34(4):665-672. http://d.old.wanfangdata.com.cn/Periodical/dqkx200904013 张凤旭, 张兴洲, 张凤琴, 等.中国东北地区重力场研究——利用改进的三方向小子域滤波划分主构造线及大地构造单元[J].地球物理学报, 2010, 53(6):1475-1485. http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201006025 肖锋, 吴燕冈, 孟令顺.位场数据处理中小子域滤波的改进[J].石油地球物理勘探, 2010, 45(1):136-139. http://d.old.wanfangdata.com.cn/Periodical/sydqwlkt201001023 蒋甫玉, 高丽坤.基于改进的小子域滤波的重力异常及重力梯度张量边缘增强[J].应用地球物理, 2012, 9(2):119-130. http://manu09.magtech.com.cn/jweb_yydqwl/CN/abstract/abstract55.shtml 马国庆, 杜晓娟, 李丽丽.优化小子域滤波方法研究及其应用[J].石油地球物理勘探, 2013, 48(4):658-662. http://d.old.wanfangdata.com.cn/Periodical/sydqwlkt201304022 段晓旭. 重磁异常小子域滤波算法的分析研究与改进[D]. 中国地质大学(北京)硕士学位论文, 2014. http://cdmd.cnki.com.cn/Article/CDMD-11415-1014233762.htm 蔡钟, 倪小东.重力异常处理中小子域滤波法的改进[J].河北工程大学学报(自然科学版), 2015, 32(4):47-51. http://d.wanfangdata.com.cn/Periodical_hbjzkjxyxb201504011.aspx 夏玲燕, 吴汉宁, 柏冠军, 等.柴达木盆地航磁资料微弱信息增强技术研究及在线性构造识别中的应用[J].地球物理学进展, 2008, 23(4):1058-1062. http://manu39.magtech.com.cn/Geoprog/CN/abstract/abstract2509.shtml 许海红, 袁炳强, 杨高印, 等.小子域滤波联合总水平导数在重力资料处理中的应用——以内蒙古银额盆地EQ区块断裂构造识别为例[J].地质通报, 2014, 33(11):1853-1860. doi: 10.3969/j.issn.1671-2552.2014.11.024 Wilhelm B, Mark J B. Principles of Digital Image Processing: Advanced Methods[M]. Springer-Verlag London, 2013.
曾华霖.重力场与重力勘探[M].北京:地质出版社, 2005:235-238.
计量
- 文章访问数: 2567
- HTML全文浏览量: 281
- PDF下载量: 2042
