1 Kunming University of Science and Technology, Kunming 650093, Yunnan, China; 2 The Third Institute of Geology and Mineral Exploration, Gansu Provincial Bureau of Geology and Mineral Exploration, Lanzhou 730050, Gansu, China 西秦岭夏河-合作地区是中国重要的金矿及多金属矿床开发地带，区内断裂构造复杂，构造控矿作用显著。文章利用分形理论方法对夏河-合作地区不同方向断裂构造作定量分析，统计计算断裂构造容量维和信息维，结合矿床空间分布分形特征和矿床Fry分析，探讨研究区构造分形特征及断裂与矿床的空间展布关系。研究结果表明，研究区全部断裂容量维1.2374，北东向断裂容量维0.8229，北西向断裂容量维1.2296，近东西向断裂容量维1.0312，全部断裂容量和北西向断裂容量维处于断裂分形临界值(1.22~1.38)区域，显示研究区断裂构造连通性较好。研究区全部断裂信息维1.1118，北东向断裂信息维1.018，北西向断裂信息维1.2296，近东西向断裂信息维1.0217，北西向断裂信息维处于断裂分形临界值区域，反映北西向对矿床的控制最强。断裂构造分形维数可以表征矿床的沉积部位，以容量维、信息维等值线图以及Fry图划分出3个有利成矿区。 The Xiahe-Hezuo area in the West Qinling Mountains is an important mining area for gold and polymetallic deposits in China, with complex fault structures and significant structural ore control. In this paper, the fractal theory method is used to quantitatively analyze the fracture structure in different directions in the Xiahe-Hezuo area, statistically calculate the fracture structure capacity and information dimension, combining with the fractal characteristics of the spatial distribution of the deposit and the Fry analysis of the deposit, the fractal characteristics of the structure of the study area and the spatial distribution relationship between the fracture and the deposit are discussed. The results show that the total fracture capacity dimension of the study area is 1.2374, the north-east fracture capacity dimension is 0.8229, the north-west fracture capacity dimension is 1.2296, and the near-east-west fracture capacity dimension is 1.0312, and the total fracture capacity and the north-west fracture capacity dimension are in the fracture fractal critical value (1.22~1.38) region, indicating that the fracture tectonic connectivity of the study area is good. The information dimension of all faults in the study area is 1.1118, the information dimension of the north-east fault is 1.018, the information dimension of the north-west fault is 1.2296, the information dimension of the near-east-west fault is 1.0217, and the information dimension of the north-west fault is in the fracture fractal critical value area, reflecting the strongest control of the deposit in the north-west direction. Fracture structure fractal dimension can be used to characterize the sedimentary site of the ore, and three favorable mineralization areas are identified by capacity dimension, information dimension contour map and Fry diagram. 分形理论是指具有自相似性的现象、图像或者物理过程（如羊齿植物、菜花、海岸线），由美籍数学家曼德尔布罗特于1975年描述碎石过程中创造出，认为局部是整体的比例缩小图像（Mandelbrot et al.,1967;1991）。对于无法用欧几里得和微积分刻画的自然界分形客体，分形理论能揭示其规律性、层次性、自相似性和标度不变性（孙文洁，2022）。近几十年，分形理论在多学科多领域得到广泛应用，特别是在地质学领域涌现大量研究成果，成矿预测（王喜生等，2000；赵少攀，2015；李飞，2016；王维，2016；孙涛等，2017；2018；崔中良，2022；Cui et al.，2022）、矿床品位分布（万丽等，2007；谢焱石等，2010）、断裂构造特征（廖家飞等，2012；崔中良，2021）、储层孔隙结构（李昊远等，2023；Gu et al.，2023）等为分形理论打下坚实的理论基础。 断裂构造是某次地质运动过程中地层错动产生位移的地质现象，在空间形态上呈现一系列不规则线状或面状分形集合体，同海岸线一样，断裂构造具有自相似性（谢和平，1994）。分形理论研究方法将断裂构造研究从定性研究阶段提升到定量研究阶段（朱剑兵等，2005）。前人研究表明，断裂构造分形维数能反映断裂发育的复杂程度，有效表征断裂发育特征，分形维数越大，断裂构造越发育（范存辉等，2017）。有学者通过岩石块体二轴压缩试验测得断裂构造分形维数临界值为1.22~1.38（Kruhl，1994），该临界值可以评价断裂系统连通性，断裂分维值低于临界值则表明断裂连通性较差，断裂孤立，分维值处于或高于临界值，则断裂连通性好，有利于成矿流体的迁移、聚集和赋存（卢新卫等，1999;谢焱石，2010）。因此，区域断裂构造分形维数值能定量表征断裂构造发育程度，有效阐明断裂构造对成矿流体的控制作用。 西秦岭成矿带是中国重要的多金属成矿带，已发现矿产44种，矿产地675处（姜寒冰等，2023）。西秦岭这种成矿优越性吸引着大量学者，在其矿产资源评价、成矿系列以及找矿预测方面做了大量研究工作（李康宁等，2020；何进忠等，2023；王雷等，2023）。位于西秦岭中部的夏河-合作地区地质构造复杂，成矿地质条件优越，矿床呈现“北西成带，北东成串”的特征（李康宁等，2020），其中的金矿床最为发育，早子沟、加甘滩、以地南等大型-特大型金矿极具代表性。区内年代学、矿床成因、找矿预测等方面针对单个矿床做了大量研究工作（Jin et al.，2017；Yu et al.，2019；Du et al.，2021；第鹏飞等，2021；Yu et al.，2022），但未见使用分形理论方法定量评价断裂构造与矿床空间分布规律的相关报道。文章以夏河-合作地区断裂构造和矿床分布特征为研究对象，定量研究区内断裂构造分形维数与矿床空间分布的关系，结合Fry分析法，以期为该区成矿预测提供理论依据。 1地质背景 夏河-合作地区位于青藏高原东北缘，隶属于秦岭弧盆系，北接中祁连地块，南邻若尔盖古陆，处于古亚洲构造域、特提斯构造域和滨太平洋构造域复合叠加的构造部位（图1a、b，张国伟等，2004;路英川，2016）。海西期地质运动造成地层逐步拉张裂陷，以及印支期的强烈挤压造山（姜寒冰，2023），造成区域构造线呈北西-南东向展布，褶皱、断裂发育，地质构造复杂，区域成矿条件优越。 区域出露地层主要有二叠系和三叠系。其中， 二叠系分布于区内东北部，呈北西向-南东向展布， 与三叠系果木沟组呈断层接触， 主要为上二叠统石关组二段及石关组三段;三叠系在区内呈北西向-南东向展布，与区域构造线方向一致， 区内出露下三叠统果木沟组和江里沟组，为区内主要的金矿赋矿地层。 区域上主体褶皱为新堡-力士山复背斜，地层强烈变形， 次级尖棱褶皱、倒转褶皱及层间板、劈理十分发育（韦良喜，2018）。背斜的核部和向斜的两翼有较多的石英脉、方解石脉及花岗闪长岩脉侵入。区域断裂构造普遍发育，已发现的金矿床和金矿（化）点的区域性展布严格受北西向断裂带及其次级断裂的控制，控矿断裂主要为夏河-合作断裂带和桑科南-格里那断裂带，2个断裂之间发育4条北西向次级主干断裂，由北向南分别是麦来梁-早仁道、扎油梁-早子沟、桑曲-地娄塘、索拉贡玛-加甘滩压性、压扭性逆冲断裂，断层延长35~50 km，平面上呈舒缓波状展布，倾向北东，倾角45°~70°。在北西向主干断裂带上局部发育北东向、近南北向断裂，断层规模较小，延长0.5~2.0 km，沿断裂带有岩脉、石英脉等充填，与成矿关系十分密切，如早子沟金矿的主矿体就产于北东、北西、近南北向断裂破碎带中。区域性深大断裂带主要起导矿的作用，其次级断裂为主要的容矿构造。 区域岩浆活动较频繁，侵入岩普遍发育，主要有桑曲、马九勒花岗闪长岩，将其那梁、塔日干、早仁道石英闪长岩，杂恰勒布、扎油梁英云闪长斑岩等规模较大的侵入体及各种岩脉（梁志录，2016）。侏罗纪岩浆活动与金矿的就位关系密切（刘春先等，2018）。 2 断裂构造分形维数特征 2.1 容量维D ₀ 2.1.1 断裂构造容量维特征 断裂构造在空间展布上呈无规则而复杂的线状或面状形态，整体呈自相似性，可用分形维数对其描述（姜颜良，2023）。分形维数存在多种计算方法，如盒计维数、滑动窗格-模型、P-A（周长-面积）模型、N-A（数量-面积）模型等，因盒计维数法简单便捷，使用该方法对研究区断裂构造系统进行维数计算，从而分析研究区断裂构造的变化规律及其对矿床的控制作用，运算过程如下：采用不同边长为r的二维正方形网格依次覆盖研究区，统计不同边长为r的网格下包含断裂的正方形数目N（r），N（r）和r存在以下幂定律关系（谢焱石，2002）： 在实际计算过程中，分别以r为50 mm（实际距离11.1111 km）、25 mm（实际距离5.5555 km）、12.5 mm（实际距离2.7777 km）、6.25 mm（实际距离1.3888 km）、3.125 mm（实际距离0.6944 km）的正方形网格覆盖研究区（图2），并分别统计包含全部断裂、北西向断裂、北东向断裂以及近东西向断裂的正方形网格数，最后做线性回归，计算出断裂构造分维值和判定系数平方（R ² ）。统计及计算结果见表1，lnN（r）-lnr最小二乘拟合直线见图3。 夏河-合作地区总体断裂构造分形维数为1.2374，处于能反映断裂系统连通性的临界值1.22~1.38区（图3），表明研究区断裂构造系统连通性较好，有利于成矿流体的迁移、聚集和成矿。判定系数平方（R ² ）能反映分形分维特征的优良性（孙文洁，2022），全区断裂构造判定系数平方（R ² ）为0.9952，大于0.99，表明研究区具备良好的分形分维特征。区内共发育3个方向的断裂构造，发育程度各不相同。其中，北西向断裂构造最为发育，全区皆有出露，分形维数为1.2385，相关系数0.9963，断裂连通性较好，有利于矿液迁移聚集，完尕-扎不浪可、和日卡加-唐四尔、来周村-早仁道、桑曲-完肯、尕恰勒布-完安囊等沿线矿床均沿北西向-南东向断裂分布；北东向断裂构造发育较少，主要出露于研究区北东地界，分形维数0.8229，判定系数0.9961，对北西向断裂相互联通起重要作用;近东西向断裂构造主要分布于研究区北东侧，少部分近东西向断裂是部分北西向断裂的延伸转向，分形维数1.0312，判定系数0.9982，对矿床的控制较差。 分区容量维计算过程如下：以r=50 mm（实际距离11.1111 km）的正方形网格覆盖研究区，将研究区平均分成40个子网格，将各子网格依次编号（图2）。采取计算容量维的相同研究标度去计算分区容量维，即用25 mm（实际距离5.5555 km）、12.5 mm（实际距离2.7777 km）、6.25 mm（实际距离1.3888 km）、3.125 mm（实际距离0.6944 km）的正方形网格覆盖研究分区子网格，统计各标度范围内子网格中包含断裂的网格数N（r）。同样的计算lnr与lnN（r），并做线性回归计算出分形维数和判定系数（R ² ）。 统计计算结果显示，除去未包含断裂构造的子网格，分形维数最小值为0.6584，最大值为1.6051（表2）。判定系数最小值为0.9561，最大值0.9995，普遍高于0.99，拟合度都比较高，具良好的分形维数特征。 2.2 信息维D ₁ 2.2.1 断裂构造信息维特征 信息维的研究方法与容量维类似，不过考虑的重点是出现在某一网格中的断裂数目（概率），也称断裂强度。某一标度为r的正方形网格将断裂构造分割成N（r）份，断裂构造出现在第i个网格中的概率为Pi（r），计算公式为式（3），信息量则为I（r），计算公式为式（4）（谢焱石，2002;崔中良，2022）。 在实际计算过程中，同样采取与容量维计算尺度相同的正方形网格覆盖研究区，正方形边长r为50 mm（实际距离11.1111 km）、25 mm（实际距离5.5555 km）、12.5 mm（实际距离2.7777 km）、6.25 mm（实际距离1.3888 km）、3.125 mm（实际距离0.6944 km）。分别计算出研究区总体断裂、北东向断裂、北西向断裂以及近东西向断裂的信息量I（r），统计数据见表3，利用excel工具作I（r）-lnr线性回归图（图4），根据直线斜率得出信息维D ₁ 。 计算结果显示，全区断裂构造信息维1.1118，判定系数0.9894，北东向断裂构造信息维1.018，判定系数0.9796，近东西向断裂构造信息维1.0217，判定系数0.9809，北西向断裂构造信息维1.2296，判定系数0.9917。北西向断裂信息维和判定系数最高，表明北西向断裂最为复杂且具有较强的相似性，该方向的断裂也是主要的控矿断裂，这与研究区以北西向断裂最为发育，且矿床沿断裂成带分布相吻合。 2.2.2 分区断裂构造信息维特征 信息维的分区计算标度与容量维标度相同，同样以边长r为50 mm（实际距离11.1111 km）的正方形覆盖研究区，共划分成40个子区，用25 mm（实际距离5.5555 km）、12.5 mm（实际距离2.7777 km）、6.25 mm（实际距离1.3888 km）、3.125 mm（实际距离0.6944 km）标度的正方形网格覆盖子区，用前面计算信息维的方法计算出分区信息维，这部分不再赘述。统计计算结果见表4，2、12等8个断裂空白分区未被统计，分区信息维0.5341~1.5817，平均值1.1081，且分区判定系数（R ² ）多数大于0.99，分区拟合度都较高，断裂有较强的统计自相似性。 三叠纪以来的印支运动使夏河-合作地区构造、岩浆活动强烈，各类矿产大量产出（李康宁等，2020），矿床成因类型主要为岩浆热液型和矽卡岩型。在分析矿床的空间分布分形特征过程中，将同一成因类型、成矿时代的矿床归为一组，分别计算其分维数。用边长r为50 mm（实际距离11.1111 km）、25 mm（实际距离5.5555 km）、12.5 mm（实际距离2.7777 km）、6.25 mm（实际距离1.3888 km）的正方形网格覆盖研究区，分别统计包含同一成因类型、成矿时代的矿床出露点的网格数N（r），分形维数计算过程同计算容量维，不过将断裂构造转变为矿床出露点。 统计结果见表5，分形维数计算见图5。岩浆热液型矿床分形维数0.4352，矽卡岩型矿床分形维数0.3967。侏罗纪矿床分形维数0.3363，三叠纪矿床分形维数0.4486。 其中，d（r）为以某个点为圆心，半径为r的圆限定的区域中的点的密度，D为密度分形维数，K为常数。本次计算过程选择研究区枣子沟金矿、以地南金矿、加甘滩矿、隆瓦寺院矿床作为原点，以半径r为6 km、12 km、18 km、24 km、30 km的圆覆盖研究区，统计圆内全区金属矿床以及金矿的数量，并计算矿床的密度（表6）。 根据计算公式（6）计算分形维数。研究区总体金属矿床和金矿床平均数量分形维数分别为1.5889、1.6384，判定系数（R ² ）分别为0.9961、0.9915;全区总体金属矿床和金矿床平均密度分形维数分别为1.589、1.638，判定系数（R ² ）分别为0.9437、0.9559，反映矿床沉积密度较高，成矿流体易聚集成超大型矿床（图6）。 Fry分析是一种对点目标进行空间自相关分析的几何方法，一开始是研究岩石中的应变和应变分配（Fry et al.，1979）。随着技术理论的深入研究，该分析方法逐步应用于探讨矿床的空间展布规律，揭示矿床沉积点位的空间位置和联系（Zhao et al.，2023），以及断裂构造与矿床的产出关系（Carranza，2009；Carranza et al.，2010；Parsa et al.，2018; Zheng et al.，2020）。计算Fry点图中每15°方向范围内矿床的数量，并绘制相应的玫瑰花图，可以用来解释矿床出露点分布的方向性，具体分析步骤为：①研究区矿床出露点被视为原始点集，选择中心部位的一点作为原点；②复制原始点集，将点集中某点覆盖在原点上，保留点集位置；③重复以上步骤，每个点都有一次机会覆盖在原点上。最后得到（n ² ~n）个点，称之为Fry图，并根据Fry图制作矿点玫瑰花图，从而分析矿床的方向性。 Fry图（图7）显示，矿点呈北西走向，与研究区断裂构造走向一致，表明北西向断裂构造对矿床的控制强烈。通过矿点趋势，划分出2个有利成矿区，该分析结果与前人分析成果吻合。 本次以研究区Fry图为底图，选择中心位置为基准点，计算统计每15°区间范围内金矿床、铜矿和其他矿床的数量，并制作玫瑰花图（图8）展示出研究区矿点方向性。研究区金矿床的沉积趋势为北西-南东、北东-南西向，但分形维数显示，研究区南西部成矿几率较低。铜矿优势方位为北西，其他金属矿床则呈各向发育。 3 研究区成矿预测 分形维数对矿床的定位有一定的指示意义，断裂构造通常是成矿流体的流通通道，也是赋存空间，结合研究区容量维和信息维作分形维数等值线图（图9a、b），判断研究区有利的成矿部位，为找矿预测提供理论指导。 断裂容量维等值线图（图9a）显示，容量维大于1.22的区域主要分布于研究区中部至东部（红色区域），整体呈北西走向串珠状，容量维高值区正好处于以地南金矿床出露部位。容量维值高低过渡区位于研究区中部，已探测到大量金矿矿床，其中以枣子沟、加甘滩等特大型金矿最为著名。容量维高值区地质构造复杂，连通性较好，能保证成矿流体远距离流动，在构造连通性减弱时沉积，该区也是有利成矿区。研究区南西侧容量维值较低，断裂构造较单一且孤立，能圈闭住成矿流体，为高值区输送来的成矿流体于容量维高低过度区沉积提供堡垒。 本次将容量维和信息维高值区、高低过度区划分为有利成矿区。研究区矿床主要分布于区内中部，整体呈“北西呈带，北东成串”的特征，分析结果符合实际勘探结果，证明断裂构造分形理论进行区域成矿预测具有一定可行性。 以分形维数临界值1.22~1.38为划分标准，综合容量维等值线图、信息维等值线图以及Fry成矿预测图相互叠加划分出有利成矿区（图10）。 区域Ⅰ：该区域为重点找矿靶区，断裂构造容量维和信息维皆处于或大于分形维数临界值，也是Fry分析有利成矿区，断裂构造在该区域发育且复杂，有利于矿床的富集，成矿潜力巨大，已探明以地南特大型金矿床及众多其他多金属矿床，根据分形维数特征，预测该区域深部、断裂交叉部位以及次级断裂部位存在隐伏矿床的可能性较大，预测结果与前人分析结果一致（李康宁等，2020）。 区域Ⅱ：大多处于Ⅰ区的边部，属于容量维值高值区、信息维值高值区、Fry分析有利成矿区两两重合区域，找矿潜力次于区域Ⅰ。 区域Ⅲ：此区域包含有容量维值高值区、信息维值高值区、Fry分析有利成矿区其中之一，该区域在研究区分布范围最大，特别是在研究区中部地带，研究区南西部分形维数低值区对成矿流体形成圈闭作用，造成区域Ⅲ同样找矿潜力巨大。 4 结 论 （1） 通过夏河-合作地区断裂构造分形维数统计分析，整体断裂构造容量维1.2374，信息维1.1118，具有较强的统计自相似性。北西向断裂容量维1.2385，信息维1.2296，反映北西向断裂对矿床的控制最强。 （2） 夏河-合作地区岩浆热液型和矽卡岩型矿床显示出单一分形分布特征，反映成矿物质来源单一。总体金属矿床和金矿床平均数量分形维数分别为1.5889、1.6384，平均密度分形维数分别为1.589、1.638，指示了矿床沉积密度高、数量多和规模大的特点，区内找矿潜力巨大。 （3） 断裂构造分形维数可以表征矿床的沉积部位，通过夏河-合作地区容量维、信息维等值线图以及Fry图综合分析，圈定出3个有利成矿区域，区域Ⅰ最有利成矿，区域Ⅱ和区域Ⅲ次之。 致 谢 研究工作得到了甘肃省合作早子沟金矿有限责任公司仲海书、殷超、田向盛等提供资料和野外工作上的支持，审稿专家提出宝贵的修改意见，在此一并致以诚挚的感谢! 主管单位：中国科学技术协会 主办单位：中国地质学会矿床地质专业委员会 中国地质科学院矿产资源研究所
地　　址： 北京市百万庄大街26号 邮编：100037 电话：010-68327284;010-68999546 E-mail： [email protected]
本系统由北京勤云科技发展有限公司设计