Geochronology and geochemistry of Mashan trachydacite, JiMo District, Shandong Province and their geological implications
-
摘要: 山东即墨马山地区出露罕见的发育柱状节理的中酸性火山岩，对于其年代学及形成的大地构造背景鲜有讨论。岩石薄片与主量元素分析结果证明马山柱状节理火山岩属于粗面英安岩。利用LA-ICP-MS对马山柱状节理火山岩中的锆石进行U-Pb测年，确定其形成时间为113.2±1.3 Ma，属于早白垩世晚期。原位Lu-Hf分析获得锆石εHf(t)值范围为−19.5～−0.4，对应的二阶段模式年龄范围约为1.97～3.26 Ga，暗示其源区岩浆成生于古老地壳的熔融。马山粗面英安岩具有SiO2和Al2O3含量较高，MgO含量低及富含碱的主量元素的特征。其稀土元素总量较低，配分模式呈轻稀土富集、重稀土亏损的右倾模式，Eu轻微负异常，Sr含量高而Y和Yb含量低。马山粗面英安岩的主微量元素显示其具有埃达克岩特征，推测其来源岩浆可能形成于燕山期华北东部区域性地壳增厚后下部老地壳的熔融。Abstract: Intermediate-acid volcanic rocks with columnar joints expose in the Mashan area, Jimo, Shandong Province. The geochronology and tectonic setting of the volcanic rocks had been rarely studied. The petrographic characteristics and major element results show that the volcanic rocks at Mashan are trachydacite. Zircon LA-ICP-MS U-Pb dating of the trachydacite shows that the volcanic rocks erupted at 113.2 ± 1.3 Ma, that is, in the late Early Cretaceous. In-situ Lu-Hf analysis suggests that the zircon εHf (t) values range from −19.5 to −0.4 and two-stage model ages of zircons range from 1.97 Ga to 3.26 Ga, indicating that the magma was sourced from the melting of ancient crust. The Mashan trachyandesite is characterized by high content of SiO2 and Al2O3, low content of MgO, enriched Na2O and K2O. Moreover, all samples are featured by lower total contents of rare earth elements, enriched light rare earth elements, deficit heavy rare earth elements, slightly negative Eu anomaly, high Sr content, low Y and Yb contents. The major and trace elements of the trachydacite in Mashan are quite similar to adakite, and the adakitic magma could have been formed by the melting of the lower old crust after the regional crust thickening in Eastern Norch China Block during the Yanshanian.
-
Keywords:
- trachydacite /
- zircon U-Pb /
- Lu-Hf /
- Early Cretaceous /
- Adakite /
- Mashan
-
沉积物粒度组成及其分布模式是沉积物的基本特征，能够揭示沉积物形成的动力条件，是判别沉积环境和沉积物源的重要指标。不同沉积环境形成的沉积物受其搬运方式、动力条件等的约束，可分选出特定的粒级组合，因此粒度分布特征分析具有重要的沉积学意义[1]。沉积物粒度频率分布曲线可以定量反映某一粒径分布含量、次总体类型、分选程度等沉积过程的细节信息[2]。这些信息的提取有助于理解沉积物的来源、成因及形成过程，尤其是粒度分布子体分离法[3]、Weibull分布函数拟合法[4]，端元模型法[5]等统计方法的应用，对于研究不同粒度组成端元或者组分的气候环境意义具有重要的推动作用。上述统计学方法重点关注某一沉积类型沉积物形成的动力过程和沉积成因，对于沉积环境复杂、沉积类型多样的沉积区，其沉积物成因相对复杂，如何快速、准确判别沉积物的沉积相显得尤为重要。聚类分析作为数据挖掘领域重要的手段之一，直接比较数据之间的性质差异，将数据性质相同或相近的聚为一类，在沉积物粒度数据分析与处理中具有较好的适用性[6]。然而，常规的聚类分析方法需要预设聚类数量，不能真实地反映根据沉积物粒度分布曲线类型反演出研究区存在的沉积动力类型，而且以往研究应用聚类分析大多以粒度参数表征粒度频率分布曲线作为研究对象，不能精确地表示粒度分布曲线形态以及各粒径分布含量[7]。值得注意的是，AP聚类算法（Affinity Propagation Clustering Algorithm）作为新兴的聚类方法，具有快速、高效、不需要预设聚类数量的优点，对于高维数据聚类效果明显，在图像、文本聚类以及基因表达数据分析等领域应用广泛[8-10]，将AP聚类算法方法引入沉积学研究中，对于粒度数据的深度挖掘与分析具有重要的推动作用。
白洋淀地区分布深厚的第四纪沉积物，沉积环境较为复杂，以冲积、洪积和湖积作用为主，是开展不同沉积环境条件下沉积物粒度分布曲线判别与划分研究的理想区域。区内主要侧重于利用气候环境代用指标进行古环境演化方面的研究，缺乏沉积环境动力条件的研究[11-13]。通过采集研究区内沉积剖面中已知沉积相的沉积物，确定不同沉积环境的沉积物粒度特征，研究AP聚类算法在沉积学研究中的可行性与适用性，尝试建立沉积物粒度特征与沉积动力条件的内在联系。
1. 自然地理概况
研究区位于华北平原白洋淀及其周边区域，是永定河冲积扇与滹沱河冲积扇的前缘低洼区，位置为38.7°～39.2°N、115.6°～116.5°E。研究区属于温带大陆性半湿润半干旱气候，年平均温度为7.3～12.7 ℃，多年平均降水量564 mm，主要集中在6—8月份[14]。区内地势平坦，由西向东地势略有降低。白洋淀有拒马河、府河、唐河等九条河流汇入以及赵王新河、大清河等出淀河流。第四纪以来，该地区不断沉降，并普遍接受沉积，而且区内水文条件的变化对第四纪沉积物的形成也具有重要的影响，逐渐形成了一套以冲洪积和湖积为主的沉积地层[15]。
2. 材料与方法
2.1 样品采集与测试
本研究对白洋淀周边地区露头剖面的沉积地层进行了详细考察，选取其中22个剖面作为研究对象（图1，剖面地层描述见表1），对各剖面典型沉积相层位进行样品采集，共获得85个沉积样品。现代河流漫滩相沉积样品来源于白沟引河，共采集2个沉积样品，沉积物为棕黄色黏土质粉砂，具有典型的二元结构。现代湖泊表层沉积样品来源于白洋淀湖泊，共采集12个沉积样品，沉积物为灰黑色粉砂质黏土、黑色黏土，有机质含量高（图2）。
表 1 采样点描述Table 1. Description of sampling sites采样剖面名称 剖面位置 剖面描述 采样数/个 新庄窠村 38°59′31.56″N，115°50′54.04″E 沉积物以灰黄色黏土质粉砂、灰黑色或灰色粉砂质黏土和浅黄色细砂为主，沉积环境主要为湖滨相、湖沼相和河床相 3 段庄村剖面 39°03′20.44″N，115°49′30.28″E 沉积物为灰黑色粉砂质黏土、灰黄色黏土质粉砂，沉积环境为湖沼相、湖滨相 3 北张村剖面 39°04′40.23″N，115°46′32.13″E 沉积物以灰黑色、棕黄色粉砂质黏土和灰黄色细砂为主，沉积环境为湖沼相以及河床相 4 向村剖面 39°00′05.42″N，115°56′46.51″E 沉积物以灰黄色、浅灰色黏土质粉砂、灰褐色粉砂质黏土、灰色细砂以及灰黄色粉砂质中砂为主，沉积环境主要为湖滨相、湖沼相和河床相、洪积相 4 大阳村 39°00′18.32"N，115°59′19.28"E 沉积物为灰黑色粉砂质黏土、浅黄色黏土质粉砂以及具有水平层理的灰黄色粉砂，沉积环境主要有湖沼相、漫滩相以及河床相 3 李庄头村 38°59′33.60"N，116°04′48.67"E 沉积物以细砂、中砂为主，沉积环境主要为河流相 3 晾马台村 39°04′14.65"N，115°59′05.65"E 沉积物具有典型的“二元结构”，沉积环境为漫滩相 1 北剧村 09°04′56.42"N，115°58′16.81"E 沉积物以浅灰色中砂为主，沉积环境为河床相 1 南文村 39°04′19.98"N，115°56′42.30"E 沉积物以黄褐色、灰黄色黏土质粉砂以及含有黑色碳屑的灰黄色粉砂质粗砂为主，沉积环境主要有洪积相、湖相 4 马家庄村 39°02′53.77"N，115°55′07.67"E 沉积物为灰黄色黏土质粉砂，沉积环境为湖滨相 2 胡家台村 39°02′25.71"N，116°07′02.92"E 沉积物以灰黄色粗粉砂、细砂为主，发育交错层理，沉积环境为漫滩相 3 大步村 39°00′31.04"N，116°10′27.99"E 沉积物为灰褐色、黑色黏土、褐黄色黏土质粉砂以及浅灰色中细砂，沉积环境为湖沼相、湖滨相、河床相 8 邢岗村 39°03′0 8.65"N，116°14′18.92"E 沉积物以棕红色、灰黄色粉砂质黏土以及灰色、灰黄色黏土质粉砂为主，沉积环境为湖沼相、湖滨相 10 袁郭村 38°51′03.40"N，116°25′15.25"E 沉积物以灰黑色粉砂质黏土、浅黄色黏土质粉砂、灰黄色中细砂为主，沉积环境为湖沼相、湖滨相 4 胡屯村 38°51′59.62"N，116°22′34.97"E 沉积物以灰褐色、灰黑色粉砂质黏土以及灰黄色粉砂为主，沉积环境主要为湖沼相、湖滨相 5 小务村 38°51′8.86"N，116°17′39.73"E 沉积物以棕黄色亚砂土、细粉砂、灰黑色粉砂质黏土为主，沉积环境为湖滨相、湖沼相 6 双塔村 38°45′26.20"N，116°1′44.33"E 沉积物以灰黄色细粉砂、灰黑色黏土质粉砂为主，沉积环境为湖滨相、湖沼相 4 东良淀村 38°44′5.70"N，
115°59′7.99"E沉积物为棕红色粉砂质黏土，沉积环境为湖相 2 北队村 38°46′22.53"N，115°57′27.94"E 沉积物为浅灰色细砂土、棕黄色黏土质粉砂，沉积环境为漫滩相、湖滨相 3 韩堡村 38°46′23.07"N，115°53′24.07"E 沉积物为灰黄色黏土质粉砂，沉积环境为湖滨相 2 东垒头村 38°52′16.08"N，115°54′02.73"E 沉积物以灰黑色、棕红色粉砂质黏土、灰黄色黏土质粉砂为主，沉积环境为漫滩相、湖滨相、湖沼相 7 董庄村 38°49′32.73"N，
115°46′06.36"E棕红色、灰黑色粉砂质黏土，沉积环境为湖沼相 3 粒度测试在中国地质科学院第四纪年代学与水文地质演变重点实验室完成，采用mastersizer2000激光粒度仪进行测试，该仪器的测试范围0.2～2 000 μm。测试前加入足量的10%双氧水并加热，去除有机质。待烧杯中不再冒泡之后，加入足量的10%盐酸，直至不再冒气泡才停止加热，然后加入蒸馏水静止24 h。利用虹吸法抽取上层清液，加入蒸馏水，静止24 h，反复几次直至呈中性。加入10 mL 10%的六偏磷酸钠，用超声波振荡10 min后上机测试。
2.2 邻近传播算法（AP聚类算法）
邻近传播聚类算法又称AP聚类算法，是B J Frey 和 D Dueck 于 2007 年首次提出的一种新的聚类算法[16]。它是一种不需要预设聚类数，通过数据点之间的“信息传递”来寻找自己的类代表点（exampler），将属于同一类代表性点的样品聚成一类，从而实现聚类过程[17]。
AP聚类算法将所有的样品点都看作为潜在的类代表性点，以数据之间的相似度系数基础，两两之间进行“信息传递”。在信息传递的过程中，利用因子图（图3）中The Max-Sμm algorithm准则进行信息的迭代更新，迭代收敛时，使得公式（1）在约束条件下网络相似度S（c）达到最大。
$$\begin{aligned} &S\left( c \right) = \mathop \sum \nolimits_{i = 1}^N s\left( {i,{c_i}} \right) + \mathop \sum \nolimits_{k = 1}^N {\delta _k}\left( c \right)\\&{\text{约束条件}}:\;\begin{array}{*{20}{l}} {{\delta _k}\left( c \right) = } \end{array}\left\{ {\begin{array}{*{20}{l}} { - \infty ,\;\;\;{\text{当}}{c_k} \ne k{\text{但是}}\exists i:{c_i} = k}\\ {0,{\text{否则}}} \end{array}} \right. \end{aligned}$$ (1) 式中，c为由类代表性点构成的向量，i为样品点，s（i，ci）表示样品点与其对应的类样品点之间的相似度；δk（c）是约束条件，表示数据点i选取k作为它的类代表性点，但是数据点k并不是类代表性点，因此使其值为−∞，用于限制没有类代表性点的类簇产生。
相似性数据矩阵s（n，n）是由n个数据之间的相似度系数构成，对角线上的元素由人为设定的偏向参数所构成，它是用于评价各数据点作为类代表性点可能性的程度。数据之间传递的信息包括：从潜在的类代表性点传递到数据点的归属度和从数据点传递到类代表性点的吸引度[16]。在因子图中，基于The Max-Sμm algorithm准则信息更新的公式：
$$ r\left( {i,k} \right) = s\left( {i,k} \right) - \mathop {\max }\limits_{k' \ne k} \left( {s\left( {i,k'} \right) + a\left( {i,k'} \right)} \right),\forall i,k $$ (2) $$ a\left( {i,k} \right) = \left\{ {\begin{array}{*{20}{l}} {\displaystyle\mathop \sum \nolimits_{i':i' \ne k} \max\left\{ {0,\left( {i',k} \right)} \right\},\forall k = i}\\ {\min \left\{ {0,r\left( {k,k} \right) + \displaystyle\mathop \sum \nolimits_{i':i' \ne i,k} \max\left\{ {0,r\left( {i',k} \right)} \right\}} \right\}},\\\forall k \ne i \end{array}} \right. $$ (3) 当信息迭代到更新收敛时，点i所对应的类代表点ci：
$$ {c_i} = \arg \max \left\{ {a\left( {i,k} \right) + r\left( {i,k} \right)} \right\},\;\forall i $$ (4) 基于AP聚类算法的基本思想及其原理，可得算法的基本步骤：
第1步：根据样品点之间的相关距离，计算相似的矩阵S以及设置各个样品点对应的偏向参数。
第2步：设置各个样品点的初始归属度
$a\left( {i,k} \right) = 0$ 。第3步：利用公式（2）、（3）进行信息更新，在更新的过程中，为了防止迭代震荡，引入阻尼因子λ，范围（0，1）。
第4步：当各个点对应的类代表性点ci在迭代过程中不再发生改变或者迭代次数达到设定的最大值时，迭代终止。将相同类代表性点的样品点归为一个类簇。
3. 沉积物粒度分布曲线的AP聚类算法
沉积物粒度分布曲线特征是判断沉积作用类型的重要手段之一。通常，研究者采用粒度参数代替粒度频率分布曲线，近似地作为沉积作用类型的代用指标，常常会丢失曲线中蕴含的重要信息并且具有多解性。为了充分描述曲线的分布特征，在−1～11.3Φ（0.4～2 000 μm）范围内按间距0.16Φ划分75个粒级。同时，考虑到某些粒级变量的含量较低并且在不同样品之间存在显著差异，使用对数比转换的数据，可能会掩盖住数据的结构信息，本文选用经过标准化的原始数据[18-20]。采用欧式距离作为沉积物粒度频率分布曲线之间的相似性的度量，则样品x、y之间的相似度：
$$ {d_{ij}} = \sqrt {\sum\limits_{k = 1}^n {{{\left( {{x_{ik}} - {x_{jk}}} \right)}^2}} } $$ (5) 式中：xik，xjk分别表示沉积物样品i，j的第
$k$ 项变量含量。经过上述数据的前处理程序，AP聚类算法将85个剖面典型沉积相沉积物的粒度频率曲线划分为11簇（图4）。
11簇聚类结果显示，每一簇的粒度频率分布曲线分布范围比较集中，各粒级所对应的含量之间的差异很小，粒度频率分布曲线相似度较高。根据AP聚类算法的聚类结果采用silhouette函数绘制每一个点的轮廓值。轮廓值作为聚类分析有效性的常用内部指标，它被定义为：
$$ Sil\left( i \right) = \frac{{b\left( i \right) - a\left( i \right)}}{{\max \left\{ {a\left( i \right),b\left( i \right)} \right\}}} $$ (7) 式中，a（i）是第i个点与同簇的其他点之间的平均距离；b为一个向量，其元素是第i个点与属于不同类的各个点之间的平均距离。轮廓值的取值范围为[−1,1]，此值越大说明该点的分类越合理[21]。本文中11簇样品点的轮廓图如图5所示。
其中，各簇内样品点的轮廓值大部分大于0.3，仅有5个点的值小于零，最小为−0.08。考虑到划分聚类的数量以及各簇之间数据差异性相对较小，采用AP聚类算法对于沉积物粒度频率分布曲线的聚类划分是合理的，并且聚类效果显著。值得一提的是，AP聚类算法可以通过调整偏向参数能够增加或者减小聚类数量进而改变各簇内数据之间的差异大小，使得根据研究目的获得适合的聚类效果。本文选用数据点相似系数最大值作为偏向系数，以此保证各簇内数据差异很小，进而将动力条件差异不明显的沉积环境形成的粒度样品聚为一类。
4. 分析与讨论
4.1 基于AP聚类算法的沉积物沉积相识别
不同沉积类型的粒度分布曲线各粒级对应的百分含量不同，沉积类型对应的动力条件、运输介质等差异越大则粒度分布曲线分布区间、形态差异越明显[22]。当沉积物所受的搬运介质、搬运方式一定，并且介质动力大小稳定时，沉积物粒度组分的总体服从于某种分布的单因子随机事件的集合，在频率分布曲线和概率累计曲线形态上表现为单峰和一段式，否则呈多峰和多段式[23-24]。如前文所述，AP聚类算法将研究区沉积物的粒度频率分布曲线划分为11类簇，各簇内曲线形态、分布范围相近。但是，考虑到沉积物颗粒对于动力条件变化十分灵敏，同一沉积环境条件下粒度频率曲线会随着流体动力条件的波动在曲线形态、分布范围上产生相应的变化。因此，需要对曲线形态相似、分布范围发生轻微变化的曲线进一步归纳分析。
从11类簇粒度分布曲线可知，每一簇曲线与该簇的类代表性点之间的差异很小，利用类代表性曲线作为每簇曲线的典型样品曲线，并按曲线峰态、分布范围以及搬运方式将11簇曲线进一步划分为4组（图6）。第一组包含Ⅳ、Ⅵ、Ⅸ三簇对应的沉积粒度频率分布曲线，曲线形式表现为单峰、一段式（图4、6a），峰值分布于5～20 μm，分布范围较为集中，仅存在悬浮搬运方式，表明沉积物形成于动力条件较弱的静水环境；第二组包含Ⅱ、Ⅴ、Ⅷ、Ⅹ、Ⅺ五簇沉积物粒度频率分布曲线，曲线形态为双峰、两段式（图4、6b），主峰在20～200 μm，次峰在5～20 μm，跳跃和悬浮总体的搬运截点在4Φ左右，各总体的含量有所差异，表明沉积环境较第一组动力强度有所增加，并且水动力相对变化。第三组包括Ⅰ、Ⅶ类簇粒度分布曲线，为双峰、两段式（图4、6c），主峰在100～500 μm，次峰在10～100 μm，主要为跳跃总体存在，搬运截点变化较大，表明沉积物形成于强动力的沉积环境；第四组为Ⅳ类簇，为双峰、两段式（图4、6d），跳跃与悬浮总体含量相当，表明沉积物形成于动力条件强并且变化明显的沉积环境。
![]()
图 6 四组代表性样品频率曲线、概率曲线Figure 6. Frequency distribution and probability cumulative curves of four groups of representative samples根据上述分析，4类粒度频率分布曲线的组合可能代表了4种较为一致的沉积环境，将各组中代表性样品粒度频率分布曲线与稳定湖相（湖心相）、湖滨相、河流相等典型沉积相粒度频率分布曲线进行对比（图7）[23]。图中可以看出聚类结果得到的曲线与已知沉积环境的粒度曲线主次峰位置、曲线形态、分布范围较为一致。结果表明，第一组沉积物形成于稳定湖相，包括动力条件相似的湖心相、湖沼相；第二组沉积物形成于湖滨相或漫滩相；第三组沉积物形成于河流相；第四组为洪水沉积物。由此可见，AP聚类算法确定的4个组合分别代表了4种沉积环境。为了进一步验证上述分析的准确性与可行性，将上述结果与野外判别的4类85个沉积样品的沉积相进行对比发现，AP聚类算法确定的4个组合可以将野外识别的沉积相类别聚集到一起，两者具有较好的一致性。
![]()
图 7 各簇代表性与典型沉积相粒度频率分布曲线对比Figure 7. Comparison of frequency distribution curves of typical sedimentary facies and representative sediments in each cluster前人对大量已知沉积相的沉积物粒度参数进行定量描述，建立了多种沉积环境之间沉积物粒度参数判别公式。利用张平等根据典型稳定湖相与典型风成黄土沉积物粒度参数建立了稳定湖相与风成黄土粒度判别式（式8）[25]。本文利用此判别式对采集的所有样品进行沉积环境的判别分析。
$$\begin{split} F =& 20.363{\rm{Mz}} - 56.371{\rm{Sd}} - 67.922{\rm{SK}} +\\& 23.516{\rm{Ku}} - 55.626 \end{split}$$ (8) 式中，F表示判别值，若
$F>0$ 为稳定湖相沉积物；Mz、Sd、SK以及Ku分别是利用矩值法计算以$\Phi $ 为单位的平均值、分选系数、偏态以及峰态。其中第Ⅳ、Ⅵ、Ⅸ簇为正值，代表稳定湖相，判别结果与对比法得出的结果一致。
4.2 研究区不同沉积环境的沉积物粒度特征
基于粒度频率分布曲线的邻近传播聚类算法将全部85个粒度样品聚集成4类，并与已知沉积相粒度频率曲线进行对比，确定4类曲线对应的沉积环境分别为稳定湖相、湖滨相或漫滩相、河流相以及洪积相四种沉积相。各类沉积物之间在粒度频率分布曲线形态、分布范围以及组成成分上存在差异，通过定量描述这些差异能够为沉积环境、类型的区分提供依据[26]。粒度频率分布曲线利用粒度参数进行定量描述，粒度参数主要有粒度平均值、分选系数、偏态和峰态。为了更加充分揭示粒度参数对于沉积环境的指示意义，选用物理意义明确的Folk-Ward图解法进行计算（表2）[27]。表3给出了聚类后类代表性样品的粒度组成成分、粒度参数统计。
表 2 Folk&Ward 图解法公式Table 2. Formulas of Folk-Ward graphic methods粒度参数 计算公式 平均粒径（Mz） ${\rm{Mz}} = \frac{{{\Phi _{16}} + {\Phi _{50}} + {\Phi _{84}}}}{3}$ 分选系数（Sd） ${\rm{Sd}} = \frac{{{\Phi _{84}} - {\Phi _{16}}}}{4} + \frac{{{\Phi _{95}} - {\Phi _5}}}{{6.6}}$ 偏态（SK） ${\rm{SK}} = \frac{1}{2}\left[ {\frac{{{\Phi _{16}} + {\Phi _{84}} + 2{\Phi _{50}}}}{{{\Phi _{84}} - {\Phi _{16}}}}{\rm{ + }}\frac{{{\Phi _5} + {\Phi _{95}} + 2{\Phi _{50}}}}{{{\Phi _{95}} - {\Phi _5}}}} \right]$ 峰态（Ku） ${\rm{Ku}} = \frac{1}{{2.44}}\left[ {\frac{{{\Phi _{95}} - {\Phi _5}}}{{{\Phi _{75}} - {\Phi _{25}}}}} \right]$ 式中，粒径为$\Phi $值粒径，$\Phi = - \log _2^d$，$d$为毫米直径；${\Phi _x}$为x%累积含量的$\Phi $值粒径。 表 3 各簇代表性样品粒度组成成分及粒度参数Table 3. The size composition and size parameters of representative samples in each cluster簇类 类代表性样品 沉积相 粒度成分/% 平均粒径/Φ 分选系数 偏态 峰态 黏土 细粉砂 粗粉砂 细砂 中砂 粗砂 Ⅰ SITE01-01 河流相 1.27 3.72 6.95 18.33 45.34 24.38 1.462 1.155 0.344 1.584 Ⅱ SITE03-03 湖滨相 12.72 25.35 40.23 21.47 0.2 0.02 5.255 1.75 0.32 0.887 Ⅲ SITE04-06 洪积相 2.85 6.15 15.43 29.62 16.16 29.8 1.582 1.976 0.319 0.908 Ⅳ SITE12-02 湖沼相 19.67 41.22 34.45 4.33 0.14 0.13 6.227 1.722 0.096 1.066 Ⅴ SITE13-08 漫滩相 4.72 7.81 34.25 52.54 0.58 0.09 4.213 1.515 0.467 1.561 Ⅵ SITE13-09 湖沼相 40.67 47.21 10.66 1.46 0 0 7.537 1.301 0.019 1.158 Ⅶ SITE14-04 河流相 2.58 5.38 12.08 72.54 7.42 0 2.919 1.108 0.343 1.797 Ⅷ SITE17-04 漫滩相 6.68 10.46 55.59 27.18 0.08 0.01 4.577 1.343 0.363 1.485 Ⅸ SITE18-01 湖心相 20.74 55.49 22.38 1.4 0 0 6.819 1.329 0.175 1.123 Ⅹ SITE21-06 湖滨相 12.75 22.99 50.4 13.71 0.13 0.02 5.507 1.799 0.374 1.014 Ⅺ SITE18-02 湖滨相 15.1 33.35 43.32 7.98 0.11 0.12 5.841 1.745 0.209 0.981 从表3中可以看出聚类结果的11簇样品中，第Ⅱ、Ⅴ、Ⅷ、Ⅹ、Ⅺ簇样品之间粒度参数和组成成分相近，第Ⅳ、Ⅵ、Ⅸ簇样品之间粒度参数和组成成分相近，第Ⅰ、Ⅶ簇样品与第Ⅲ簇样品粒度分布形成于高能环境。4类典型沉积相的粒度组成与粒度参数之间存在显著差异。
第一类包括Ⅳ、Ⅵ、Ⅸ簇沉积物，粒度组分以细粉砂为主，黏土、细粉砂和粗粉砂三者含量之和占95%以上，平均粒径6.227～7.537Φ，在4类沉积环境中粒度成分最细，分选系数为1.301～1.722，分选差，除Ⅸ簇中部分样品呈正偏外，其余全为近对称分布，峰态变化于1.066～1.158，第Ⅳ 簇呈常峰态，其余全为尖峰态。
第二类包括Ⅱ、Ⅴ、Ⅷ、Ⅹ、Ⅺ簇沉积物，其中Ⅱ、Ⅷ、Ⅹ、Ⅺ簇沉积物以粗粉砂为主，占40%以上，黏土、细粉砂、细砂含量均大于10%，四者之和占98%以上，平均粒径为5.255～5.841Φ ，第Ⅴ簇沉积物以细砂（52%）为主，粗粉砂（34%）次之，平均粒径为4.213Φ，较Ⅱ、Ⅷ、Ⅹ、Ⅺ簇沉积物颗粒稍粗，表明强度略大于Ⅱ、Ⅷ、Ⅹ、Ⅺ簇沉积物对应的动力环境。该组分选系数为1.343～1.799，分选差，偏态为0.209～0.467，属于正偏、极正偏，峰态为0.887～1.561，第Ⅱ簇沉积物粒度分布曲线呈宽峰态，第Ⅴ、Ⅷ簇粒度分布曲线呈尖峰态或极尖峰态，其余两簇为常峰态。
第三类包括Ⅰ、Ⅶ簇沉积物，该组粒度组分较粗并且含量差异相对较大。其中，第Ⅰ簇沉积物以中砂（45%）为主，细砂（18%）、粗砂（22%）次之，平均粒径为1.462，样品分选差，呈极正偏，为尖峰态，代表强动力沉积环境；第Ⅶ簇沉积物中细砂含量占70%以上，平均粒径为2.919Φ，颗粒较第Ⅰ簇沉积物偏细，表明介质动力条件强并且稳定，相比于Ⅰ簇沉积物对应的沉积环境动力条件较低。
第四类中第Ⅲ簇沉积物以细砂为主，粗粉砂、中砂、粗砂含量均大于15%，平均粒径为1.582Φ，分选系数1.976，样品分选差，极正偏，宽峰态，表明沉积动力复杂多变且动力较强。
4.3 AP聚类算法在沉积环境识别中应用的可行性
根据上述方法对粒度分布曲线直接进行聚类分析，避免人为干扰聚类效果，有效地揭示研究区内的沉积环境差异是此方法的优点，也是能够利用此方法对粒度分布曲线进行研究区所存在的沉积环境类型反演的出发点。本文以白洋淀地区作为研究对象，尝试利用AP聚类算法建立粒度分布曲线与沉积环境之间的内在联系，通过识别结果与野外识别进行对比，发现所得结果较为一致，85个样品中有6个不一致。如下给出了沉积环境识别过程中存在差异的沉积物样品信息（表4）以及样品对应的粒度分布曲线（图8）。
表 4 沉积环境识别过程中存在差异的沉积物样品Table 4. Sediment samples with differences in the process of sedimentary environment identification样品编号 簇类 野外人为识别 粒度分布曲线对比识别 SITE01-03 Ⅱ 灰黑色粉砂质黏土，湖沼相沉积 湖滨相或漫滩相 SITE03-04 Ⅱ 灰黑色黏土，湖沼相沉积 湖滨相或漫滩相 SITE05-01 Ⅺ 灰黄色粉砂，河流相沉积 湖滨相或漫滩相 SITE13-02 Ⅺ 灰黄色粉砂，河流相沉积 湖滨相或漫滩相 SITE14-01 Ⅺ 灰黑色粉砂质黏土，湖相沉积 湖滨相或漫滩相 SITE14-02 Ⅱ 湖沼相沉积 湖滨相或漫滩相 ![]()
图 8 沉积环境识别过程中存在差异的样品粒度曲线Figure 8. Sediment samples with differences in sedimentary environment identification其中，SITE01-03、SITE14-01、SITE14-02三个样品的粒度曲线具有明显双峰，分布范围和曲线形态与典型湖滨相或漫滩相较为符合，野外观察结果存在识别错误（图8）。SITE03-04、SITE05-01、SITE13-02样品峰值粒径比典型湖沼相粒径较粗，但是没有明显次峰（图8）。对比其他簇曲线形态及分布范围，此类样品更为可能形成于沉积环境相对稳定，但动力强度较强的湖沼相沉积环境。AP聚类算法进行聚类分析时，由于样品相对典型湖沼相样品粒径较粗，被划分到湖滨相沉积物一类，出现了偏差。因此，对于85个样品中有3个样品沉积环境识别出现误差，错误率为3.53%。
鉴于聚类方法是尽可能将相似的数据点聚集成一类，沉积物的形成受物源、介质等多个因素影响，同一沉积环境条件下，粒度分布曲线会在一定范围内发生波动，因此采用邻近传播聚类分析时需要合理选取偏向系数控制各簇之间数据的差异。本文为了充分挖掘粒度数据中蕴含的沉积环境类型，选取了相似系数最大值作为偏向系数，尽可能减小数据之间差异，导致出现同一沉积环境条件下形成的粒度曲线可能被划分到不同簇中的现象。
5. 结论
（1）AP聚类算法具有不需要预设聚类数量，能够准确地确定研究区沉积环境中存在的动力类型，将曲线形态、分布范围相近的沉积物聚成一类。每一簇的曲线形态和分布范围都能反演出研究区对应的沉积环境存在的动力类型，并且建立沉积物粒度特征与沉积动力学之间的内在联系，为沉积物类型划分提供依据。
（2）根据粒度概率累积曲线以及分布曲线特征，将11簇曲线中形态、分布范围相近的曲线进一步合并，得到4类曲线，并与已知沉积环境样品曲线进行对比，确定出白洋淀地区剖面地层中出露着湖沼相、湖相、河流相以及洪积相四种沉积相。
（3）白洋淀地区存在的4类典型沉积环境中湖心相、湖沼相的动力环境相近，以细粉砂为主，分选差，动力强度最低；湖滨相以粗粉砂为主，双峰明显，分选差。河流相颗粒较粗，以细、中砂为主，沉积物粒度特征反映漫滩相沉积环境与湖滨相相似；洪积相具有明显双峰，含量相近，分选差，水动力条件动荡。
由此可见，基于粒度频率分布曲线的邻近聚类传播算法可以作为挖掘研究区存在的沉积类型的新方法，为沉积环境动力条件反演、分区等提供潜在的新手段。
-
![]()
图 1 研究区地质概况、马山柱状节理火山岩露头与镜下特征
a. 研究区地质简图（据文献[13]），b. 马山石林剖面，c. 粗面英安岩中的晶洞构造，d. 马山粗面英安岩镜下显微特征，左侧为单偏光（-），右侧为正交光（+）。
Figure 1. Geological map of Mashan region with pictures of outcrops and photomicrographs of the volcanic rocks with columnar joints
a. sketched geological map of study area （According to reference [13]）, b. section of columnar joints developed in Mashan, c. geode structure of the trachydacite, d. petrographic characteristics of Mashan trachydacite under microscope, with single polarized light （-） on the left and crossed light （+） on the right.
![]()
图 2 马山粗面英安岩锆石特征与测年结果
a. 马山粗面英安岩的锆石阴极发光图像及对应的Pb206/U238年龄与εHf（t）值，b. 加权平均值，c. U-Pb年龄谐和图。
Figure 2. Zircon characteristics and dating results of the trachyandesite in Mashan
a. cathodoluminescence （CL） images, yellow numbers=U-Pb ages and red number= εHf（t） values; b. weighted mean age, c. concordia diagrams of zircons from the Mashan trachydacite.
![]()
![]()
图 4 马山粗面英安岩中锆石的εHf（t）值与对应的二阶段模式年龄
Figure 4. εHf（t） values and two stage Hf model ages of zircons from trachydacite in Mashan
![]()
表 1 马山粗面英安岩锆石U-Pb测试结果
Table 1 Zircon U-Pb dating results of trachydacite at Mashan
点 号 Th /10−6 U/10−6 Th/U 比值 年龄/Ma 谐和度*/% 207Pb/206Pb ± 1σ 207Pb/235U ± 1σ 206Pb/238U ± 1σ 207Pb/206Pb ± 1σ 207Pb/235U ± 1σ 206Pb/238U ± 1σ MS-01 3.2 370.2 0.01 0.04995 0.00297 0.13378 0.00741 0.01984 0.00039 195 134 127 7 127 2 99 MS-02 165.0 235.1 0.70 0.10304 0.00387 2.14901 0.07840 0.15008 0.00276 1679 69 1165 25 901 15 74 MS-03 0.0 0.45 0.01 0.03465 0.01682 4.37117 2.07928 0.03412 0.01459 -- -- 1707 393 216 91 56 MS-04 1327.3 935.3 1.42 0.04782 0.00197 0.11810 0.00484 0.01785 0.00027 100 87 113 4 114 2 99 MS-05 376.9 294.5 1.28 0.04951 0.00329 0.11991 0.00749 0.01787 0.00040 172 156 115 7 114 3 99 MS-06 100.5 115.3 0.87 0.05643 0.00525 0.13687 0.01219 0.01767 0.00053 478 205 130 11 113 3 85 MS-07 899.6 806.2 1.12 0.05167 0.00246 0.12889 0.00600 0.01787 0.00037 333 105 123 5 114 2 92 MS-08 96.9 148.0 0.65 0.11618 0.00370 5.07929 0.16091 0.31262 0.00570 1898 57 1833 27 1754 28 95 MS-09 1121.3 809.3 1.39 0.05088 0.00258 0.12693 0.00628 0.01795 0.00036 235 114 121 6 115 2 94 MS-10 80.6 120.9 0.67 0.08139 0.00443 0.71790 0.04065 0.06321 0.00167 1231 140 549 24 395 10 67 MS-11 0.0 6.2 0.00 0.09050 0.02278 0.45048 0.12001 0.03243 0.00361 1436 496 378 84 206 23 41 MS-12 510.8 1110.0 0.46 0.05150 0.00206 0.12502 0.00506 0.01727 0.00030 265 93 120 5 110 2 91 MS-13 167.7 2466.9 0.07 0.16904 0.00406 9.84273 0.23879 0.41252 0.00633 2548 40 2420 22 2226 29 91 MS-14 184.5 166.1 1.11 0.05229 0.00471 0.12991 0.01039 0.01795 0.00046 298 210 124 9 115 3 92 MS-15 69.1 246.8 0.28 0.05578 0.00257 0.26734 0.01174 0.03436 0.00062 443 102 241 9 218 4 90 MS-16 109.4 315.4 0.35 0.05255 0.00237 0.24168 0.01064 0.03271 0.00060 309 99 220 9 208 4 94 MS-17 937.5 625.8 1.50 0.05214 0.00239 0.13070 0.00511 0.01794 0.00034 300 106 125 5 115 2 91 MS-18 396.1 555.5 0.71 0.05303 0.00222 0.13004 0.00486 0.01749 0.00032 332 94 124 4 112 2 90 MS-19 52.9 1325.4 0.04 0.05007 0.00146 0.22608 0.00622 0.03209 0.00044 198 69 207 5 204 3 98 MS-20 614.1 290.2 2.12 0.04650 0.00531 0.11427 0.01412 0.01715 0.00048 33 243 110 13 110 3 99 MS-21 3467.1 1542.5 2.25 0.04698 0.00156 0.11737 0.00414 0.01766 0.00029 56 72 113 4 113 2 99 MS-22 1321.0 881.1 1.50 0.04750 0.00189 0.12030 0.00482 0.01797 0.00034 76 89 115 4 115 2 99 MS-23 412.3 496.5 0.83 0.04682 0.00251 0.11249 0.00555 0.01728 0.00040 39 126 108 5 110 3 97 注：*谐和度=100${\rm{ \times \{ 1 - [}}{{\rm{(}}^{{\rm{206}}}}{\rm{Pb}}{{\rm{/}}^{{\rm{238}}}}{\rm{U}}\;{{\rm{ - }}^{{\rm{207}}}}{\rm{Pb}}{{\rm{/}}^{{\rm{235}}}}{\rm{U}}\;{\rm{)/(}}{{\rm{(}}^{{\rm{206}}}}{\rm{Pb}}{{\rm{/}}^{{\rm{238}}}}{\rm{U}}\;{\rm{}}{{\rm{ + }}^{{\rm{207}}}}{\rm{Pb}}{{\rm{/}}^{{\rm{235}}}}{\rm{U}}\;{\rm{)/2)]\} }}$ 
下载: 导出CSV
表 2 马山粗面英安岩锆石原位Lu-Hf测试结果
Table 2 Zircon in-situ Lu-Hf results of trachydacite at Mashan
点 号 176Yb/177Hf ±2σ 176Lu/177Hf ±2σ 176Hf/177Hf ±2σ t (Ma)① T C ② DM ±2σ εHf(t)③ ±2σ MS-01 0.010770 0.000113 0.000468 0.000004 0.282340 0.000014 127 1971 32 −12.5 0.51 MS-04 0.092622 0.002184 0.003462 0.000074 0.282242 0.000021 114 2209 47 −16.5 0.76 MS-05 0.080054 0.001520 0.003027 0.000039 0.282239 0.000021 114 2213 46 −16.6 0.74 MS-07 0.145124 0.004230 0.005205 0.000141 0.282221 0.000020 114 2263 44 −17.4 0.71 MS-08 0.004921 0.000209 0.000125 0.000005 0.281265 0.000015 1898 3263 33 −11.1 0.54 MS-09 0.105545 0.000874 0.003823 0.000037 0.282214 0.000017 115 2272 38 −17.5 0.61 MS-12 0.133218 0.001571 0.004997 0.000055 0.282232 0.000033 110 2239 72 −17.0 1.17 MS-13 0.022171 0.000979 0.000714 0.000025 0.281177 0.000027 2548 3087 58 −0.4 0.96 MS-14 0.064095 0.001647 0.002542 0.000062 0.282191 0.000029 115 2317 63 −18.2 1.01 MS-15 0.022432 0.000587 0.000875 0.000019 0.282176 0.000016 218 2287 36 −16.4 0.58 MS-16 0.013703 0.000135 0.000627 0.000005 0.282294 0.000012 208 2030 28 −12.5 0.44 MS-17 0.110210 0.000395 0.004199 0.000016 0.282159 0.000017 115 2395 38 −19.5 0.61 MS-19 0.024116 0.000554 0.001112 0.000024 0.282102 0.000012 204 2458 27 −19.4 0.44 MS-20 0.082807 0.001461 0.003316 0.000056 0.282225 0.000025 110 2248 55 −17.2 0.88 MS-21 0.150046 0.004172 0.005299 0.000146 0.282263 0.000023 113 2172 51 −15.9 0.82 MS-22 0.132288 0.002429 0.005110 0.000093 0.282254 0.000022 115 2191 48 −16.2 0.77 MS-23 0.074151 0.004070 0.002805 0.000147 0.282187 0.000018 110 2330 39 −18.5 0.64 注：①为锆石表观年龄；②为基于DM的二阶段模式年龄计算使用平均大陆地壳176Lu/177Hf值0.015。表中数据的误差为2σ（标准误差），不确定度为最后一位数字；③为εHf(t)值的计算使用锆石的表观年龄，该年龄为锆石结晶年龄的最小估计值，λ176Lu=1.867×10−11/a[17]，CHUR现今的176Lu/177Hf及176Hf/177Hf值分别为0.0332和0.282772[18]，DM现今的176Lu/177Hf及176Hf/177Hf值分别为0.0384和0.28325[19]。
下载: 导出CSV
表 3 马山粗面英安岩全岩主量元素分析结果
Table 3 Whole-rock major element compositions of trachydacite at Mashan
% 样品号 M-01 M-02 M-03 M-04 平均值 SiO2 65.44 67.26 67.06 67.18 66.74 TiO2 0.46 0.45 0.44 0.40 0.44 Al2O3 14.58 13.37 13.43 13.38 13.69 Fe2O3 3.61 3.89 3.65 3.62 3.69 MnO 0.08 0.10 0.08 0.09 0.09 MgO 2.10 1.95 2.14 2.36 2.14 CaO 3.40 2.67 2.76 2.58 2.85 Na2O 3.77 4.48 4.71 4.41 4.34 K2O 3.83 3.76 3.58 3.98 3.79 P2O5 0.18 0.19 0.21 0.17 0.19 L.O.I 2.12 1.46 1.86 1.65 1.77
下载: 导出CSV
表 4 马山粗面英安岩微量元素分析结果
Table 4 Trace element compositions of trachydacite at Mashan
10−6 样品号 M-01 M-02 M-03 M-04 平均值 Li 37.5 41.4 36.4 38.6 38.5 Be 1.75 1.66 1.70 1.68 1.70 Sc 18.30 7.80 7.91 8.10 10.50 Ti 2823 2700 2640 2400 2641 V 53.1 55.1 56.1 54.4 54.7 Cr 68.5 63.5 65.6 81.6 69.8 Co 10.80 11.00 11.00 11.40 11.10 Ni 32.1 29.5 29.6 39.2 32.6 Cu 4.76 45.80 28.60 28.20 26.80 Zn 48.5 72.2 68.0 67.2 64.0 Ga 17.9 16.7 17.2 15.9 16.9 Rb 91.6 85.0 87.1 99.2 90.7 Sr 402 341 422 489 414 Y 13.4 13.6 13.9 13.0 13.5 Zr 130 141 134 126 133 Nb 9.5 11.0 10.8 10.6 10.5 Mo 0.61 0.88 0.58 0.66 0.68 Cd 0.036 0.061 0.066 0.038 0.051 Cs 0.88 0.62 0.65 1.00 0.79 Ba 2381 1461 1939 1727 1877 Hf 3.73 4.35 4.40 4.00 4.12 Ta 0.73 1.20 1.00 1.15 1.02 Tl 0.38 0.38 0.38 0.44 0.40 Pb 30.3 54.6 39.2 43.6 41.9 Th 6.46 5.56 5.38 5.26 5.67 U 1.69 1.50 1.54 1.54 1.57
下载: 导出CSV
表 5 马山粗面英安岩稀土元素分析结果
Table 5 REE compositions of trachydacite at Mashan
10−6 样品号 M-01 M-02 M-03 M-04 平均值 La 29.6 29.9 30.8 27.8 29.5 Ce 52.4 53.1 54.6 49.5 52.4 Pr 6.34 6.18 6.34 5.77 6.16 Nd 22.9 23.3 24.0 21.8 23.0 Sm 3.96 3.86 3.92 3.64 3.85 Eu 0.94 1.08 1.10 1.00 1.03 Gd 3.53 3.28 3.38 3.10 3.32 Tb 0.50 0.47 0.48 0.46 0.48 Dy 2.50 2.53 2.60 2.46 2.52 Ho 0.50 0.51 0.52 0.48 0.50 Er 1.39 1.42 1.44 1.38 1.41 Tm 0.24 0.22 0.22 0.22 0.23 Yb 1.45 1.46 1.44 1.40 1.44 Lu 0.23 0.23 0.23 0.22 0.23 Eu/Eu* 0.77 0.93 0.92 0.91 0.88 LaN/YbN 13.89 13.94 14.56 13.51 13.98 ΣREE 126.48 127.54 131.07 119.23 126.08
下载: 导出CSV
-
[1] Cope T. Phanerozoic magmatic tempos of North China [J]. Earth and Planetary Science Letters, 2017, 468: 1-10. doi: 10.1016/j.jpgl.2017.03.022
[2] 任纪舜, 陈廷愚, 牛宝贵, 等. 中国东部及邻区大陆岩石圈的构造演化与成矿[M]. 北京: 科学出版社, 1990: 1-205. REN Jishun, CHEN Tingyu, NIU Baogui, et al. Tectonic Evolution of the Continental Lithosphere and Metallogeny in Eastern China and Adjacent Areas[M]. Beijing: Science Press, 1990: 1-205.
[3] 邓晋福, 刘厚祥, 赵海玲, 等. 燕辽地区燕山期火成岩与造山模型[J]. 现代地质, 1996, 10(2):137-148 DENG Jinfu, LIU Houxiang, ZHAO Hailing, et al. Yanshanian igneous rocks and orogeny model in Yanshan-Liaoning area [J]. Geoscience, 1996, 10(2): 137-148.
[4] 董树文, 吴锡浩, 吴珍汉, 等. 论东亚大陆的构造翘变: 燕山运动的全球意义[J]. 地质论评, 2000, 46(1):8-13 doi: 10.3321/j.issn:0371-5736.2000.01.002 DONG Shuwen, WU Xihao, WU Zhenhan, et al. On tectonic seesawing of the East Asia Continent: Global implication of the Yanshanian Movement [J]. Geological Review, 2000, 46(1): 8-13. doi: 10.3321/j.issn:0371-5736.2000.01.002
[5] 董树文, 张岳桥, 李海龙, 等. “燕山运动”与东亚大陆晚中生代多板块汇聚构造: 纪念“燕山运动”90周年[J]. 中国科学: 地球科学, 2018, 61(6):1888-1909 DONG Shuwen, ZHANG Yueqiao, LI Hailong, et al. The Yanshan orogeny and late Mesozoic multi-plate convergence in East Asia: Commemorating 90th years of the “Yanshan Orogeny” [J]. Science China Earth Sciences, 2018, 61(6): 1888-1909.
[6] 张旗, 钱青, 王二七, 等. 燕山中晚期的中国东部高原: 埃达克岩的启示[J]. 地质科学, 2001, 36(2):248-255 doi: 10.3321/j.issn:0563-5020.2001.02.014 ZHANG Qi, QIAN Qing, WANG Erqi, et al. An East China Plateau in Mid-Late Yanshanian Period: Implication from adakites [J]. Chinese Journal of Geology, 2001, 36(2): 248-255. doi: 10.3321/j.issn:0563-5020.2001.02.014
[7] 张旗, 许继峰, 王焰, 等. 埃达克岩的多样性[J]. 地质通报, 2004, 23(9-10):959-965 ZHANG Qi, XU Jifeng, WANG Yan, et al. Diversity of adakite [J]. Geological Bulletin of China, 2004, 23(9-10): 959-965.
[8] 张旗, 焦守涛. 埃达克岩来自高压背景: 一个科学的、可靠的、有预见性的科学发现[J]. 岩石学报, 2020, 36(6):1675-1683 doi: 10.18654/1000-0569/2020.06.02 ZHANG Qi, JIAO Shoutao. Adakite comes from a high-pressure background: A scientific, reliable, predictable scientific discovery [J]. Acta Petrologica Sinica, 2020, 36(6): 1675-1683. doi: 10.18654/1000-0569/2020.06.02
[9] 汪洋, 程素华. “C型埃达克岩”: 一个基于误解的概念?[J]. 矿物岩石地球化学通报, 2010, 29(3):284-292, 298 doi: 10.3969/j.issn.1007-2802.2010.03.012 WANG Yang, CHENG Suhua. The “C-type Adakite”: A concept based on misunderstanding? [J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2010, 29(3): 284-292, 298. doi: 10.3969/j.issn.1007-2802.2010.03.012
[10] 汪洋. 涉及埃达克岩实验的若干基本概念与问题辨析[J]. 岩石矿物学杂志, 2016, 35(1):162-176 doi: 10.3969/j.issn.1000-6524.2016.01.012 WANG Yang. A discussion on some basic conceptions and problems related to the experimental study of adakite [J]. Acta Petrologica Et Mineralogica, 2016, 35(1): 162-176. doi: 10.3969/j.issn.1000-6524.2016.01.012
[11] Martin H, Smithies R H, Rapp R, et al. D. An overview of adakite, tonalite–trondhjemite–granodiorite (TTG), and sanukitoid: relationships and some implications for crustal evolution [J]. Lithos, 2005, 79(1-2): 1-24. doi: 10.1016/j.lithos.2004.04.048
[12] 王强, 郝露露, 张修政, 等. 汇聚板块边缘的埃达克质岩: 成分和成因[J]. 中国科学: 地球科学, 2020, 63(12):1992-2016 doi: 10.1007/s11430-020-9678-y WANG Qiang, HAO Lulu, ZHANG Xiuzheng, et al. Adakitic rocks at convergent plate boundaries: Compositions and petrogenesis [J]. Science China Earth Sciences, 2020, 63(12): 1992-2016. doi: 10.1007/s11430-020-9678-y
[13] 山东省地质调查院. 1: 25万青岛市幅区域地质调查报告与地质图[R]. 2004. Shandong Institute of Geological Survey. 1: 250000 Qingdao Regional Geological Survey Report and Geological Map[R]. 2004.
[14] 韩宗珠, 衣伟虹, 余芝华, 等. 青岛马山地区安山玢岩地球化学特征及成因机理研究[J]. 中国海洋大学学报, 2011, 41(S1):305-311 HAN Zongzhu, YI Weihong, YU Zhihua, et al. Geochemistry and petrogenesis of andesite-prophyrite in Mashan area, Qingdao [J]. Periodical of Ocean University of China, 2011, 41(S1): 305-311.
[15] 何登洋, 邱昆峰, 于皓丞, 等. 华北克拉通胶莱盆地马山地区早白垩世粗面英安岩岩石成因[J]. 岩石学报, 2020, 36(12):3705-3720 doi: 10.18654/1000-0569/2020.12.09 HE Dengyang, QIU Kunfeng, YU Haocheng, et al. Petrogenesis of the Early Cretaceous trachy-dacite from Mashan in the Jiaolai Basin, North China Craton [J]. Acta Petrologica Sinica, 2020, 36(12): 3705-3720. doi: 10.18654/1000-0569/2020.12.09
[16] Pearce N J G, Perkins W T, Westgate J A, et al. A compilation of new and published major and trace element data for NIST SRM 610 and NIST SRM 612 glass reference materials [J]. Geostandards Newsletter, 1997, 21(1): 115-144. doi: 10.1111/j.1751-908X.1997.tb00538.x
[17] Söderlund U, Patchett P J, Vervoort J D, et al. The 176Lu decay constant determined by Lu-Hf and U-Pb isotope systematics of Precambrian mafic intrusions [J]. Earth and Planetary Science Letters, 2004, 219(3-4): 311-324. doi: 10.1016/S0012-821X(04)00012-3
[18] Bichert-Toft J, Albarède F. The Lu-Hf isotope geochemistry of chondrites and the evolution of the mantle-crust system [J]. Earth and Planetary Science Letters, 1997, 148(1-2): 243-258. doi: 10.1016/S0012-821X(97)00040-X
[19] Griffin W L, Pearson N J, Belousova E, et al. The Hf isotope composition of cratonic mantle: LAM-MC-ICPMS analysis of zircon megacrysts in kimberlites [J]. Geochimica et Cosmochimica Acta, 2000, 64(1): 133-147. doi: 10.1016/S0016-7037(99)00343-9
[20] 吴元保, 郑永飞. 锆石成因矿物学研究及其对U-Pb年龄解释的制约[J]. 科学通报, 2004, 49(15):1554-1569 doi: 10.1007/BF03184122 WU Yuanbao, ZHENG Yongfei. Genesis of zircon and its constraints on interpretation of U-Pb age [J]. Chinese Science Bulletin, 2004, 49(15): 1554-1569. doi: 10.1007/BF03184122
[21] McDonough W F, Sun S S. The composition of the Earth [J]. Chemical Geology, 1995, 120(3-4): 223-253. doi: 10.1016/0009-2541(94)00140-4
[22] Le Bas M J, Le Maitre R W, Streckeisen A, et al. A chemical classification of volcanic rocks based on the total alkali-silica diagram [J]. Journal of Petrology, 1986, 27(3): 745-750. doi: 10.1093/petrology/27.3.745
[23] Kay R W. Aleutian magnesian andesites: Melts from subducted Pacific ocean crust [J]. Journal of Volcanology and Geothermal Research, 1978, 4(1-2): 117-132. doi: 10.1016/0377-0273(78)90032-X
[24] Defant M J, Drummond M S. Derivation of some modern arc magmas by melting of young subducted lithosphere [J]. Nature, 1990, 347(6294): 662-665. doi: 10.1038/347662a0
[25] Defant M J, Xu J F, Kepezhinskas P, et al. Adakites: some variations on a theme [J]. Acta Petrologica Sinica, 2002, 18(2): 129-142.
[26] 王焰, 张旗, 钱青. 埃达克岩(adakite)的地球化学特征及其构造意义[J]. 地质科学, 2000, 35(2):251-256 doi: 10.3321/j.issn:0563-5020.2000.02.016 WANG Yan, ZHANG Qi, QIAN Qing. Adakite: geochemical characteristics and tectonic significances [J]. Scientia Geologica Sinica, 2000, 35(2): 251-256. doi: 10.3321/j.issn:0563-5020.2000.02.016
[27] 王强, 许继锋, 赵振华. 一种新的火成岩: 埃达克岩的研究综述[J]. 地球科学进展, 2001, 16(2):201-208 doi: 10.3321/j.issn:1001-8166.2001.02.010 WANG Qiang, XU Jifen, ZHAO Zhenghua. The summary and comment on research on a new kind of igneous rock-Adakite [J]. Advance in Earth Sciences, 2001, 16(2): 201-208. doi: 10.3321/j.issn:1001-8166.2001.02.010
[28] Rapp R P, Xiao L, Shim Izu N. Experimental constraints on the origin of potassium-rich adakites in eastern China [J]. Acta Petrologica Sinica, 2002, 18(3): 293-302.
[29] 罗照华, 柯珊, 谌宏伟. 埃达克岩的特征、成因及构造意义[J]. 地球科学进展, 2002, 21(7):436-440 LUO Zhaohua, KE Shan, CHEN Hongwei. Characteristics, petrogenesis and tectonic implications of adakite [J]. Geological Bulletin of China, 2002, 21(7): 436-440.
[30] 翟明国. 埃达克岩和大陆下地壳重熔的花岗岩类[J]. 岩石学报, 2004, 20(2):193-194 doi: 10.3969/j.issn.1000-0569.2004.02.001 ZHAI Mingguo. Adakite and related granitoids from partial melting of continental lower crust [J]. Acta Petrologica Sinica, 2004, 20(2): 193-194. doi: 10.3969/j.issn.1000-0569.2004.02.001
[31] 张旗, 王元龙, 金惟俊, 等. 晚中生代的中国东部高原: 证据、问题和启示[J]. 地质通报, 2008, 27(9):1404-1430 doi: 10.3969/j.issn.1671-2552.2008.09.004 ZHANG Qi, WANG Yuanlong, JI Weijun, et al. Eastern China Plateau during the Late Mesozoic: evidence, problems and implication [J]. Geological Bulletin of China, 2008, 27(9): 1404-1430. doi: 10.3969/j.issn.1671-2552.2008.09.004
[32] Martin H. Effect of steeper Archean geothermal gradient on geochemistry of subduction-zone magmas [J]. Geology, 1986, 14(9): 753-756. doi: 10.1130/0091-7613(1986)14<753:EOSAGG>2.0.CO;2
[33] 邓晋福, 苏尚国, 刘翠, 等. 关于华北克拉通燕山期岩石圈减薄的机制与过程的讨论: 是拆沉, 还是热侵蚀和化学交代?[J]. 地学前缘, 2006, 13(2):105-119 doi: 10.3321/j.issn:1005-2321.2006.02.009 DENG Jinfu, SU Shangguo, LIU Cui, et al. Discussion on the lithospheric thinning of the North China craton: delamination? Or thermal erosion and chemical metasomatism? [J]. Earth Science Frontiers, 2006, 13(2): 105-119. doi: 10.3321/j.issn:1005-2321.2006.02.009
[34] Best M G. Igneous and Metamorphic Petrology[M]. 2nd ed. Malden, Mass: Blackwell Science, 2003: 340-341.
[35] Taylor S R, McClennan S M. The Continental Crust: Its Composition and Evolution[M]. Oxford: Blackwell Scientific Publication, 1985: 312.
[36] 张旗, 王焰, 熊小林, 等. 埃达克岩和花岗岩: 挑战与机遇[M]. 北京: 中国大地出版社, 2008: 1-344. ZHANG Qi, WANG Yan, XIONG Xiaolin, et al. Adakite and Granite: Challenges and Opportunities[M]. Beijing: China Land Press, 2008: 1-344.
[37] 赵越, 徐刚, 张拴宏, 等. 燕山运动与东亚构造体制的转变[J]. 地学前缘, 2004, 11(3):319-328 doi: 10.3321/j.issn:1005-2321.2004.03.030 ZHAO Yue, XU Gang, ZHANG Shuanhong, et al. Yanshanian movement and conversion of tectonic regimes in East Asia [J]. Earth Science Frontiers, 2004, 11(3): 319-328. doi: 10.3321/j.issn:1005-2321.2004.03.030
[38] 孟庆任. 中国东部燕山期沉积盆地演化与资源效应[J]. 矿物岩石地球化学通报, 2017, 36(4):567-569 doi: 10.3969/j.issn.1007-2802.2017.04.007 MENG Qingren. Development of sedimentary basins in Eastern China during the Yanshanian period [J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2017, 36(4): 567-569. doi: 10.3969/j.issn.1007-2802.2017.04.007
[39] 邵济安, 刘福田, 陈辉, 等. 大兴安岭—燕山晚中生代岩浆活动与俯冲作用关系[J]. 地质学报, 2001, 75(1):56-63 doi: 10.3321/j.issn:0001-5717.2001.01.006 SHAO Ji’an, LIU Futian, CHEN Hui, et al. Relationship between Mesozoic Magmatism and Subduction in Da Hinggan-Yanshan Area [J]. Acta Geologica Sinica, 2001, 75(1): 56-63. doi: 10.3321/j.issn:0001-5717.2001.01.006
[40] 胡受奚, 赵乙英, 胡志宏, 等. 中国东部中—新生代活动大陆边缘构造: 岩浆作用演化和发展[J]. 岩石学报, 1994, 10(4):370-381 doi: 10.3321/j.issn:1000-0569.1994.04.003 HU Shouxi, ZHAO Yiying, HU Zhihong, et al. Evolution and development of tectonics and magmatism at the active continental margin of the East China (E106°) during Mesozoic and Cenozoic [J]. Acta Petrologica Sinica, 1994, 10(4): 370-381. doi: 10.3321/j.issn:1000-0569.1994.04.003
[41] 王德滋, 任启江, 邱检生, 等. 中国东部橄榄安粗岩省的火山岩特征及其成矿作用[J]. 地质学报, 1996, 70(1):23-34 WANG Dezi, REN Qijiang, QIU Jiansheng, et al. Characteristics of volcanic rocks in the Shoshonite province, Eastern China, and their metallogenesis [J]. Acta Geologica Sinica, 1996, 70(1): 23-34.
[42] 万天丰, 赵庆乐. 中国东部构造-岩浆作用的成因[J]. 中国科学: 地球科学, 2012, 55(3):347-354 doi: 10.1007/s11430-011-4361-4 WAN Tianfeng, ZHAO Qingle. The genesis of tectono-magmatism in eastern China [J]. Science China Earth Sciences, 2012, 55(3): 347-354. doi: 10.1007/s11430-011-4361-4
[43] 张旗. 中国东部中生代岩浆活动与太平洋板块向西俯冲有关吗?[J]. 岩石矿物学杂志, 2013, 32(1):113-128 doi: 10.3969/j.issn.1000-6524.2013.01.010 ZHANG Qi. Is the Mesozoic magmatism in eastern China related to westward subduction of the Pacific plate? [J]. Acta Petrologica ET Mineralogica, 2013, 32(1): 113-128. doi: 10.3969/j.issn.1000-6524.2013.01.010
[44] Zhu X Q, Zhu W B, Ge R F, et al. Late Paleozoic provenance shift in the south-central North China Craton: Implications for tectonic evolution and crustal growth [J]. Gondwana Research, 2014, 25(1): 383-400. doi: 10.1016/j.gr.2013.04.009
[45] Li S Z, Zhao S J, Liu X, et al. Closure of the Proto-Tethys Ocean and Early Paleozoic amalgamation of microcontinental blocks in East Asia [J]. Earth-Science Reviews, 2018, 186: 37-75. doi: 10.1016/j.earscirev.2017.01.011
[46] 张岳桥, 李金良, 张田, 等. 胶莱盆地及其邻区白垩纪—古新世沉积构造演化历史及其区域动力学意义[J]. 地质学报, 2008, 82(9):1229-1257 doi: 10.3321/j.issn:0001-5717.2008.09.007 ZHANG Yueqiao, LI Jinliang, ZHANG Tian, et al. Cretaceous to Paleocene Tectono-sedimentary evolution of the Jiaolai Basin and the contiguous areas of the Shandong Peninsula (North China) and Its geodynamic implications [J]. Acta Geologica Sinica, 2008, 82(9): 1229-1257. doi: 10.3321/j.issn:0001-5717.2008.09.007
[47] 张田, 张岳桥. 胶北隆起晚中生代构造-岩浆演化历史[J]. 地质学报, 2008, 82(9):1210-1228 doi: 10.3321/j.issn:0001-5717.2008.09.006 ZHANG Tian, ZHANG Yueqiao. Late Mesozoic tectono-magmatic evolution history of the Jiaobei Uplift, Shandong Peninsula [J]. Acta Geologica Sinica, 2008, 82(9): 1210-1228. doi: 10.3321/j.issn:0001-5717.2008.09.006
[48] Guo Z X, Yang Y T, Zhao X Z, et al. Early Cretaceous tectonostratigraphic evolution of the Erlian Basin, NE China: A record of Late Mesozoic intraplate deformation in East Asia [J]. Marine and Petroleum Geology, 2019, 110: 539-564. doi: 10.1016/j.marpetgeo.2019.07.043
[49] 吕洪波, 王俊, 张海春. 山东灵山岛晚中生代滑塌沉积层的发现及区域构造意义初探[J]. 地质学报, 2011, 85(6):938-946 LÜ Hongbo, WANG Jun, ZHANG Haichun. Discovery of the Late Mesozoic slump beds in Lingshan Island, Shandong, and a pilot research on the regional tectonics [J]. Acta Geologica Sinica, 2011, 85(6): 938-946.
[50] 孙天柱, 武斌. 临沂方城盆地中生代青山群火山岩地球化学特征及地质意义[J]. 山东国土资源, 2020, 36(6):14-22 SUN Tianzhu, WU Bin. Geochemical characteristics and geological significance of volcanic rocks in Mesozoic Qingshan Group in Fangcheng Basin of Linyi City [J]. Shandong Land and Resources, 2020, 36(6): 14-22.
[51] 王宇林, 郭强, 李喜海, 等. 辽西喀左盆地九佛堂组的划分及相关问题[J]. 地层学杂志, 2009, 33(2):147-154 doi: 10.3969/j.issn.0253-4959.2009.02.004 WANG Yulin, GUO Qiang, LI Xihai, et al. Subdivision of the Jiufotang Formation and related issues in the Kazuo Basin [J]. Journal of Stratigraphy, 2009, 33(2): 147-154. doi: 10.3969/j.issn.0253-4959.2009.02.004
[52] 刘为付, 刘双龙, 孙立新, 等. 包日温都构造带九佛堂组高分辨率层序结构及储层分析[J]. 地层学杂志, 2015, 39(3):310-318 LIU Weifu, LIU Shuanglong, SUN Lixin, et al. High resolution sequence stratigraphy configurations and reservoir analysis of the Jiufotang Formation in the Baoriwendu tectonic belt [J]. Journal of Stratigraphy, 2015, 39(3): 310-318.
[53] 姜同海. 胶莱盆地构造特征及演化[D]. 山东科技大学硕士学位论文, 2005: 1-45. JIANG Tonghai. The structural characteristics and evolution of Jiaolai Basin[D]. Master Dissertation of Shandong University of Science and Technology, 2005: 1-45.
[54] Wu L, Monié P, Wang F, et al. Multi-phase cooling of Early Cretaceous granites on the Jiaodong Peninsula, East China: Evidence from 40Ar/39Ar and (U-Th)/He thermochronology [J]. Journal of Asian Earth Sciences, 2018, 160: 334-347. doi: 10.1016/j.jseaes.2017.11.014
[55] Zhang L M, Wang C S, Cao K, et al. High elevation of Jiaolai Basin during the Late Cretaceous: Implication for the coastal mountains along the East Asian margin [J]. Earth and Planetary Science Letters, 2016, 456: 112-123. doi: 10.1016/j.jpgl.2016.09.034
[56] Yang Y T. An unrecognized major collision of the Okhotomorsk Block with East Asia during the Late Cretaceous, constraints on the plate reorganization of the Northwest Pacific [J]. Earth-Science Reviews, 2013, 126: 96-115. doi: 10.1016/j.earscirev.2013.07.010
[57] Chen P J. Paleoenvironmental changes during the Cretaceous in eastern China[M]//Okada H, Mateer N J. Developments in Palaeontology and Stratigraphy. Amsterdam: Elsevier, 2000, 17: 81-90.
-
期刊类型引用(4)
1. 袁瑞，张昌民，赵芸，张莉，陈哲，张宝进，黄若鑫. 基于偏正态概率分布的粒度分布次总体分离及其沉积环境指示意义. 地质论评. 2022(03): 1033-1047 . 
百度学术
2. 王燕校，王永，姚培毅，田飞，袁路朋，叶梦旎. 白洋淀漕河全新世早期古洪水事件的沉积特征及气候背景. 岩石矿物学杂志. 2022(05): 916-928 . 百度学术
3. 曹胜伟，刘春雷，李亚松，李静，郝奇琛，高婕，董岩，陆晨明. 福建泉州湾近岸海域沉积物重金属来源分析与生态风险评价. 中国地质. 2022(05): 1481-1496 . 百度学术
4. 董续舒. 粒度分析在沉积岩中的应用. 西部资源. 2021(04): 4-5+8 . 百度学术
其他类型引用(6)
计量
- 文章访问数: 2029
- HTML全文浏览量: 514
- PDF下载量: 46
- 被引次数: 10
