郯庐断裂带是纵贯亚洲东部的一条巨型深大断裂带，自1957年航磁调查发现以来，至今已有60多年的研究历史。断裂带南起于长江，途径大别造山带东缘，经庐江、肥东地区，向北横穿苏鲁造山带，经郯城、渤海过沈阳后分为东西两支延伸至俄罗斯境内。郯庐断裂带开始活动的时间，被认为是中国东部中生代特提斯构造域向滨太平洋构造域发生转换的起点(董树文等, 2000, 2007；张岳桥等，2007；宋传中等, 2010a, b ，2011，2014；Zhu et al., 2015)。该断裂带的活动记录了中国东部中生代以来多期、复杂的构造活动与演变(万天丰，1995；万天丰等，1996；陈宣华等，2000；Zhu et al., 2009, 2010, 2015, 2018；Zhao et al., 2016；朱光等，2016；Gu et al., 2017)，一直深受国内外学者的关注。
郯庐断裂带最显著的特征就是走滑平移，而由于走滑作用产生了多处大型韧性剪切。其中，最典型的韧性剪切带分布在其中南段，即苏鲁造山带西缘、肥东地区(张八岭隆起带)、大别山东缘。剪切带内发育大量的糜棱岩，很好地记录了郯庐断裂带的变形特征。自二十世纪九十年代以来，对郯庐断裂带中南段的研究非常深入，由此也获得了一系列丰硕的成果(徐嘉炜等, 1984, 1995；万天丰等，1996；王小凤等，1998；陈宣华等，2000；王小凤，2000；朱光等, 2001, 2002, 2003, 2004, 2005a, b ，2006a, b ；牛漫兰等, 2002, 2005, 2008；汤加富和徐卫，2002；侯明金等, 2003, 2006；宋传中等，2003；汤加富等，2003；王勇生等, 2004, 2005；Lin et al., 2005, 2009；Zhu et al., 2005, 2009, 2010, 2018；Wang，2006；牛漫兰，2006；Zhang et al., 2007；张岳桥和董树文，2008；石永红等，2009；康涛，2014；赵田等, 2014a, b ，2016；Zhao et al., 2016；林伟等，2016)。
肥东地区处于郯庐断裂带的中南部，是张八岭隆起南段的重要组成部分，也是郯庐断裂带内韧性剪切带出露的主体区域之一(宋传中等，2003；朱光等，2005b；韩雨和牛漫兰，2016；王微等, 2016a, b )。肥东地区韧性剪切带总体呈NNE向展布(图 1)，长约40km，宽7~10km，区内发育大量NNE-SSW走向左行韧性剪切带，整体面理倾向SE，倾角较陡，线理倾伏向NE-SW，倾伏角多小于10°或在10°~20°之间。区内浮槎山、西韦、桃源一带韧性剪切带出露较好，研究成果也较多(刘德良等，1996；童劲松和许卫，2000；许卫等，2001；宋传中等，2003；Lin et al., 2005；朱光等, 2005b, 2006b；吴小奇等，2006；Zhang et al., 2007；马涛, 2014；韩雨和牛漫兰，2016；王微等, 2015, 2016a, b )。
图 1
Fig. 1
图 1
郯庐断裂带肥东地区地质简图(据王微等，2016a修改)
Fig. 1
The geological sketch map of Feidong area in Tan-Lu Fault Zone (modified after Wang et al., 2016a)
肥东地区是郯庐深断裂变形最强烈的地段，也是研究郯庐断裂带深部变形条件和动力学过程的良好场所(宋传中等，2003)，前人对该地区的深入研究已经取得了一系列成果。但由于受研究手段的制约，目前对肥东地区构造样式和构造变形分析方面的研究仍不全面、系统，特别是缺乏精细化构造变形研究。为此，笔者尝试从肥东北部文集地区入手，选择并精细测制了两条构造剖面，通过对剖面上构造变形的详细野外观测结合室内显微构造变形分析，来确定其涡度、有限应变、分维度、古应力差、应变速率等参数，从几何学、运动学及动力学等方面展开精细化构造变形研究，以期反应出郯庐断裂带肥东段的构造样式和运动学特征，探讨郯庐断裂带的起源及演化过程。
1 地质概况 位于郯庐断裂带中南部的肥东地区，是整个郯庐断裂带内韧性剪切带出露最好的地区之一。区内发育多条次级韧性剪切带，并在地表出露了一系列糜棱岩类岩石，很好地记录了郯庐断裂带的变形特征，是研究郯庐断裂中深层次构造变形的理想场所。
研究区位于肥东北部文集地区，龙山南麓，区内主要出露角闪岩相肥东变质杂岩。构造岩以强直片麻岩、糜棱岩、超糜棱岩及同构造伟晶岩脉为主，岩石变形强烈。面理向SE陡倾，矿物拉伸线理主要向NNE缓倾。笔者在对肥东北部文集地区两条构造剖面的野外调查中发现，该地区面理呈一背形构造，局部可见鞘褶皱发育。岩性方面呈现出一定的对称性，核部物质以深色的黑云角闪斜长质糜棱岩及超糜棱岩为主，夹同构造伟晶岩脉，向两翼物质以浅色的长英质糜棱岩、片麻岩为主。且核部岩石变质、变形强，两翼渐弱，近似对称性。本文在研究区南北两端选择两条近平行的构造剖面来探讨郯庐断裂带的构造特征，它们横穿了肥东杂岩，岩性出露较全，具备详细研究郯庐断裂带肥东段的构造变形的有利条件。
2 地质特征与几何学形态分析 本次横切肥东北部文集地区两条近平行的构造剖面(图 2a)，即文集-龚刘高剖面(A-B)和文集-郭杨剖面(C-D)，岩石类型复杂、变质变形强烈，为精细研究肥东北部郯庐断裂带的变质变形特征和形成环境提供了难得的第一手资料(图 3)。
图 2
Fig. 2
图 2
肥东北部文集地区地质简图和应力场分布图
Fig. 2
The geological sketch map and the stress field distribution diagrams of Wenji area in the northern Feidong
2.1 构造剖面及其地质特征 2.1.1 文集-龚刘高(A-B)剖面 沿剖面自NW向SE，出露的岩性主要为不同程度糜棱岩化的长英质片麻岩、黑云角闪斜长片麻岩、斜长角闪片麻岩，以及大量伟晶岩脉。面理以片麻理、糜棱面理为主。经过野外观测和室内研究发现，自剖面NW端至SE端，面理产状一致向SE倾，但倾角有规律变化，NW端(NO.000-NO.005)倾角较SE端(NO.007-NO.013)略陡，这一变化在图 3中的面理线理赤平投影(本文均采用吴氏网下半球投影)图中(图 3a, b)可以看出来。另外，剖面中还发育许多次级褶皱，笔者对NO.006点发育的次级褶皱两翼产状也进行了赤平投影(图 3c)：两翼产状分别为：118°∠55°、100°∠15°，赤平投影图中两翼的投影交点即为次级褶皱枢纽(B)的投影，产状为：30°∠6°；过枢纽法线的投影点做一大圆弧(π圆)，将两翼截π圆所在的弧一分为二得到轴迹，使枢纽和轴迹位于同一大圆弧即可得到轴面(A)产状：113°∠32°。可以发现，次级鞘褶皱(图 4a)的枢纽走向NNE-SSW，产状较缓，与区域性线理产状基本一致；轴面产状倾向SE，与区域性面理产状基本一致。
图 3
Fig. 3
图 3
肥东北部文集地区A-B及C-D构造剖面图和应变速率图
(a) A-B剖面NW部面理线理赤平投影；(b) NO.006褶皱两翼产状赤平投影；(c) A-B剖面SE部面理线理赤平投影；(d) C-D剖面NW部面理线理赤平投影；(e) NO.019褶皱两翼产状赤平投影；(f) C-D剖面SE部面理赤平投影
Fig. 3
The structural maps and strain rate diagrams of A-B and C-D sections of Wenji area in the northern Feidong
(a) northwestern foliation and lineation stereographic projection in A-B section; (b) NO. 006 two wings of the fold attitude stereographic projection; (c) southeastern foliation and lineation stereographic projection in A-B section; (d) northwestern foliation and lineation stereographic projection in C-D section; (e) NO. 019 two wings of the fold attitude stereographic projection; (f) southeastern foliation and lineation stereographic projection in C-D section
2.1.2 文集-郭杨(C-D)剖面 沿剖面自NW向SE，构造岩主要有强直片麻岩、糜棱岩及伟晶岩脉。面理以片麻理、糜棱面理为主。较北侧的文集-龚刘高剖面(A-B剖面)，本条剖面岩性出露不是很完全，但仍有规律可循。自剖面NW端至SE端，面理产状基本倾向SE，部分地区褶皱发育，面理产状有所变化，NW端(NO.014-NO.018)面理倾角较陡(图 3d)，而SE端(NO.020-NO.022)面理产状较缓(图 3f)。剖面中也发育许多次级褶皱，同样对NO.019点发育的次级褶皱两翼产状进行了赤平投影(图 3e)，两翼产状分别为：330°∠60°、130°∠50°，赤平投影图中两翼的投影交点即为次级褶皱枢纽(B)的投影，产状为：52°∠14°；过枢纽法线的投影点做一大圆弧(π圆)，将两翼截π圆所在的弧一分为二得到轴迹，使枢纽和轴迹位于同一大圆弧即可得到轴面(A)产状：141°∠84°。可以发现，次级鞘褶皱的枢纽走向NE-SW，产状较缓，与区域性线理产状基本一致；轴面产状倾向SE，与区域性面理产状基本一致。
2.2 几何学形态分析 2.2.1 几何学形态 野外观察和发现，肥东北部文集地区沿剖面方向上看形态为一大背形，其枢纽走向NE-SW，轴面倾向SE，与区域性线理、面理产状基本一致。枢纽与轴面同在C面理上，稳定于郯庐断裂带的总体产状。通过大量野外观测发现，郯庐断裂带肥东段的构造样式具有A型褶皱的特征，其褶轴与拉伸线理方向基本一致。成因可能是郯庐断裂带剪切活动在大别造山带背形构造的基础上变形加强的结果。从构造剖面图中可以看出，文集地区大背形的核部在龙山附近，转折端部位宽而缓，自NO.006至NO.011部分可能都位于转折端内部，以对称出现同构造伟晶岩脉为依据。面理产状局部西倾，内部发育许多背形褶皱(图 4b，d，f)、杆状构造(图 4c)及次级小褶皱(图 4e)。
图 4
Fig. 4
图 4
肥东北部文集地区野外构造变形特征
(a) NO.002糜棱岩中长石残斑；(b) NO.006背形褶皱；(c) NO.007杆状构造；(d) NO.008背形褶皱；(e) NO.019次级小褶皱；(f) NO.019背形褶皱；(g) NO.019片麻岩中的S-C组构；(h) NO.022伟晶岩中的长石残斑
Fig. 4
The field structural deformation characteristics in Wenji area in the northern Feidong
(a) NO.002 feldspar porphyroclast in the mylonite; (b) NO. 006 antiformal fold; (c) NO.007 rodding structure; (d) NO.008 antiformal fold; (e) NO.019 secondary fold; (f) NO.019 antiformal fold; (g) NO.019 S-C fabric in gneiss; (h) NO.022 feldspar porphyroclast in the pegmatite
2.2.2 物质分布规律 研究区内岩性出露的有黑云角闪片麻岩、黑云斜长片麻岩、花岗片麻岩等，构造岩类型以糜棱岩、超糜棱岩及大量同构造伟晶岩脉为主。由于构造活动影响，片麻岩几乎都受到了不同程度的糜棱岩化。分布规律呈现出一定的对称性，核部物质以中-基性糜棱岩为主，夹同构造伟晶岩脉，两翼物质渐为中-酸性，近似对称性分布。
3 构造变形的运动学分析 3.1 应力状态分析 由构造作用造成的应力场为构造应力场，它可以直接的反映地质体所受外力作用的状态，并帮助推导出构造应力的方向及来源。郯庐断裂带作为中国东部一条主要的构造带，在不同时期分别受到特提斯与环太平洋两大构造体制的作用力的作用(王微等，2016a)。通过对肥东北部的应力场分析，可以就郯庐断裂带的形成、演化时的构造应力有更加直观的了解。
目前，推求古应力场的方法主要应用于脆性断层，即利用断层面及其上的擦痕产状进行构造古应力场分析。而同一断裂体系中的脆、韧性断裂带是在相同动力机制下不同构造域形成的不同变形特征的反映，其中，韧性剪切带中的C面理往往是由片状矿物构成的流劈理。本文利用韧性剪切带的剪切C面理类比脆性断层的断层面、韧性剪切带的矿物线理类比脆性断层内的擦痕，借鉴脆性断层分析古应力的方法来半定量分析研究区的古应力场。具体来说，就是通过对文集地区两条剖面内的岩石面理、线理的产状及其组合关系的观测，运用Win-Tensor软件分析处理，计算出研究区的应力方向，并进行制图(图 2b)。区内岩石面理向SE陡倾，矿物拉伸线理多为SW向缓倾伏，产状十分稳定。用面理模拟断层面，矿物线理模拟擦痕，为保证数据的准确性，每个点的测量数据为8组，进行构造应力场的反演(表 1)。
表 1
Table 1
表 1(Table 1)
表 1
肥东北部文集地区应力场分析结果
Table 1
The analysis results of stress field of Wenji area in the northern Feidong
点号
位置
σ1(az/pl)
σ2(az/pl)
σ3(az/pl)
NO000
N31°57′31″、E117°43′20″
349°∠2°
82°∠47°
258°∠43°
NO003
N31°57′19″、E117°43′42″
8°∠8°
107°∠52°
272°∠37°
NO006
N31°57′23″、E117°43′53″
345°∠14°
97°∠56°
247°∠30°
NO014
N31°56′36″、E117°42′45″
350°∠10°
106°∠67°
257°∠21°
NO015
N31°56′35″、E117°42′44″
345°∠16°
103°∠58°
247°∠27°
NO016
N31°56′39″、E117°42′41″
345°∠16°
103°∠59°
247°∠26°
NO017
N31°56′38″、E117°42′41″
354°∠10°
100°∠58°
258°∠30°
NO021
N31°55′57″、E117°43′59″
17°∠125°
120°∠26°
250°∠52°
NO022
N31°55′58″、E117°43′57″
15°∠10°
117°∠50°
277°∠38°
注：σ1-最大主应力；σ2-中间主应力；σ3-最小主应力；az-倾伏向；pl-倾伏角
表 1 肥东北部文集地区应力场分析结果
Table 1 The analysis results of stress field of Wenji area in the northern Feidong
结果显示，郯庐断裂带肥东北部文集地区的主压应力方向(σ1)非常稳定，倾伏向在345°~354°之间，倾伏角近水平，表明构造体主要受到近南北方向的挤压，是由于该时期内中国东部构造主要受环太平洋构造控制。
3.2 剪切类型分析 野外可见肥东北部地区韧性剪切带内左旋剪切的运动学方向十分清晰，垂直于面理、平行于线理方向上(即XZ面)，岩石中发育着许多指示左旋剪切的小构造，如残斑(图 4a，h)，S-C组构(图 4g)等。从剖面的NW端到SE端都有这些运动学特征标志的保留。
剪切带内物质流变的运动学特征是阐明地壳构造演化的关键，也是论释变形岩石中组构运动学意义的重要方面(Xypolias，2009)。综合前人研究来看，郯庐断裂带不同区段的走滑运动是以简单剪切为主。为了进一步探查肥东北部剪切带的运动学特征，笔者结合有限应变分析对肥东北部文集地区韧性剪切带进行运动学涡度分析(Wk)，进而确定该剪切带的剪切作用类型。
一般剪切带中运动学涡度的确定具有重要的构造意义，运动学涡度的值实际上反映了剪切带中简单剪切和纯剪切的相对份量。它还可以表征剪切作用类型，进而解释区域构造带特别是韧性剪切带的运动学特征。
涡度的概念来源于流体力学，是指某一运动形式所拥有的旋转量(王勇生等，2004)。Truesdell(1954)首次提出了运动学涡度(Wk)的概念，Means et al.(1980)首次把Wk引入地质学，成功地实现了递进变形过程非共轴度的定量化。随后，Passchier(1987)等构造地质学家陆续对天然变形岩石进行了涡度分析，并提出了许多分析方法。运动学涡度分析已成为构造定量化必不可少而又非常有效的工具之一。
自Passchier第一次对天然变形糜棱岩进行涡度分析以来，许多构造地质学家相继提出了多种有效的分析方法。目前，估算运动学涡度的方法已有多种，如变形脉体法(Passchier，1990；Short and Johnson, 2006)，刚性颗粒或旋转碎斑法(Xypolias and Koukouvelas, 2001；Jessup et al., 2007)，Rs-θ法(张波等，2008)，有限应变轴率RXZ/石英C轴组构法(Xypolias and Koukouvelas, 2001；Xypolias，2009)，石英C轴组构法(郑亚东等，2008)，极摩尔圆法(张进江和郑亚东, 1995, 1997；王新社等，2005；刘江等，2012)等。
极摩尔圆法是目前涡度计算方法中定量化程度最高的方法(王勇生等，2004)。综合分析肥东北部剪切带内变形岩石宏观和显微等变形特征，笔者选用了极摩尔圆法对研究区内20个变形岩石样品进行了运动学涡度估算。
极摩尔圆的构建方法由Simpson and De Poar(1993)提出，张进江和郑亚东(1995)改进了这一方法，提出了三种构建极摩尔圆的方法。通过剪切带长度轴比Rs和最大拉伸方向与剪切方向的夹角θ两个参数建立相应的极摩尔圆来计算运动学涡度，Rs为XZ面的有限应变。
本文采用长短轴法，显微镜下对20个变形岩石样品中的变形长石进行测量统计，每个样品至少量30组。测出XZ面上的有限应变椭圆长短轴比Rs和椭圆长轴与剪切带方向的夹角θ，并取平均值，利用以上两个参数构建极摩尔圆，得到流脊间的夹角V，即可求得运动学涡度Wk=cosV(表 2、图 5)。
表 2
Table 2
表 2(Table 2)
表 2
肥东北部文集地区岩石涡度分析结果
Table 2
The vorticity analysis results of rocks in Wenji area in the northern Feidong
样品号
θ(°)
Rs
V(°)
Wk
WZ000
10.176
3.305
39.82
0.768
WZ002-1
9.899
3.444
40.1
0.765
WZ003-3
12.917
2.928
34.08
0.828
WZ003-6
15.439
2.732
29.56
0.870
WZ004-2
11.077
3.307
38.92
0.778
WZ008-1
13.554
3.324
32.45
0.844
WZ017-1
13.075
3.921
33.93
0.830
WZ018-2
13.210
3.280
33.79
0.831
WZ019-5
14.072
3.703
31.93
0.849
WZ020-2
13.134
3.281
33.87
0.830
WZ001-1
9.731
3.356
40.27
0.763
WZ002-2
8.841
3.231
43.16
0.729
WZ003-5
9.753
3.195
40.25
0.763
WZ003-8
10.857
3.342
38.14
0.787
WZ006-1
11.627
3.449
36.37
0.805
WZ008-3
8.209
3.947
42.79
0.734
WZ018-1
12.226
3.246
35.77
0.811
WZ019-1
11.627
3.953
37.37
0.795
WZ020-1
13.500
3.226
32.50
0.843
WZ021-1
11.547
3.339
36.45
0.804
表 2 肥东北部文集地区岩石涡度分析结果
Table 2 The vorticity analysis results of rocks in Wenji area in the northern Feidong
图 5
Fig. 5
图 5
肥东北部文集地区韧性剪切带运动学涡度极摩尔圆法图解
Rs-岩石有限应变；V-流脊间夹角；Wk-运动学涡度
Fig. 5
The kinematic vorticity polar Mohr diagram of the ductile shear zone in Wenji area in the northern Feidong
Rs-rock finite strain; V-angle between the ridges; Wk-kinematic vorticity
计算结果显示，郯庐断裂带肥东北部文集地区构造带中的岩石有限应变Rs在2.732~3.953范围内，涡度值Wk在0.729~0.870范围内，岩石变形强，涡度值很稳定，均大于0.71，小于1。结合应力场分析，说明肥东北部的韧性剪切带是以简单剪切为主，纯剪切为辅，具压扭性质的左行韧性剪切带。
4 显微变形特征与显微构造变形分析 4.1 显微变形特征 通过对肥东北部文集地区两条构造剖面进行观测发现，区内出露的变形岩石主要有糜棱岩和糜棱岩化的片麻岩，岩石变形强烈。对野外采集的定向标本进行室内磨片，切片方位沿XZ面，即沿垂直于面理、平行于线理方向切片，而后在显微镜下进行观察鉴定。
糜棱岩：主要以长英质为主，夹少量基性(角闪斜长质)糜棱岩。镜下观察：残斑主要为长石(Fl)、石英(Q)，基质主要为石英、长石以及少量的黑云母(Bt)；共生矿物组合主要为石英+长石+黑云母，具有典型的糜棱结构、流动构造；石英多发生边界迁移式(GBM)动态重结晶(图 6c-e)，颗粒边界呈锯齿状，部分为亚颗粒旋转式(SR)向颗粒边界迁移式(GBM)转变(图 6f)，粒径较小，0.02~0.3mm，波状消光，呈他形粒状、多晶条带状(图 6g)，变形强者呈拔丝状，颗粒长宽比很大，含量约50%~55%；长石残斑多呈眼球状、枣核状，偶见拖尾，指示运动学方向为左旋(图 6a, b)，波状消光，明显发生塑性变形，在残斑边缘位置可见新晶出现，长石发生膨凸式(BLG)动态重结晶(图 6h)，粒径大小0.3~2.5mm，含量约30%~35%；黑云母含量较少，多呈针状，以基质形式存在，围绕残斑沿糜棱面理分布，粒径较小，0.01~0.03mm，含量约8%~15%。
图 6
Fig. 6
图 6
肥东北部文集地区岩石显微构造图片
(a、b)糜棱岩的流动构造及长石残斑；(c-e)石英颗粒边界迁移式(GBM)动态重结晶；(f)石英亚颗粒旋转(SR)向边界迁移式(GBM)动态重结晶转变，可看到核幔构造，开始出现锯齿状边界；(g)石英的塑性流动及单晶条带；(h)长石的膨凸式(BLG)动态重结晶，出现新生亚颗粒.Q-石英；Fl-长石；Bt-黑云母
Fig. 6
The microscopic structure photos of rocks in Wenji area in the northern Feidong
(a, b) feldspar porphyroclast and flow structure of the mylonite; (c-e) grain boundary migration (GBM) dynamic recrystallization of quartz; (f) subgrain rotation (SR) transformed into grain boundry migration (GBM) dynamic recrystallization of quartz, core and mantle structures can be seen, beginning to appear serrated boundaries; (g) plastic flow and single crystal band of quartz; (h) bulging (BLG) dynamic recrystallization of feldspar, with new occurrence of subgranules. Q-quartz; Fl-feldspar; Bt-biotite
糜棱岩化片麻岩：区内主要出露庙山组(Mgn)片麻岩，以花岗质、长英质为主。镜下观察：岩石的矿物组成与糜棱岩基本一致，共生矿物组合主要为石英+长石+黑云母，具粒状变晶结构，片麻状构造；石英具有波状消光，发生塑性变形，部分颗粒边界呈锯齿状，主要发生亚颗粒旋转式(SR)动态重结晶，部分开始向颗粒边界迁移式(GBM)转变，动态重结晶作用较糜棱岩弱，粒径0.05~0.3mm之间，含量约55%；长石残斑表面较破碎，出现波状消光，发生塑性变形，部分呈眼球状、枣核状，发生膨凸式(BLG)动态重结晶，粒径0.1~2mm之间，含量约30%；黑云母含量较少(约10%)，多呈针状，以基质形式存在。
4.2 变质变形温度的估计 温压条件对矿物的变质变形机制以及岩石的变形类型具有重要影响，而温度对于矿物的变质变形程度影响最大，因此笔者根据变形岩石中石英、长石矿物的动态重结晶形式来估计变质变形温度。
据Sibson et al.(1979)、Passchier(1990)、van Daalen et al.(1999)及向必伟等(2007)的研究表明，石英的动态重结晶作用开始于280℃左右，在280~380℃时，表现为膨凸式重结晶(BLG)，颗粒边缘亚颗粒化；380~420℃之间，由膨凸式动态重结晶(BLG)向亚颗粒旋转式动态重结晶(SR)转变；420~480℃之间以独立的亚颗粒旋转式动态重结晶(SR)形式存在，且亚颗粒与残斑形成核幔构造；480~530℃之间由亚颗粒旋转式动态重结晶(SR)向颗粒边界迁移式动态重结晶(GBM)转变；530~630℃之间以独立的颗粒边界迁移式动态重结晶(GBM)形式存在，颗粒大小不规则，细粒化明显，颗粒边界呈港湾状、锯齿状。长石的动态重结晶形式与石英类似，指示温度范围有所差异。长石在低于300℃时以脆性破裂为主；300~400℃之间以显微破裂、波状消光为主；400~500℃之间呈塑性拉长及颗粒边缘亚颗粒化；500~600℃之间表现为BLG；650~700℃由BLG向SR转变；700~800℃以独立的SR存在；800~850℃由SR向GBM转变；850℃以上以独立的GBM存在。
研究区内出露的糜棱岩、片麻岩变形强烈，石英和长石矿物均发生动态重结晶，镜下清晰可见，便于估算温度。笔者对肥东北部文集地区所采集的岩石样品进行显微镜下观察分析，挑选出14块岩石标本，根据石英与长石矿物的动态重结晶形式对相应的变质变形温度做了估计，结果见表 3。
表 3
Table 3
表 3(Table 3)
表 3
肥东北部文集地区岩石变质变形特征与温度估计
Table 3
Characteristics of metamorphism and deformation and temperature estimation of rocks in Wenji area in the northern Feidong
样品号
岩石类型
矿物组合
动态重结晶类型
变质变形温度估计(℃)
石英
长石
WZ000
糜棱岩
Q+Fl+Bt
SR~GBM
BLG
500~600
WZ001-1
糜棱岩
Q+Fl+Bt
SR~GBM
BLG
500~600
WZ002-1
糜棱岩
Q+Fl+Bt
SR~GBM
BLG
500~600
WZ003-1
糜棱岩
Q+Fl+Bt
SR~GBM
BLG
500~600
WZ006-1
糜棱岩
Q+Fl+Bt
GBM
BLG~SR
530~630
WZ008-3
糜棱岩
Q+Fl+Amp+Bt
GBM
BLG~SR
530~630
WZ009-2
糜棱岩
Q+Fl+Amp+Bt
GBM
BLG~SR
530~630
WZ014-1
片麻岩
Q+Fl+Bt
SR~GBM
BLG
500~600
WZ016-1
片麻岩
Q+Fl+Bt
SR~GBM
BLG
500~600
WZ017-1
片麻岩
Q+Fl+Bt
SR~GBM
BLG
500~600
WZ018-2
片麻岩
Q+Fl+Amp+Bt
SR~GBM
BLG
500~600
WZ019-5
糜棱岩
Q+Fl
GBM
BLG~SR
530~630
WZ020-1
糜棱岩
Q+Fl+Bt
SR~GBM
BLG
500~600
WZ021-1
片麻岩
Q+Fl+Amp
SR~GBM
BLG
500~600
注：Q-石英；Fl-长石；Bt-黑云母；Amp-角闪石
表 3 肥东北部文集地区岩石变质变形特征与温度估计
Table 3 Characteristics of metamorphism and deformation and temperature estimation of rocks in Wenji area in the northern Feidong
结果表明，糜棱岩的变质变形温度高于片麻岩，综合剖面特征可知，肥东北部文集地区大背形核部岩石的变质变形温度高于两翼。
4.3 构造变形分析 韧性剪切带是变形和流体活动集中的地区，其内的变形岩石更是构造变形行为的载体。精细的显微构造变形研究包括分形、差异应力及应变速率等，不仅对认识剪切带内的构造变形环境具有约束意义，对研究区的构造演化、形成机制也具有重要的地质意义。
4.3.1 石英的分形特征 分形理论自引入地质科学以来，作为一种简单实用的工具被广泛应用于地质学的各个分支领域，在构造地质学领域主要用于确定韧性剪切带变形岩石的变形温度和应变速率等。Kruhl and Nega(1996)对不同级别变质岩和花岗岩中石英颗粒边界形态研究发现，石英动态重结晶的颗粒边界具有自相似性，表现出分形特征(王新社等，2001)，不同温度范围的石英颗粒边界的分形具有不同的分维数D。动态重结晶石英颗粒边界的分维数随颗粒边界形成温度的增大而减小，随应变速率的增大而增大(Kruhl and Nega, 1996；Takahashi et al., 1998)。
本次分形分析的研究对象为肥东北部文集地区剪切带内的糜棱岩及糜棱岩化片麻岩，采用面积-周长法(Lovejoy，1982)对石英动态重结晶颗粒进行分维数的计算。首先对重结晶石英颗粒进行显微拍照，利用Photoshop在显微照片上测定石英颗粒的真实周长P及面积A。每个样品统计的石英颗粒数均超过30个，通过A计算出等面积圆的直径d，即为粒径。由于周长和粒径遵循幂次定律，因此把每个样品的数据以周长的对数logP为纵轴，粒径的对数logd为横轴投于双对数图上(图 7)，对投影点进行最小二乘法拟合成一条直线，其斜率即为分维数D(王新社等，2001；吴小奇等，2006)。结果见表 4、图 7。
表 4
Table 4
表 4(Table 4)
表 4
肥东北部文集地区岩石动态重结晶石英颗粒边界分形特征
Table 4
Fractal characteristics of dynamically recrystallized quartz grain boundary of the rocks in Wenji area in the northern Feidong
样品号
测量数(N)
粒径分布d(μm)
平均粒径(μm)
周长分布P(μm)
平均周长(μm)
分维数(D)
相关系数(R)
WZ000
34
19.98~65.98
40.83
66.80~260.09
152.75
1.149
0.941
WZ001-1
39
38.63~118.48
72.96
133.23~517.45
275.28
1.143
0.985
WZ002-1
36
61.22~153.52
95.05
11.27~679.69
379.65
1.160
0.965
WZ003-1
37
60.93~302.37
131.02
205.75~1271.87
493.62
1.146
0.996
WZ006-1
41
44.33~126.57
71.62
155.56~690.49
303.52
1.169
0.945
WZ008-3
40
29.93~71.60
51.60
102.73~320.28
200.20
1.170
0.941
WZ009-2
34
43.14~82.03
63.81
149.82~340.44
242.48
1.151
0.953
WZ014-1
44
51.12~112.41
78.07
193.60~564.81
305.23
1.208
0.946
WZ016-1
34
93.24~252.10
158.96
317.55~1001.78
599.77
1.147
0.978
WZ017-1
41
133.99~336.50
212.39
468.86~1403.30
800.87
1.159
0.968
WZ018-2
35
124.61~281.00
193.15
424.84~1375.09
709.04
1.160
0.941
WZ019-5
34
29.50~96.32
69.97
102.88~425.88
263.45
1.147
0.970
WZ020-1
32
39.12~123.01
72.03
134.51~506.17
264.67
1.151
0.983
WZ021-1
38
142.19~333.01
204.11
508.67~1248.51
761.02
1.148
0.958
表 4 肥东北部文集地区岩石动态重结晶石英颗粒边界分形特征
Table 4 Fractal characteristics of dynamically recrystallized quartz grain boundary of the rocks in Wenji area in the northern Feidong
图 7
Fig. 7
图 7
肥东北部文集地区韧性剪切带动态重结晶石英颗粒的周长-粒径双对数图
d-重结晶石英颗粒粒径；P-重结晶石英颗粒周长；N-测量数；R-相关系数
Fig. 7
Log-Log plot of dynamically recrystallized quartz grain size against perimeter of the ductile shear zone in Wenji area in the northern Feidong
d-diameter of dynamically recrystallized quartz; P-perimeter of dynamically recrystallized quartz; N-number of measurements; R-correlation coefficient
结果显示，14个样品中动态重结晶石英颗粒边界统计数据的相关系数R均大于0.9，可见研究区内动态重结晶石英颗粒具有统计意义上的自相似性，且1 < D < 2(Kruhl and Nega, 1996)。分维数在1.143~1.208之间，并且随肥东北部几何学形态有规律变化。在大背形核部分维数D值较大，两翼渐小，近似对称性。颗粒粒径也有一定规律，大背形核部颗粒粒径较小，两翼渐大，近似对称性。这表明大背形核部的变质变形强度较强，由核部向两翼，变质变形作用减弱，岩石颗粒粒径变大，分维数D值变小。
值得注意的是：Takahashi et al.(1998)表明韧性剪切过程中变形温度不变的情况下，分维数随应变速率的增加而增大，因此，由石英的分形特征可以看出剪切带的变质变形强度变化特征。但Takahashi提出的高应变速率计所得数值较大，难以适合自然界较深层次形成的韧性剪切带(吴小奇等，2006)。不同变质程度的石英颗粒其分维数D值在1.05~1.30范围内，动态重结晶石英颗粒边界的分维数随着颗粒边界温度的升高而减小，可作为地质温度计使用：低绿片岩相变质岩中石英颗粒边界的分维数在1.23~1.31之间变化，高绿片岩相到低角闪岩相的分维数为1.14~1.23，麻粒岩相及同构造花岗岩中的分维数则为1.05~1.14(Kruhl and Nega, 1996)。分析结果发现，本次研究的分维数均位于高绿片岩相到低角闪岩相的分维数范围内，对应变质-变形温度在520~610℃之间，与利用矿物动态重结晶程度估算的变质变形温度在误差范围内基本一致。
4.3.2 差异应力 古差异应力值的计算对恢复变形构造层次具有重要的参考价值(吕承训等，2011)。显微构造测量差应力方法包括：自由位错密度统计法、亚颗粒大小法、机械双晶法及动态重结晶颗粒粒径法等(胡玲，2009；李海龙等，2017)。目前测量差异应力常用的方法是动态重结晶颗粒粒径法。前人通过对动态重结晶颗粒的分析发现，重结晶石英的颗粒粒径与差异应力存在对应关系，从而建立了推算古差异应力的重结晶颗粒粒径压力计。本文选用Twiss(1980)提出的石英动态重结晶颗粒粒径计算差异应力的公式：
(1)
其中，Δσ差异应力，单位为MPa；D为动态重结晶石英颗粒的平均粒径，单位mm；A和m是常数6.1和0.68。将表 4中计算分维值所测量的平均粒径带入上式，即可得到差异应力值为17.86~55.18MPa，结果见表 5。
表 5
Table 5
表 5(Table 5)
表 5
肥东北部文集地区韧性剪切带差异应力及应变速率
Table 5
Differential stress and strain rate of the ductile shear zone in Wenji area in the northern Feidong
样品号
测量数(N)
粒径分布(μm)
平均粒径(μm)
差异应力(MPa)
应变速率(s-1)
WZ000
34
19.98~65.98
40.83
55.18
3.470×10-12
WZ001-1
39
38.63~118.48
72.96
37.42
1.272×10-12
WZ002-1
36
61.22~153.52
95.05
31.55
8.200×10-13
WZ003-1
37
60.93~302.37
131.02
25.75
4.940×10-13
WZ006-1
41
44.33~126.57
71.62
38.00
4.330×10-12
WZ008-3
40
29.93~71.60
51.60
46.84
7.330×10-12
WZ009-2
34
43.14~82.03
63.81
40.10
4.790×10-12
WZ014-1
44
51.12~112.41
78.07
35.34
1.071×10-12
WZ016-1
34
93.24~252.10
158.96
22.01
3.177×10-13
WZ017-1
41
133.99~336.50
212.39
17.86
1.810×10-13
WZ018-2
35
124.61~281.00
193.15
19.06
2.145×10-13
WZ019-5
34
29.50~96.32
69.97
38.80
4.756×10-12
WZ020-1
32
39.12~123.01
72.03
38.04
1.331×10-12
WZ021-1
38
142.19~333.01
204.11
18.40
1.960×10-13
表 5 肥东北部文集地区韧性剪切带差异应力及应变速率
Table 5 Differential stress and strain rate of the ductile shear zone in Wenji area in the northern Feidong
结果表明，肥东北部文集地区剪切带内岩石的差异应力同样具有规律，与分维数类似，靠近大背形核部，差异应力值越大，两翼渐小，近似对称性。这表明，靠近大背形核部，岩石所受古差异应力越大，石英颗粒越小，变形越强。
4.3.3 应变速率 应变速率指的是单位时间内发生的线应变或剪应变。构造运动的应变范围很宽，应变速率可反映韧性剪切带的应变状态，同时会对矿物的动态重结晶机制产生影响。应变速率的估算理论上可以用分维数的应变速率计来计算，即Takahashi et al.(1998)的高应变速率计，但已有学者证明其计算结果普遍偏高(王新社等，2001；张秉良等，2008；梁琛岳等，2015)，自然界一般很难达到如此高的应变速率。传统的计算应变速率的方法为结合差异应力的方法，即Parrish et al.(1976)的湿石英流变速率公式最为可靠，公式如下：
(2)
其中，ε为应变速率，单位s-1；Δσ为差应力值，单位为MPa；T为温度，单位为K，温度采用表 3中各个样品变形温度估计值的中值温度，计算结果见表 5。
计算结果显示：肥东北部文集地区剪切带的应变速率ε范围为1.960×10-13~7.330×10-12s-1之间，一般区域性韧性剪切带的应变速率范围为10-13~10-15s-1，所以，肥东北部的应变速率与世界上大多数韧性剪切带应变速率一致，是缓慢变形的结果。对比研究区内应变强度、差异应力大小和应变速率与几何学形态关系可以发现，差异应力与应变速率呈正相关性，在大背形核部附近应变计速率大，两翼较小，近似对称性(图 3)。
5 结论 肥东北部文集地区位于郯庐断裂带中南段，岩石变形强烈，通过对研究区详细的野外观测和室内变形研究分析，本文主要得出以下几个结论：
(1) 肥东北部文集地区的构造形态为一走向北东-南西，枢纽向北东倾伏，轴面倾向南东的大型背形构造，其枢纽、轴面产状与区域性线理、面理产状基本一致，且面理产状西陡东缓。物质分布规律呈现出一定的对称性，核部物质以中-基性糜棱岩为主，夹同构造伟晶岩脉，两翼岩性渐中-酸性，呈现出一定的层状结构。郯庐断裂带肥东段具有A型褶皱的特征，其成因可能是郯庐断裂带剪切活动在大别造山带背形构造的基础上变形加强的结果。
(2) 应力场统计结果表明，肥东地区由于古太平洋俯冲作用影响，构造体主要受到近南北方向的挤压。根据涡度分析可知，剪切带是以简单剪切为主，纯剪切为辅，具压扭性质的左行韧性剪切带，所受力与剪切带有一定交角。
(3) 通过精细、系统的显微构造变形分析发现，肥东北部文集地区随几何学形态变化在变形温度、变形强度等方面也呈现出一定的对称性和规律：(a)区内岩石变形强烈，石英和长石矿物均发生动态重结晶，部分长石残斑可指示运动学方向为左旋。对石英和长石动态重结晶形式进行变形温度估计发现，肥东北部大背形的核部附近，岩石变形强烈，变形温度较高，两翼较低，近似对称性。(b)动态重结晶石英颗粒边界具有统计意义上的自相似性，分维数在1.143~1.208之间，位于高绿片岩相到低角闪岩相的分维数范围内，对应变质变形温度在520~610℃之间。自大背形的核部向两翼靠近，颗粒粒径和分维值也有规律的变化。核部分维数D值较大，两翼较小，核部颗粒粒径较小，两翼较大，近似对称性。这表明大背形核部的变质变形强度较强，两翼变质变形作用减弱。(c)差异应力和应变速率表现出一定的正相关性，石英粒径则与差异应力和应变速率表现出负相关性。在大背形核部附近差异应力和应变计速率大，两翼较小，近似对称性。表明从核部到两翼，差异应力和应变速率变小，差异应力越大，导致岩石受力越强，变形也越强，石英颗粒粒径越小。