Quaternary strata lithological characteristics and environmental geological significance of borehole WKZK01 in the eastern Songnen Plain
-
摘要: 松嫩平原东部第四纪地层特征研究对于了解区域气候变化和沉积环境演化具有重要意义。然而，目前关于该区的第四纪地层研究较少，从而限制了对该区气候环境变化机理的理解。文章选择松嫩平原东部望奎WKZK01钻孔岩芯为研究对象，在岩性特征分析的基础上，利用光释光、孢粉和粒度等测试方法，对比分析哈尔滨荒山剖面的相关成果，对钻孔第四纪地层划分及所反映的沉积环境进行了探讨。将望奎县第四纪地层由老至新划分为：中更新统荒山组、中更新统哈尔滨组、上更新统顾乡屯组，无全新统坦途组。钻孔揭露荒山组厚13.05 m，为湖泊相，孢粉带I−1记录了中更新世早期凉干气候；哈尔滨组厚24.05 m，为湖泊相，孢粉带I−2记录了中更新世晚期凉湿气候；顾乡屯组厚15.6 m，为河流相，孢粉带II记录了晚更新世冷湿气候。望奎WKZK01钻孔岩芯记录了松嫩古大湖周边面积较小的湖泊由湖心到湖滨最后消亡形成河流的环境演化过程。Abstract:
Objective Study of the Quaternary stratigraphical characteristics of the eastern Songnen Plain is of particular importance for gaining an understanding of regional climate change and the evolution of sedimentary environments. To data, however, there has been comparatively little research focusing on the Quaternary strata in this area, which has accordingly limited our comprehension of the mechanisms underlying changes in the regional climate and environment. Methods In this study we focused on the borehole WKZK01 in the eastern Songnen Plain. On the basis of lithological characteristic analysis, we obtained dating, sedimentary environment evolution, and climatic evolution information using optically stimulated luminescence, palynology, and laser grain size analytical methods, and based on comparisons with the results obtained from a typical Huangshan section, Harbin, we discuss the Quaternary stratigraphy of the borehole and the sedimentary environments it reflects. Results The Quaternary strata in Wangkui county are divided from old to new as follows: the middle Pleistocene Huangshan Formation, the middle Pleistocene Harbin Formation, and the upper Pleistocene Guxiangtun Formation, with an absence of a Holocene Tantu Formation. The Huangshan Formation is 13.05 m thick, with a median grain size (Md) ranging from 7.77 μm to 11.84 μm, and an average value of 9.62 μm. Pollen zone I-1 is dominated by xerophytic herbaceous pollen, with small amounts of mixed coniferous and broadleaf woody plant pollen, although lacks pollen derived from hydrophilic herbaceous plant. The Harbin Formation is 24.05 m thick, with a median grain size (Md) ranging from 9.67 μm to 14.98 μm, and an average value of 11.74 μm. Pollen zone I-2 is still dominated by xerophytic herbaceous pollen, with that derived from mixed coniferous and broadleaf woody plants as the secondary component, and the appearance of hydrophilic herbaceous plant pollen. The Guxiangtun Formation is 15.6 m thick, with a median grain size (Md) ranging from 11.61 μm to 19.91 μm, and an average value of 14.12 μm. Compared with the earlier stages, pollen zone II has a less abundant pollen content, although is characterized by an increase in hydrophilic herbaceous pollen, and a marked decline in the content of pollen derived from xerophytic woody and herbaceous plants. Conclusion The Huangshan Formation represents lacustrine deposits, with pollen zone I-1 providing evidence of a cool and dry climate during the early Middle Pleistocene. The Harbin Formation also comprises lacustrine deposits, with the pollen record in zone I-2 indicating a cool and humid climate during the late Middle Pleistocene. The Guxiangtun Formation comprises fluvial deposits, with pollen zone II recording a cold and wet climate during the Late Pleistocene. Significance The borehole WKZK01 provides a record of the environmental evolution of a small lake adjacent to the large Songnen paleolake, revealing a transition from the lake center to the margin and eventually disappearing with the formation of rivers. -
Key words:
- paleoenvironment /
- paleoclimate /
- optically stimulated luminescence /
- grain size /
- pollen /
- Songnen Plain
-
0. 引言
近年来随着西部地区的基础工程建设数量及规模不断增加，西部高原地区的季节性冻融循环效应的影响也随之增强，这对工程建设，尤其是隧道围岩稳定性提出了更高的要求。因此，开展冻融循环作用下花岗岩的细观特性及强度劣化性质研究对指导西部寒区基础工程建设至关重要。
岩石在冻融循环效应的作用下，细观结构会发生一系列变化，而在寒区工程建设中，细观结构对于岩体强度有着极大的影响。高精度CT扫描作为一种获取试样内部细观结构的基本技术手段，得到了广泛的应用，但在实际研究当中，仅从CT扫描结果无法直观地获得研究结论，往往需要与其他技术手段相结合进行多维度的分析论证。针对此问题，国内外学者对于岩石细观结构的演化规律以及相关的分析方法进行了诸多研究（Tian and Han，2017；延恺等，2018； Xia et al.， 2019；刘惠民等，2020；戚利荣等，2021；张艳博等，2021；朱楚雄等，2021；Shan et al.，2022； 戴俊等，2022；张慧梅等，2022；张慧梅和王云飞，2022；张光明，2023），研究构建出以高精度CT扫描为主要技术手段，以CT扫描为基础，与神经网络、数值分析、3D打印、分形理论等技术手段结合，探讨岩石在各物理场、化学场拟合条件下的细观结构演化规律的方法体系。
岩石的强度降低与结构变化存在着明确的联系。张超等（2023）、郝志远（2022）在室内试验基础上，结合统计损伤理论，建立了一套岩石脆延转化统计损伤本构模型，以优化损伤变量参数的形式将微观结构的变化与岩石力学强度结合。在物理场耦合条件下岩石的微观结构的损伤研究方向上，于洪雯等（2021）、缪澄宇等（2021）分别从微波辐射、含水率等耦合方向对岩石结构与强度之间的关系做出了研究（蔡沛辰等，2021）。杨鸿锐（2022）指出岩石冻融劣化的主要原因是未冻水迁移、积聚，水分结晶膨胀导致岩石孔隙扩展、联通，提出温度区间是影响冻融损伤的关键性因素。
以上研究大多是在常规条件下针对岩石结构或者力学性质展开探讨，分布在高寒地带的花岗岩存在普遍的冻融循环劣化效应，近年来花岗岩等岩石对于冻融循环效应敏感性的研究需求逐步体现，花岗岩具有原始结构致密，岩石强度较高等特点，而现有对于冻融循环条件下岩石力学特性的研究多以砂岩等冻融循环敏感性较高的岩石种类为主，针对高寒地区的花岗岩缺乏针对性较强的冻融循环力学特性研究及微观特征分析，其冻融条件下结构演化与强度之间的关系还不够明确。
因此，文章以川西地区某隧道花岗岩为研究对象，利用CT扫描技术，对冻融后的花岗岩进行扫描，对扫描图层利用阈值分割进行二值化处理，堆叠得到样品内外结构的高分辨3D数据及影像。结合分形理论计算图像计盒维数并由此对图像复杂度做出量化判断，由此对冻融循环对花岗岩内部结构演化分布特点进行分析。在此基础上，结合传统力学试验，揭示其强度演化规律，探究结构演化与强度之间的关系。
1. 测试仪器及试验方案
在青藏高原地区，随着海拔的攀升，季节交替导致的冻融循坏效应十分显著。考虑到研究对象在冻融环境下的代表性，在康定地区某隧道采集了二长花岗岩样品，研究冻融循环对细观结构和力学特征的影响，样品取自隧道施工开挖出的新鲜岩块，整体结构致密，完整性较好，风化程度低，原始结构面不发育。
依照《水利水电工程岩石试验规程（SL/T264-2020）》 ，岩样加工为50 mm×100 mm标准圆柱试样。
1.1 岩石薄片鉴定
参照《地质矿产实验室测试质量管理规范》（DZ /T 0130.9—2006)，选用现场新鲜岩样切割制成薄片，通过偏光显微镜，鉴定岩石样本矿物的组成及其特征，显微照片和鉴定结果如图1所示。
图 1 花岗岩薄片及矿物组成Q—石英； Pl—斜长石； Bi—黑云母； Mic—微斜长石Figure 1. Thin section and mineral composition of the granite sample (a) Microscopic image of thin section; (b) Pie chart of mineral compositionQ–quartz; Pl–plagioclase; Bi–biotite; Mic–microcline由图1可知，岩石呈块状构造，具有似斑状粗粒不等粒花岗结构，局部见交代蠕虫结构。似斑晶矿物主要为碱性长石；其他矿物粒径0.25~4.0 mm为主，矿物成分主要为石英、斜长石、碱性长石（微斜长石、条纹长石），次要矿物为黑云母、绿帘石，副矿物有磷灰石、锆石、黄铁矿。岩石总体较为完整，但可见有一组平行发育的微裂隙，为岩石中的劈理面，基本未充填；定性为：似斑状粗粒不等粒黑云二长花岗岩。
1.2 冻融循环试验方案
依据《工程岩体试验方法标准》（GBT 50266—2013），试样冻融前，先在真空条件下对试样进行48小时的强制饱水。饱水后，利用DBD-5型全自动冻融循环机进行冻融循环试验（图2），冻融循环次数设定为20次。冻融循环过程中，试样处于无约束状态。样品采集地的年最低气温能达到−18 ℃，最高温度可达到23 ℃，因此，在确定温度区间时，为保证试验效果且能达到研究区温度要求，初步设定冻融循环温度区间为−25 ℃~25 ℃，每次循环9小时，降温0.5小时，保持4小时，接着升温0.5小时，保持4小时，达到循环次数后，机器回到常温，试验结束，循环路径如图3所示。
从宏观上看，对比未进行冻融的花岗岩试样，20次冻融的花岗岩试样整体形态上无太大变化，试样整体完好，表面无裂痕产生，仅存在极少许片状剥落。
1.3 计算机层析扫描（CT）试验
CT扫描实验采用nanoVoxel-4000系列微纳米高分辨CT扫描仪，结构如图4所示，其最高分辨率可达到0.5 μm。为了对试样细观结构进行详细的观察，此次扫描试验预设层数为2500层。试样分层如图5，底层编号为0，顶部编号为2500。
CT扫描试验中，样品内部结构的相对密度与CT图像的灰度成正相关，扫描后得到试样的灰度四视图，通过灰度能直观地观察样品内部孔隙及颗粒分布特征，密度大的岩石颗粒以高亮点显示，密度较小的孔隙结构则由暗点显示，由于肉眼无法准确察觉灰度图像中的微小变化，因此对得到的灰度图像进行RGB通道渲染处理，使得深色部分更为突出，得到处理后的图像如图6所示。
图 6 试样灰度四视图a—俯视图；b—正视图；c—左侧视图；d—立体图Figure 6. Four-view grayscale of the sample(a) Top view; (b) Front view; (c) Left view; (d) 3D stereogram通过观察图6，在试样经历20次冻融循环试验后，试样有一定程度的孔隙分布，在XY平面，孔隙分布较为均匀，整个圆形面上皆有分布，孔隙直径0.1~0.6 mm。相对较大的孔隙主要分布于图像XY面中上部。从YZ以及XZ平面中，也存在少许孔隙分布，花岗岩在冻融循环过程中出现小规模表面剥落以及颗粒脱落现象。
1.4 岩石力学试验
张尧等（2008）认为岩石蠕变是围岩变形失稳的重要原因之一，室内蠕变试验具有能够长期观测、较严格控制试验条件、重复次数多等优点。
对冻融前后的花岗岩试样分别进行了单轴压缩试验，三轴压缩试验以及基于三轴抗压强度进行的卸荷蠕变试验以探究冻融循环对花岗岩力学性质的影响。试验平台为成都理工大学地质灾害防治与地质环境保护国家重点试验室自主研发的YSJ-01-00型岩石三轴蠕变试验机（图7），其轴向荷载范围0~600 kN，围压范围0~30 MPa，满足此次试验的需求。
结合单轴压缩试验结果得到花岗岩冻融前后基本力学参数（表1）。可以看出，在经历20次冻融循环后，花岗岩的各物理参数均存在一定程度的降低，单轴抗压强度降幅为16.9%。
表 1 花岗岩物理力学参数Table 1. Physical and mechanical parameters of granite冻融循环次数 泊松比 弹性模量/GPa 单轴抗压强度/MPa 0 0.19 11.2 106 20 0.22 10.4 88 2. 基于CT图像的岩石细观及孔隙特征分析
2.1 CT图像二值化图像处理
为了能够对试样内部孔隙与岩石颗粒的分布特点进行定量分析，参考现有研究成果（王慧琴，2021；刘天舒等，2022；李茂民和邹臣嵩，2022），采用阈值分割法对CT 扫描模型进行二值化处理。其原理为：设定灰度值阈值$\alpha $，若图像中某像素点的灰度值小于该阈值，则将该像素灰度值设置为0或225，否则设置为225或0。二值化处理函数表达式如式1。
f(X,Y)={0X<α225X<α (1) 对每一层CT层面利用阈值分割方法进行二值化处理，再按照层面编号由下至上堆叠，便得到二值化处理后的CT扫描模型，对模型中孔隙部分进行染色，得到图8冻融循环20次后试样二值化处理后CT立体模型，蓝色部分对应模型中的孔隙。
通过二值化CT模型，可以看出试样内部整体孔隙分布较为广泛，孔隙直径处于较低水平，孔隙在试样内部集中度高于试样外围，试样内部分布较多微孔隙。考虑到试样由整体均一性较好的岩块切割而成，该现象的出现只与后续试验相关，因此，该现象与冻融循环效应有着紧密联系。
2.2 基于CT图像的岩石孔隙率分析
基于图8的立体模型，利用计算机进行计算得到试样逐层孔隙率，冻融循环20次后的花岗岩试样中，孔隙占试样的总体积为0.134%。整体孔隙率处于较低水平，且从CT立体模型来看，几乎不存在孔隙的贯通。为了进一步分析Z轴方向孔隙的分布规律，对每一CT扫描层进行孔隙率统计，得到Z方向逐层面孔隙率，示于图9。
冻融循环后试样渗透率为0.023×10−3 μm2，岩石渗透水平小幅增加但仍处于较低水平，因此在0至20次冻融循环期间岩石内部孔隙贯程度较低，同时冻融循环后试样内部孔隙平均孔喉半径为0.015 μm，最大孔喉半径为0.053 μm，与CT模型显示的孔隙以微小孔隙的结论一致。为了进一步分析Z轴方向孔隙的分布规律，对每一CT扫描层进行孔隙率统计，得到Z方向逐层面孔隙率，示于图9。
分析图9数据，逐层孔隙率随着层面编号的增加，整体呈现出逐渐上升的趋势，这与CT立体模型中蓝色标记密集度分布规律保持较好的一致性。试样孔隙率最高层出现在2436层，达到 0.254%。此外，曲线走势在24层、1223层、2436层出现三次峰值。据此，结合图7、图8、表1可以推断，冻融循环效应在影响花岗岩细观结构时，会导致花岗岩内部孔隙率的整体上升，与此同时，在花岗岩原始孔隙或者结构面相对密集的部位，冻融循环效应的影响更为显著。随着冻融循环次数的增加，最有可能于孔隙率较高的部位率先形成贯通的结构面（田彦德，2019），在此次试验中体现于试样底部，以及中上部最可能率先发生结构面的贯通。
2.3 基于分形理论的岩石细观结构分析
2.3.1 岩石计盒维数的计算
分形理论运用连续维度的思路对物体属性及状态进行描述，分形维数已发展成为描述物体复杂度的重要参数，而分形维数又可分Hausdorff维数DH、相似维数DS、信息维数Di等（李后强和汪富泉，1992；乔力江汉和何克焓，2021；张子涵等，2022）。此次研究使用的计盒维数为Hausdorff维数的等价位数。
计盒维数的计算思路为选取半径为D的盒子将图像中的非孔隙部分全部覆盖，所需要的最小盒子数量为Nd，对二者进行公式拟合，使其满足如下关系（张慧梅等，2022）。
Nd=A×DB (2) 其中，A、B为参数。
对式（2）进行变形，得到
B=−lnNdlnD+lnAlnD (3) 其中，$ \dfrac{\ln{N}_{d}}{\ln D} $ 以及 $ \dfrac{\ln A}{\ln D} $ 即为计盒维数。
2.3.2 分形维数分析
参考Matlab图形处理领域相关研究成果（郝凯越等，2020），进行Matlab程序编程，利用程序提取图像二值化后的二维矩阵数据，利用fraclab插件进行不同半径D下的盒数Nd的数目统计，根据D值的不同，可以得到多个不同的Nd数目，然后在双对数坐标系下画出数据对，利用最小二乘法对数据进行拟合，得到一元一次函数的斜率 k，其相反数即为图像的计盒维数。利用fraclab拟合得到的20次冻融循环花岗岩试样的Nd -D曲线如图10所示。
图 10 20次冻融循环花岗岩试样的Nd -D曲线Figure 10. Nd -D curve of granite samples after 20 freeze-thaw cycles在图10中的红色曲线为二分法拟合的一次函数，其斜率的相反数为1.88，两条曲线拟合度较高，说明花岗岩表面均具有典型的分形特征，该冻融试样图像的计盒维数为1.88。横向对比现有对花岗岩的分形研究（谢韬，2013；李德建，2010；樊水龙，2020），完整度较好的花岗岩计盒维数一般处于1.85~1.95之间。此次试验冻融循环20次的花岗岩试样，其内部结构有轻微劣化，复杂度出现小幅降低，但试样整体完整度仍然较好。CT图像表明冻融过后，试样内部新萌生较多微孔隙，但岩石分形维数虽出现小幅降低但仍保持在较高水平，说明微孔隙的萌生，对岩石结构复杂度的影响较小。由此推断，低次数冻融循环会使得花岗岩类致密岩石内部萌生微小孔隙，但对整体结构复杂度影响较小，随着冻融循环的继续进行，出现微孔隙不断聚集连通现象时，试样内部复杂度可能会出现较大幅度的下降。
3. 冻融循环作用下强度特征演化
3.1 花岗岩卸荷蠕变特征
赵娜等（2023）认为岩石蠕变破裂是由岩石内部围观破裂逐渐向宏观发展的过程。因此，为了探究研究区花岗岩结构演化与强度变化的特点，开展花岗岩卸荷蠕变试验。
在常规三轴试验的基础上，取三轴抗压强度的75%作为固定轴向荷载开展卸荷蠕变试验，以契合边坡及隧道的卸荷工况，为保证试样能在3~5级荷载时破坏，且围压卸荷梯度有一定幅度，确定试验初始围压为7.5 MPa，围压卸荷梯度为Δσ =1.5 MPa，每级荷载保持24 h，直到试样发生破坏。得到20次冻融前后花岗岩试样卸荷蠕变曲线如图11所示。
图 11 冻融前后花岗岩卸荷蠕变时间−应变曲线Figure 11. Time-strain curves of unloading creep of granite before and after freeze-thawing(a) Before freeze-thawing; (b) After 20 freeze-thaw cycles对比冻融前后的花岗岩卸荷蠕变曲线，冻融前试样在第四级荷载发生破坏，冻融后则是在第五级荷载下发生破坏。冻融前试样在轴向应变为2.12%时进入加速蠕变，而冻融后试样则是在轴向应变达到2.52%时才进入加速蠕变。冻融后的岩石CT模型显示，试样在经历冻融循环作用后，内部微孔隙发育较多，使试样在压缩过程具有更大的位移空间，表现为黏性增加，因此进入加速蠕变阶段时相对于冻融前的试样具有较大的应变阈值。
3.2 长期强度分析
长期强度常作为衡量工程是否长期稳定的重要指标，对于边坡及隧洞施工有着重大的意义。在蠕变试验中，当附加应力水平大于或等于某一定值时，试样的蠕变变形会随时间逐渐增大，且岩石的蠕变速率会在某一阶段迅速增加，试样最终发生破坏。这个应力水平的临界值就认为是岩石的长期强度（Innocente et al.， 2021；Liu et al.， 2022； 唐志强等，2022）。
结合岩石蠕变试验数据，利用Bolziman线性叠加的方式确定冻融前后花岗岩的长期强度，得到冻融前后花岗岩偏应力−应变曲线如图12所示。
图 12 冻融前后花岗岩偏应力−应变曲线a—冻融前；b—20次冻融循环后Figure 12. Deviatoric stress–strain curve of granite before and after freeze-thawing(a) Before freeze-thawing; (b) After 20 freeze-thaw cycles通过观察图12，可以看出图像走势可分为两部分，一部分是图像左侧斜率较大的部分，另一部分是右侧斜率较小的部分。取曲线走势转折点点对应的偏应力水平即为该条件下花岗岩的长期强度。因此，得到冻融前花岗岩的长期强度分别为128.21 MPa，冻融后为102.50 MPa。由此经过20次冻融循环后的花岗岩的长期强度会出现一定幅度的衰减，衰减幅度为20%。
从孔隙率上看，20次冻融循环后，花岗岩孔隙率虽有小幅上升，但分形维数显示岩石整体完整性仍较好，结构复杂度变化不大。结合岩石CT模型中的微孔隙发育特征可以推断，在经历过20次冻融循环后，花岗岩内部结构的变化以微小孔隙的发育为主，不伴随大规模的贯通，因此岩石整体孔隙率小幅上升但复杂度不发生较大变化。与此同时，岩石长期强度有着较大幅度的变化，说明微孔隙的产生对花岗岩的力学性质有着较大的影响。
综合上述分析，经历20次冻融循环后，岩石试样在更低的偏应力水平下即进入加速蠕变阶段，与此同时，在进入加速蠕变时具有更大的变形量，应变阈值提升19.3%，即在强度更低的同时会发生更大的变形。而在自然条件下，冻融循环往往是长时间、多次数的，结合 CT试验，冻融循环后，花岗岩内部会萌生较多对强度影响较大的微小孔隙，因此在几十次甚至上百次的冻融循环条件下，花岗岩的力学性质必然会出现较大幅度的折减，其加速蠕变偏应力阈值会出现较大幅度下降且变形量会明显增大，体现在边坡上即应该警惕冻融循环效应引起的滑坡及危岩体灾害的发生；体现在围岩上，应警惕围岩大变形甚至结构失稳的发生。在冻融循环前期，由于微孔隙的萌生在结构复杂度上未产生宏观影响，但强度会发生劣化，此类边坡及围岩应从岩体结构及强度两方面做好监测工作，预防加速蠕变阶段产生的安全隐患。
4. 结论
对采集自川西地区典型地段的花岗岩进行了20次冻融循环，利用岩石CT层析扫描技术对每次试样分2500层进行了层析扫描，利用阈值分割方法对扫描的图像进行二值化处理后，由下到上进行堆叠，建立了精度0.3 μm的岩石CT扫描三维模型，为研究冻融循环对岩石细观结构的影响提供了基础，主要取得了以下认识。
（1）通过分析CT模型内部孔隙分布特征及压汞试验结果，发现冻融循环效应在影响花岗岩细观结构时，会导致花岗岩内部孔隙率的整体上升，但岩石渗透性变化不大，岩石渗透率仅上升0.003×10−3 μm2。内部孔隙发育不均匀，试样整体结构改变以萌生较多新的微孔隙为主。
（2）冻融循环后岩石图像分形计算所得的计盒维数为1.88，内部结构复杂度有所下降，但岩石整体完整性仍然较好，分形维数仍保持在较高水平。微孔隙的萌生，对岩石结构复杂度的影响较小。对于低次数短时冻融岩石内部结构的分析有一定指导意义。
（3）分形研究显示，20次冻融循环并未导致花岗岩的结构复杂度发生较大变化，同时试样整体力学特性出现下降，黏性增加以及长期强度出现较大幅度的衰减，进入蠕变试验阶段的应变阈值提高。因此，在评价此类原生结构较致密的岩石的安全性时，仅从结构上进行考量与实际情况往往会出现偏差，应结合必要的强度指标综合评估。
（4）此次研究对象二长花岗岩在经历冻融循环后，在强度更低的同时会发生更大的变形，施工方应从岩体结构及强度两方面做好监测工作，预防加速蠕变阶段可能出现的安全隐患。
-
图 1 松嫩平原东部钻孔WKZK01区域概况图与岩性柱状图
a—区域概况图；b—WKZK01钻孔岩性柱状图
Figure 1. General map and lithological histogram of borehole WKZK01 in the eastern Songnen Plain
(a) General map; (b) Lithological histogram of borehole WKZK01
图 2 松嫩平原东部WKZK01钻孔岩性三角图
Figure 2. Lithological triangle map of borehole WKZK01 in the eastern Songnen Plain
图 4 WKZK01钻孔粒度参数及不同粒级含量变化曲线
Figure 4. Variation curves of grain size parameters and contents of different grades in borehole WKZK01
图 5 WKZK01钻孔顾乡屯组、哈尔滨组和荒山组地层样品的粒度参数散点图、C−M图
a—偏度−平均粒径散点图；b—标准偏差−平均粒径散点图；c—峰度−平均粒径散点图；d—C−M图
Figure 5. Grain size parameters of formation samples from the Guxiangtun, Harbin and Huangshan Formations in borehole WKZK01
(a) Skewness-mean scatter plot; (b) Standard deviation-mean scatter plot; (c) Kurtosis-mean scatter plot; (d) C(First percentile)-M(Median) plot
图 6 WKZK01钻孔岩芯粒度参数特征及其与哈尔滨荒山剖面粒度特征对比
Figure 6. Grain size parameter characteristics of borehole WKZK01 and comparison with the grain size characteristics of the Huangshan section in Harbin
表 1 钻孔WKZK01钻孔分层及岩性特征
Table 1. Lithological and stratigraphical characteristics of borehole WKZK01
层号 层顶深度/m 层底深度/m 层厚/m 岩性特征 1 0 0.6 0.6 粉砂，棕黑色，团粒状结构，略为紧实，略为湿润，少量孔隙，少量根系 2 0.6 1.0 0.4 砂质粉砂，灰黑色，团粒状结构，略为紧实，略为湿润，少量孔隙，少量根系 3 1.0 7.8 6.8 粉砂，暗棕色，少量孔隙，少量根系，含铁染锈斑；2.7~6.5 m处为灰黄色，未见铁锰结核，颜色较为单一；6.5~7.4 m处为黄褐色；其中，6.5~6.8 m处为暗棕色，6.8~7.0处为黄褐色，含铁锰结核与锈斑，往下铁锰结核减少，颜色略带青灰色，反映氧化−还原环境 4 7.8 8.2 0.4 黏土质粉砂，棕黑色，块状结构，略为湿润，略为紧实 5 8.2 15.6 9.4 粉砂，黄褐色，少见铁锰结核与锈斑；该层从上往下颜色由深变浅。11.7~13.1 m处为棕黑色，见较多铁锰结核；该地层与上覆地层颜色界限较为清楚，可能为古土壤层；13.1~15.6 m处为暗棕色，块状结构，略为紧实，略为湿润；该地层较上覆地层颜色偏浅，但变化不明显 6 15.6 17.6 2.0 粉砂，暗棕色与棕黑色条带相间分布，块状结构，略为紧实，略为湿润 7 17.6 18.4 0.8 黏土质粉砂，黄褐色，块状结构，略为紧实，略为湿润，见少量红褐色铁结核 8 18.4 19.0 0.6 粉砂，黄褐色，见少量红褐色铁结核 9 19.0 19.6 0.6 黏土质粉砂，黄褐色，块状结构，略为紧实，略为湿润，见少量红褐色铁结核 10 19.6 20.2 0.6 粉砂，黄褐色，含黑色锰结核；该层与上覆地层区别为含较多黑色锰结核，见少量红褐色铁结核 11 20.2 20.8 0.6 黏土质粉砂，该层与上覆地层区别为含较多黑色锰结核，见少量红褐色铁结核 12 20.8 27.0 6.2 粉砂，黄褐色，含黑色锰结核；23.1~23.4 m处为棕黑色，可能为古土壤，往下颜色逐渐过渡变化 13 27.0 28.2 1.2 黏土质粉砂，黄褐色，含黑色锰结核；该层含较多黑色斑点状锰结核，见少量红褐色铁结核 14 28.2 39.6 11.4 粉砂，黄褐色，块状结构，紧实，略为湿润，含较多黑色斑点状锰结核；38.1~38.7 m处颜色偏暗，为棕黑色，发育灰白色斑点状、角砾状钙盐粉末，基质为灰黑色网纹状填充 15 39.6 40.2 0.6 黏土质粉砂，暗棕色，含黑色斑点状锰结核和红褐色铁结核，微细层理较为发育 16 40.2 40.6 0.4 粉砂，暗棕色，含黑色斑点状锰结核和红褐色铁结核，40.2~40.4 m颜色偏暗，微细层理较为发育 17 40.6 41.0 0.4 黏土质粉砂，暗棕色，含黑色斑点状锰结核和红褐色铁结核，微细层理较为发育 18 41.0 42.8 1.8 粉砂，暗棕色，含黑色斑点状锰结核和红褐色铁结核，微细层理较为发育 19 42.8 43.2 0.4 黏土质粉砂，暗棕色，含黑色斑点状锰结核和红褐色铁结核，微细层理较为发育 20 43.2 43.6 0.4 粉砂，暗棕色，含黑色斑点状锰结核和红褐色铁结核，微细层理较为发育 21 43.6 50.0 6.4 黏土质粉砂，暗棕色，含黑色斑点状锰结核和红褐色铁结核，微细层理较为发育；47.2~47.4 m处颜色偏暗，为棕黑色，发育灰白色斑点状、角砾状钙盐粉末，基质为灰黑色网纹状填充 22 50.0 50.6 0.6 粉砂，暗棕色，块状结构。含黑色斑点状锰结核和红褐色铁结核，微细层理较为发育；底部50.3~50.6 m处颜色偏暗，为棕黑色 23 50.6 51.4 0.8 黏土质粉砂，暗棕色，块状结构，含黑色斑点状锰结核和红褐色铁结核，微细层理较为发育；底部颜色偏暗，为棕黑色，发育灰白色斑点状、角砾状钙盐粉末，基质为灰黑色网纹状填充 24 51.4 52.2 0.8 粉砂，暗棕色，块状结构，含黑色斑点状锰结核和红褐色铁结核，微细层理较为发育 25 52.2 52.7 0.5 黏土质粉砂，暗棕色，含黑色斑点状锰结核和红褐色铁结核，微细层理较为发育；底部颜色偏暗，为棕黑色，发育灰白色斑点状、角砾状钙盐粉末，基质为灰黑色网纹状填充 
下载: 导出CSV
表 2 WKZK01钻孔光释光测年样品信息及测年数据
Table 2. Sample information and dating data of optical stimulated luminescence (OSL) of borehole WKZK01
样品编号 埋深/m U/（μg/g） Th/（μg/g） K/% 等效剂量/Gy 误差/Gy 年剂量/ （Gy/ka） 含水量/% 年龄/ka WKZK01−1 2.33 1.68 8.32 1.85 183.81 7.90 2.82 24.91 65.2±2.8 WKZK01−2 5.08 1.22 7.25 1.90 199.13 6.13 2.51 28.66 79.5±2.4 WKZK01−3 8.84 1.31 7.17 1.61 215.21 9.27 2.50 17.80 86.0±3.7 WKZK01−4 12.19 1.43 6.01 1.96 236.06 12.77 2.60 22.52 90.7±4.9
下载: 导出CSV
-
[1] Bureau of Geology and Mineral Resoures of Heilongjiang Province, 1993. Regional geology of Heilongjiang province[M]. Beijing: Geology Press. (in Chinese) [2] Editorial Board of Chinese Vegetation, 1980. Vegetation of China[M]. Beijing: Science Press. (in Chinese) [3] FENG Z Z, 1993. Sedimentary petrology[M]. 2nd ed. Beijing: Petroleum Industry Press. (in Chinese) [4] GAO H Y, XIE Y Y, CHI Y P, et al., 2023. Sedimentary composition and detrital zircon U-Pb dating of the sand-gravel profile in Binxian County, Harbin: Indications for sedimentary environment and regional tectonic evolution[J]. Chinese Journal of Geology, 58(4): 1354-1372. (in Chinese with English abstract [5] GAO W, JIA D C, LI T L, et al., 2013. Stratigraphic boundary between Middle and Upper Pleistocene Series in Eastern Songliao Plain: age of optical luminescennce and element changes[J]. Journal of Jilin University (Earth Science Edition), 43(6): 1889-1896. (in Chinese with English abstract [6] HE H C, DING H Y, ZHANG Z K, et al., 2005. Grain-size characteristics and their environmental significance of Hongze Lake Sediments[J]. Scientia Geographica Sinica, 25(5): 590-596. (in Chinese with English abstract [7] HOU X R, XIE Y Y, KANG C G, et al., 2023. Sedimentological characteristics of the Baitushan Formation in the Eastern Foothills of the Great Xing’an Range: implications for stratigraphic divisions[J]. Acta Sedimentologica Sinica, 41(3): 720-734. (in Chinese with English abstract [8] LIN N F, TANG J, 2005. Study on the environment evolution and the analysis of causes to land salinization and desertification in Songnen Plain[J]. Quaternary Sciences, 25(4): 474-483. (in Chinese with English abstract [9] LENG Y K, XIE Y Y, KANG C G, et al., 2023. Sedimentary characteristics and environmental significance of the Juren Sandy Gravel Profile in Harbin[J]. Acta Sedimentologica Sinica, 41(2): 472-484. (in Chinese with English abstract [10] LI W Y, YAO Z J, 1990. A study on the quantitative relationship between pinus pollen in surface sample and Pinus vegetation[J]. Acta Botanica Sinica, 32(12): 943-950. (in Chinese with English abstract [11] LI X X, MAO D L, LAI F B, et al., 2024. Grain size characteristics and sand source analysis of three aeolian landforms in the lower reaches of the Qira River floodplain[J]. Arid Zone Research, 41(8): 1413-1422. (in Chinese with English abstract [12] LI Z Z, LING Z Y, CHEN X L, et al., 2010. Late holocene climate changes revealed by Grain-size cycles in takemukul desert in Yili of Xinjiang[J]. Scientia Geographica Sinica, 30(4): 613-619. (in Chinese with English abstract [13] LINK A G, 1966. Textural classification of sediments[J]. Sedimentology, 7(3): 249-254. doi: 10.1111/j.1365-3091.1966.tb01598.x [14] LIU D Y, LI W R, PENG S S, et al., 2010. Current application of grain size analysis in Chinese loess paleoclimatic study[J]. Periodical of Ocean University of China, 40(2): 79-84. (in Chinese with English abstract). [15] LIU Z R, XUE H Y, WANG C S, 2021. Late Quaternary depositional characteristics and environment significance of the Xibozhang section in Baoding, central Hebei Plain, China[J]. Journal of Geomechanics, 27(6): 1011-1023. (in Chinese with English abstract [16] MA Y F, ZHAN T, YANG Y, et al., 2021. The indication of chroma characteristics and its palaeoclimatic significance in the Tianhengshan (THS) core from the eastern part of the Northeast China Plain to the evolution of Songnen paleo-lake[J]. Acta Geologica Sinica, 95(11): 3519-3531. (in Chinese with English abstract [17] QIU S W, WANG X K, ZHANG S Q, et al., 2012. The evolution of the Large Paleolake in Songliao plain and its formation[J]. Quaternary Sciences, 32(5): 1011-1021. (in Chinese with English abstract [18] QIU S W, WANG X K, MAKHINOV A N, et al., 2014. Summary of the paleodrainage pattern changes in the Northeast China Plain and its neighboring areas[J]. Acta Geographica Sinica, 69(11): 1604-1614. (in Chinese with English abstract [19] SHEPARD F P, 1954. Nomenclature based on sand-silt-clay ratios[J]. Journal of Sedimentary Research, 24(3): 151-158. [20] SUN J H, XIE Y Y, KANG C G, et al., 2022. Stratigraphic properties of the baitushan formation in Ping’an town, the eastern foot of the Great Hinggan Mountains: an indication of provenance and sedimentary environment[J]. Journal of Stratigraphy, 46(2): 196-208. (in Chinese with English abstract [21] SUN J Z, WANG Y Z, ZHANG Q Y, 1982. Classification of Quaternary stratigraphy in Songliao Plain-Application of several chronological methods[J]. Journal of Chang'an University (Earth Science Edition)(2): 79-91, 10. (in Chinese [22] SUN Y, XIE Y Y, CHI Y P, et al., 2022. Stratigraphic characteristics of the Baitushan Formation in Longjiang county, eastern foothills of the great Xing'an Range, China: chemical weathering, sediment cycling, source-sink system and sedimentary environment[J]. Mountain Research, 40(1): 14-28. (in Chinese with English abstract [23] TIAN M Z, CHENG J, 2020. Quaternary geology and geomorphology[M]. 2nd ed. Beijing: Geology Press. (in Chinese) [24] WANG E B, 2012. Quaternary sedimentary sequences and their environmental significance in Harbin Area[D]. Changchun: Jilin University. (in Chinese with English abstract [25] WANG J X, XIE Y Y, KANG C G, et al., 2020. The indication of the heavy mineral characteristics of the core in Harbin Huangshan to the Quaternary drainage evolution of Songhua River[J]. Quaternary Sciences, 40(1): 79-94. (in Chinese with English abstract [26] WANG X Y, WU L, ZHANG G S, et al., 2008. Characteristics and environmental significance of magnetic susceptibility and grain size of lake sediments since Holocene in Chaohu Lake, Anhui Province[J]. Scientia Geographica Sinica, 28(4): 548-553. (in Chinese with English abstract [27] WANG Y, DONG J, YANG J S, 2020. Quaternary Stratigraphy of the Huangshan Section in Harbin[J]. Earth Science, 45(7): 2662-2672. (in Chinese with English abstract [28] WU P, XIE Y Y, KANG C G, et al., 2020a. The genesis of Huangshan loess in Harbin: integrated evidence from grain size, geochemistry, magnetization, sedimentation and landform[J]. Acta Geoscientica Sinica, 41(3): 420-430. (in Chinese with English abstract [29] WU P, XIE Y Y, KANG C G, et al., 2020b. The capture of the Songhua River system in the late Early Pleistocene: geochemical and sedimentological records[J]. Acta Geologica Sinica, 94(10): 3144-3160. (in Chinese with English abstract [30] WU X H, PU Q Y, QIAN F, et al., 1984. Preliminary study on the quaternary magnetostratigraphy of the Songliao plain in North-East China[J]. Marine Geology & Quaternary Geology, 4(2): 1-13. (in Chinese with English abstract [31] XIE Y Y, YUAN F, ZHAN T, et al., 2018. Geochemistry of loess deposits in northeastern China: constraint on provenance and implication for disappearance of the large Songliao palaeolake[J]. Journal of the Geological Society, 175(1): 146-162. doi: 10.1144/jgs2017-032 [32] XU Y Y, XIE Y Y, KANG C G, et al., 2022. Drainage evolution of the Songhua River system in the Early Pleistocene: evidence from TIMA minerals and geochemistry[J]. Chinese Journal of Geology, 57(1): 190-206. (in Chinese with English abstract [33] ZENG J L, CHI Y P, XIE Y Y, et al. , 2024. Stratigraphic attribute characteristics of sand and gravel accumulation of Sanjiedi Profile in Harbin, China[J/OL]. Acta Sedimentologica Sinica, 1-21[2024-11-27]. https://doi.org/10.14027/j.issn.1000-0550.2024.026. (in Chinese with English abstract [34] ZHAN T, ZENG F M, XIE Y Y, et al., 2018. Grain size characteristics of Tianhengshan core and their indications for stratigraphic division in the eastern part of the northeast plain of China[J]. Journal of Geomechanics, 24(4): 515-521. (in Chinese with English abstract [35] ZHAN T, ZENG F M, XIE Y Y, et al., 2019. Magnetostratigraphic dating of a drill core from the Northeast Plain of China: implications for the evolution of Songnen paleo-lake[J]. Chinese Science Bulletin, 64(11): 1179-1190. (in Chinese with English abstract doi: 10.1360/N972018-01212 [36] ZHANG M, XIE Y Y, KANG C G, et al., 2021. Chemical weathering characteristics of the Loess-Paleosol Sequences in Harbin Huangshan Rock Core: implication for formation environment of the Paleosol[J]. Acta Pedologica Sinica, 58(3): 673-684. (in Chinese with English abstract [37] ZHANG X, ZHENG Z, HUANG K Y, et al., 2020. Sensitivity of altitudinal vegetation in Southwest China to changes in the Indian summer monsoon during the past 68000 years[J]. Quaternary Science Reviews, 239: 106359. doi: 10.1016/j.quascirev.2020.106359 [38] ZHANG Y X, CHI Y P, XIE Y Y, et al., 2020. Organic carbon isotope composition of Harbin loess since the Mid-Pleistocene and its paleoclimatic significance[J]. Acta Geoscientica Sinica, 41(4): 525-534. (in Chinese with English abstract [39] 冯增昭,1993. 沉积岩石学[M]. 2版. 北京:石油工业出版社. [40] 高宏宇,谢远云,迟云平,等,2023. 哈尔滨宾县砂砾石剖面沉积物组成及碎屑锆石U-Pb测年:对沉积环境和区域构造演化的指示[J]. 地质科学,58(4):1354-1372. doi: 10.12017/dzkx.2023.075 [41] 高文,贾大成,李桐林,等,2013. 松嫩平原东部中、晚更新世地层界限:光释光年龄及元素变化[J]. 吉林大学学报(地球科学版),43(6):1889-1896. [42] 何华春,丁海燕,张振克,等,2005. 淮河中下游洪泽湖湖泊沉积物粒度特征及其沉积环境意义[J]. 地理科学,25(5):590-596. doi: 10.3969/j.issn.1000-0690.2005.05.013 [43] 黑龙江省地质矿产局,1993. 黑龙江省区域地质志[M]. 北京:地质出版社. [44] 侯心茹,谢远云,康春国,等,2023. 大兴安岭东麓白土山组地层的沉积学特征:对地层划分的指示[J]. 沉积学报,41(3):720-734. [45] 林年丰,汤洁,2005. 松嫩平原环境演变与土地盐碱化、荒漠化的成因分析[J]. 第四纪研究,25(4):474-483. doi: 10.3321/j.issn:1001-7410.2005.04.011 [46] 冷宇坤,谢远云,康春国,等,2023. 哈尔滨居仁砂砾石剖面沉积特征及其环境意义[J]. 沉积学报,41(2):472-484. [47] 李文漪,姚祖驹,1990. 表土中松属花粉与植物间数量关系的研究[J]. 植物学报,32(12):943-950. [48] 李鑫鑫,毛东雷,来风兵,等,2024. 策勒河下游风积地貌沉积物粒度特征及沙源分析[J]. 干旱区研究,41(8):1413-1422. [49] 李志忠,凌智永,陈秀玲,等,2010. 新疆伊犁河谷晚全新世风沙沉积粒度旋回与气候变化[J]. 地理科学,30(4):613-619. [50] 刘冬雁,李巍然,彭莎莎,等,2010. 粒度分析在中国第四纪黄土古气候研究中的应用现状[J]. 中国海洋大学学报,40(2):79-84. [51] 刘智荣,薛怀宇,王昌盛,2021. 河北平原中部保定西伯章剖面晚第四纪沉积特征及其环境意义[J]. 地质力学学报,27(6):1011-1023. doi: 10.12090/j.issn.1006-6616.2021.27.06.082 [52] 马永法,詹涛,杨业,等,2021. 东北平原东部天恒山钻孔色度变化特征与古气候意义及对松嫩古湖演化的指示[J]. 地质学报,95(11):3519-3531. doi: 10.3969/j.issn.0001-5717.2021.11.024 [53] 裘善文,王锡魁,张淑芹,等,2012. 松辽平原古大湖演变及其平原的形成[J]. 第四纪研究,32(5):1011-1021. doi: 10.3969/j.issn.1001-7410.2012.05.17 [54] 裘善文,王锡魁,MAKHINOV A N,等,2014. 中国东北平原及毗邻地区古水文网变迁研究综述[J]. 地理学报,69(11):1604-1614. doi: 10.11821/dlxb201411002 [55] 孙建华,谢远云,康春国,等,2022. 大兴安岭东麓平安镇白土山组的地层属性:对物源和沉积环境的指示[J]. 地层学杂志,46(2):196-208. [56] 孙建中,王雨灼,张庆云,1982. 松辽平原第四纪地层的划分:几种年代学方法的应用[J]. 长安大学学报(地球科学版)(2):79-91,10. [57] 孙杨,谢远云,迟云平,等,2022. 大兴安岭东麓龙江县白土山组地层特征:化学风化、沉积循环、源-汇体系和沉积环境[J]. 山地学报,40(1):14-28. doi: 10.3969/j.issn.1008-2786.2022.1.sdxb202201002 [58] 田明中,程捷,2020. 第四纪地质学与地貌学[M]. 2版. 北京:地质出版社. [59] 王恩宝,2012. 哈尔滨地区第四纪沉积序列及其环境意义[D]. 长春:吉林大学. [60] 王嘉新,谢远云,康春国,等,2020. 哈尔滨荒山岩芯重矿物特征对松花江第四纪水系演化的指示[J]. 第四纪研究,40(1):79-94. doi: 10.11928/j.issn.1001-7410.2020.01.08 [61] 王心源,吴立,张广胜,等,2008. 安徽巢湖全新世湖泊沉积物磁化率与粒度组合的变化特征及其环境意义[J]. 地理科学,28(4):548-553. doi: 10.3969/j.issn.1000-0690.2008.04.016 [62] 王永,董进,杨劲松,2020. 哈尔滨荒山剖面第四纪地层研究[J]. 地球科学,45(7):2662-2672. [63] 吴鹏,谢远云,康春国,等,2020a. 哈尔滨荒山黄土的成因:粒度、地球化学、磁化率、沉积和地貌特征的整合记录[J]. 地球学报,41(3):420-430. [64] 吴鹏,谢远云,康春国,等,2020b. 早更新世晚期松花江水系袭夺:地球化学和沉积学记录[J]. 地质学报,94(10):3144-3160. [65] 吴锡浩,浦庆余,钱方,等,1984. 松辽平原第四纪磁性地层的初步研究[J]. 海洋地质与第四纪地质,4(2):1-13. [66] 徐园园,谢远云,康春国,等,2022. 松花江早更新世水系演化:来自TIMA矿物和地球化学的证据[J]. 地质科学,57(1):190-206. doi: 10.12017/dzkx.2022.012 [67] 曾吉莉,迟云平,谢远云,等,2024. 哈尔滨三截地砂砾石堆积地层属性特征:对地层区域划分和构造的指示[J/OL]. 沉积学报,1-21[2024-11-27]. https://doi.org/10.14027/j.issn.1000-0550.2024.026. [68] 詹涛,曾方明,谢远云,等,2018. 东北平原东部天恒山钻孔的粒度特征及其对地层划分的指示[J]. 地质力学学报,24(4):515-521. doi: 10.12090/j.issn.1006-6616.2018.24.04.054 [69] 詹涛,曾方明,谢远云,等,2019. 东北平原钻孔的磁性地层定年及松嫩古湖演化[J]. 科学通报,64(11):1179-1190. [70] 张曼,谢远云,康春国,等,2021. 哈尔滨荒山岩芯黄土-古土壤的化学风化特征:对古土壤形成环境指示[J]. 土壤学报,58(3):673-684. [71] 张月馨,迟云平,谢远云,等,2020. 中更新世以来哈尔滨黄土有机碳同位素组成及其古气候意义[J]. 地球学报,41(4):525-534. doi: 10.3975/cagsb.2020.040602 [72] 中国植被编辑委员会,1980. 中国植被[M]. 北京:科学出版社. 期刊类型引用(3)
1. 孔凯，戚盛名，李鹏杰，黄书翔，崔超超，方坤. 浸润作用及动载损伤后砂岩力学性能研究. 采矿与安全工程学报. 2025(01): 40-49 . 
百度学术2. 吕远强，叶万军，赵金刚，姜海波，蒋蓓茹，张雪丽. 冻融循环作用下裂隙岩石损伤扩展演化特征. 河北大学学报(自然科学版). 2024(06): 581-592 . 百度学术3. 叶美娜，仲海书，张震，李保珠，田向盛，张英. 甘肃早子沟金矿典型岩石微观结构特征及其渗透性研究. 地质与勘探. 2024(06): 1142-1152 . 百度学术其他类型引用(0)
-
下载:
下载:
百度学术
计量
- 文章访问数: 285
- HTML全文浏览量: 55
- PDF下载量: 57
- 被引次数: 3
