河北海兴小山CK3孔磁性地层与第四纪火山活动记录
胡云壮1, 胥勤勉1, 袁桂邦1, 裴军令2, 杨吉龙1, 张玉发1, 王强1
1 中国地质调查局天津地质矿产研究所,天津 300170
2 中国地质科学院地质力学研究所,北京 100081

通讯作者简介 王强,男,1945年生,天津地质矿产研究所研究员,研究方向为第四纪地质与海洋地质。E-mail:tjwq1945@163.com

第一作者简介 胡云壮,男,1983年生,天津地质矿产研究所助理研究员,研究方向为第四纪地质与水文地质。E-mail:huyunzhuang@163.com

摘要

基于对渤海湾西岸河北省海兴小山火山丘附近 500.25 m CK3孔进行的磁性地层学研究,确认该钻孔揭示了吉尔伯特( Gilbert)极性时科奇蒂( Cochiti)亚时以来的地层,底部层位年龄约为 4.38 Ma;结合 14C 测年和地层关系,可见本孔中所见 4期火山活动分别发生在 2.12 Ma 1.34 Ma 35~70 ka 10~15 ka;与渤海湾滨海平原几个质量较高的钻孔古地磁研究成果对比显示,黄骅坳陷第四纪期间整体向 NE方向倾斜,沉积物平均粒径变化与测井曲线有较好的对比关系和同步性。依据沉积物特征、沉积构造以及测井相、粒度分析资料综合分析,确定全孔揭露地层大多发育泛滥平原沉积,即使是分支间湾亚相也少见湖泊—湖沼相沉积,沉积物多呈氧化状态。晚更新世以来 3期海侵地层皆为潮坪沉积,其底板埋深与整个渤海湾西岸所见相符,可以按海洋氧同位素分期划分。

关键词: 河北省; 海兴小山; 火山口地区; 磁性地层; 第四纪火山活动; 海侵地层; 测井相
中图分类号:P534.63 文献标志码:A 文章编号:1671-1505(2014)03-0411-16
Magnetostratigraphy of Borehole CK3 and record of the Quaternary volcanic activities in Xiaoshan of Haixing, Hebei Province
Hu Yunzhuang1, Xu Qinmian1, Yuan Guibang1, Pei Junling2, Yang Jilong1, Zhang Yufa1, Wang Qiang1
1 Tianjin Institute of Geology and Mineral Resources,China Geological Survey Bureau,Tianjin 300170
2 Institute of Geomechanics,Chinese Academy of Geological Sciences,Beijing 100081

About the corresponding author Wang Qiang,born in 1945,is a professor at Tianjin Institute of Geology and Mineral Resources.He is mainly engaged in Quaternary geology and marine geology.E-mail: tjwq1945@163.com.

About the first author Hu Yunzhuang,born in 1983,is an associate researcher at Tianjin Institute of Geology and Mineral Resources. He is mainly engaged in Quaternary geology and hydrogeology.E-mail:huyunzhuang@163.com.

Abstract

After the palaeomagenetism study on the 500.25 m deep Borehole CK3,which is located near the Xiaoshan volcanic dome of Haixing County on the western coast of Bohai Bay,it can be determined that this borehole has drilled through the strata since the Cochiti subchron of Gilbert polarity chron,with the basal horizon having an age of 4.38 Ma. Combining with14C dating data as well as stratigraphic relations,the total four periods of volcanic activities occurred at 2.12 Ma,1.34 Ma,35~70 ka and 10~15 ka,respectively. The comparison with the high-quality palaeomagenetism-studied boreholes on the coastal plain of the Bohai Bay indicated that the Huanghua Depression was overall tilted northeastwards during the Quaternary. The average grain-size change of sediments and the logging curves showed a good contrastive relationship and synchronicity. The comprehensive analysis of the sediment characteristics,sedimentary structures,log facies and particle sizes suggested that most of the drilled strata were floodplain deposits,and had rare lacustrine-limnetic sediments even in the interdistributary subfacies,while the deposits were mostly in a state of oxidation. All the three periods of transgressive strata since the Late Pleistocene were tidal flat deposits with the burial depth of basal horizons conformed to the entire western Bohai Bay,and could be subdivided based on marine oxygen isotope stages.

Key words: Hebei Province; Xiaoshan of Haixing; crater area; magnetostratigraphy; Quaternary volcanic activity; transgressive stratum; log facies

渤海湾西岸作为现代沉降海岸堆积区, 第四系厚度一般达300~400 m(李应培, 1965)。上世纪70年代中后期, 古地磁技术开始应用于钻孔地层研究, 得以确定地层年代。除去北京山前平原凹陷区, 总结上世纪末华北平原已有的古地磁资料(王强等, 2003), 结合近10年来天津、冀东地区钻孔磁性地层学研究进展(赵长荣等, 2003; 姚政权等, 2006; 肖国桥等, 2008; 施林峰等, 2009; 胡云壮等, 2014), 可以确认冀东平原滨海的南堡凹陷地区(李华梅和王俊达, 1983; 袁桂邦等, 2014), 是目前中国北方沿海平原沉降最大的地点。

与此同时, 渤海湾也是新构造运动相对活跃的地区, 除差异构造沉降外, 在部分地区第四系松散沉积物中亦发现火山岩类堆积。渤海湾盆地在新生代拉伸应力场成盆过程中, 可能有来自深层的动力过程(李三忠等, 2010)。由此, 河北海兴小山火山活动受到了关注。

本次研究在小山凝灰岩丘西南翼施工的500 m 深CK3孔, 系1︰25万黄骅幅区调项目标准孔, 揭露了4期火山堆积物, 其中一层为玄武岩; 60 m 以上见3期海侵沉积地层, 构成了沉降海岸区有玄武岩层位的洪泛— 湖泊— 火山— 海侵— 海退层序, 为磁性地层研究和区域古地理演化研究提供了理想材料。

1 区域地质背景
1.1 研究区概况

渤海湾西岸在大地构造位置上属华北坳陷区, 其中包括冀中坳陷、沧县隆起、黄骅坳陷、埕宁隆起和济阳坳陷等三级构造单元(图 1)(武汉地质学院北京研究生院石油地质岩相组和大港油田石油勘探开发研究院勘探室岩相组, 1987; 河北省地质矿产局, 1989; 大港油田地质志编辑委员会, 1991); 黄骅坳陷、埕宁隆起、济阳坳陷与渤中坳陷、下辽河坳陷又构成中新生代裂谷型的渤海湾盆地(翟光明, 1990); 进入第四纪后, 该地区仍保持差异沉降的特征。

图1 渤海湾西岸主要构造单元及CK3钻孔位置示意图(据陈望和和倪明云, 1987, 重绘)Fig.1 Map of main tectonic units and location of Borehole CK3 in western Bohai Bay(modified from Chen and Ni, 1987)

小山位于海兴县城东北3 km 远处滨海平原上, 东距渤海湾约23.5 km, 火山残丘高出地表约30 m, 山顶高程35.1 m(冯金良, 1997), 呈向南开口的椅状隆丘, 东翼、北翼高于西翼, 且出露较好, 西翼大部伏于地表之下; 东翼、北翼地表所见火山堆积主要为岩屑凝灰岩、沉凝灰岩、凝灰质砂层与黏土层等, 仅局部略见玄武岩体。结合区域钻孔地层所见, 邵时雄等(1983, 1984)对河北平原东部火山活动期次和火山岩类型、特征以及相关的新构造运动进行了总结; 但是由于没有专门钻孔的磁性地层学和其他更详细的内容, 故而小山火山未纳入《中国火山》(刘嘉麒, 1999)专著。20年前小山椅状隆丘西南翼钻孔7.3~7.4 m 凝灰岩上覆泥炭沉积获14C测年6300± 80 a BP 数据, 近期在该火山口东翼露头与钻孔剖面岩屑凝灰岩中获得37± 6.7 ka至46± 16 ka 的6个ESR测年数据(尹功明等, 2013), 从而显示火山隆丘西南翼地表下是正常沉积地层, 其东北方向1 km 余的东翼隆丘是MIS3时期火山喷发形成的。

1.2 第四纪磁性地层学研究

华北平原区东部钻孔磁性地层研究始于1970年代末河北省境内衡水钻孔工作(李华梅等, 1977)。随后发表有河北省和天津市等地钻孔研究成果(陈望和和倪明云, 1987; 张志良, 1988; 王淑芳和王云生, 1991; 王强等, 1992); 其中海兴地区古地磁学研究, 包括小山近正南方向约 10 km 的辛集沧13孔(张宏才等, 1978)以及再向西南约 12 km 的高湾7-17孔(陈望和和倪明云, 1987), 但这2个钻孔样品采样密度过低, 很难使用其结论。

与海兴地区所在埕宁隆起相邻的黄骅坳陷钻孔古地磁工作较多, 冀东柏各庄农场Ba3(柏3)孔研究(李华梅和王俊达, 1983), 显示南堡凹陷地区构造沉降较大; 随后天津汉沽沿海大神堂HG1孔测试, 确定B/M界线在200 m, M/G界线在410 m, 显示黄骅坳陷内部第四纪期间向北东方向倾斜(王强等, 1992); 近期南堡凹陷曹妃甸生态城Bg10孔研究, 确定B/M界线在162.3 m, M/G界线在477.7 m(袁桂邦等, 2014), 证实了此前认识。

近10年来沿海地区又有多个钻孔古磁性地层学研究的报道, 无论是样品采样密度还是测试质量都有了很大进展。

天津市区以南BZ2孔位于沧县隆起上, 孔深203.60 m, 56.2 m 即为0.78 Ma的B/M界线, 第四系底界M/G界线仅在164 m 深度(姚政权等, 2006); 沧东断裂以东、黄骅坳陷西斜坡上、大港区北部中塘BZ1孔, 孔深204.41 m, B/M界线为102.43 m(肖国桥等, 2008)。依据这2个钻孔构成的北西— 南东向剖面线对比, 判断其B/M 界线埋深差异, 是2个构造单元之间的沧东断裂第四纪以来活动所致(陈宇坤等, 2008)。同样位于黄骅坳陷西斜坡的大港区南部小王庄CQJ4孔, 孔深500 m, B/M界线位于74.0 m, M/G界线位于340.0 m(施林峰等, 2009)。大港区东南端歧口凹陷马棚口村西BQ1孔, 孔深95.16 m, 未钻遇B/M界线(赵长荣等, 2003)。在渤海海域渤中凹陷-25 m 水深处施工的BH08钻孔中, B/M界线位于125.66 m(Yao et al., 2014), 成为与渤海湾西岸对比的新基点。

1.3 海侵层序列研究

渤海湾西岸在新近纪晚期以来就是河湖相为主的沉积环境, 第四纪晚期发生3期海侵(李应培, 1965)。20世纪70年代, 借助河北省沧州地区与天津市大规模农田供水钻探的几十个钻孔第四纪研究, 证实了前述华北新近纪以来河湖相沉积为主的认识, 确定大规模海侵发生在间冰期, 由此提出海侵期次的划分; 嗣后, 随着工作的深入、认识的提升以及与沉积学的密切配合, 将海相微体生物标本丰富、沉积特征确切的3期海侵地层划归入晚更新世以来(王强等, 1986), 结合局部下切河谷地区综合总结, 可见其底板埋深大致在60 m 上下或60~70 m 深度内(王强和李从先, 2009), 仅冀东曹妃甸地区区域构造沉降较大、晚更新世以来海侵层埋深较大(汪品先等, 1981; 李元芳等, 1982; 胥勤勉等, 2011)。

2 CK3孔基本情况及沉积特征

河北海兴县小山区域构造位置处于埕宁隆起之上, CK3孔位于火山丘西翼南端, 孔深500.25 m, 岩心采取率为95.7%; 砂层累计厚度仅有83.9 m, 以致全孔测井曲线少有起伏, 其余多为黄棕色、棕褐色泛滥平原相黏土质粉砂、粉砂质黏土或黏土沉积, 少见棕灰色— 灰绿色湖相细粒沉积, 连续且最厚的浅灰色沉积物仅出现在108.9~110.0 m, 由此判断该孔在地质时期长期处于较高地貌部位。孔深约60 m 以内自上而下见有相当区域MIS(海洋氧同位素)1、3、5的3期海侵层, 深度分别为 3.21~13.14 m、21.20~30.89 m 和56.12~59.10 m, 以及其间的泛滥平原沉积。全孔见有4期火山活动堆积, 深度分别为13.14~16.18 m, 33.37~48.02 m, 145.80~154.05 m 及222.16~233.80 m。依据岩性特征, 结合自然电位(SP)、侧向视电阻率(LRt)和自然伽马(γ )测井曲线、平均粒度曲线反映的粒序、沉积旋回, 可将该孔从上到下大致划分为22大层:

1)0~3.21 m:以浊黄橙色、亮黄棕色粉砂为主; 上部见有大量现生芦苇根系、洞穴及氧化斑点, 系近代洼地沉积; 中、下部夹薄层状浊黄棕色黏土, 黏土中含少量铁锰质结核, 为泛滥平原沉积。

2)3.21~13.14 m:自上向下依次为浊黄棕色— 黄灰色— 灰色有机质黏土, 多见生物扰动构造, 中部(图 2-a)多见海相贝壳碎片(中华青蛤Cyclina sinensis(Gmelin)、中国不等蛤Anomia chinensis Philippi、光滑河篮蛤Potamocorbula laevis(Hinds)等), 为淤泥质潮坪沉积; 下部9.05~13.14 m 以灰色有机质黏土为主, 其中9.2~9.3 m和11.90~11.97 m为基底泥炭层(图2-b), 下伏滨海湖沼沉积。

图2 CK3孔岩心照片Fig.2 Core photographs of Borehole CK3

3)13.14~16.18 m:灰橄榄色玄武质沉凝灰岩(图 3-a), 半固结, 发育水平层理。

图3 CK3孔典型火山堆积物岩心照片Fig.3 Photographs of volcanic accumulation in Borehole CK3

4)16.18~31.89 m:上部16.18~21.20 m以黄棕色— 黄灰色黏土、黏土质粉砂、粉砂为主, 粉砂层发育水平层理, 为河流泛滥沉积, 相当末次盛冰期地层; 中部21.20~30.89 m灰色黏土中发育潜穴, 见少量海相贝壳碎片, 系高潮坪沉积(图 2-c); 下部29.80~30.89 m见夹海相贝壳碎片的砂质潮下带沉积(图 2-d), 下伏滨海湖沼沉积。

5)31.89~33.37 m:浊棕色黏土夹浊黄橙色粉砂, 发育少量铁锰质结核和钙质结核, 砂层发育小型交错层理, 为极浅水、略有坡度、地下水活动频繁的泛滥平原越岸沉积。

6)33.37~48.02 m:暗灰色玄武质角砾岩、玄武质沉凝灰岩、玄武质火山碎屑岩(图 3-b, 3-c), 其中47.78~48.02 m褐色玄武质角砾岩夹褐色黏土, 为下伏松散沉积地层的烘烤层。

7)48.02~59.74 m:上部48.02~56.12 m 以浊棕色粉砂质黏土、黏土为主, 含极少量钙质结核、铁质结核; 局部夹薄层粉砂(图 2-e), 为河间洼地或边滩沉积; 中部56.12~59.10 m 为灰黄棕色黏土与粉砂互层, 富含有机质, 多见生物扰动构造, 58.80~59.10 m 见灰色含大量海相贝壳碎片粉砂, 为砂质潮下带沉积; 59.10~59.74 m 见冲刷界面, 下伏灰棕色中营养湖泥质沉积, 见有生物潜穴(图 2-f)。

8)59.74~72.83 m:顶部为杂色黏土, 发育潜育化和钙质淀积; 上部为暗灰黄色粉砂质黏土, 较硬, 含较多钙质结核; 下部65.03~72.83 m 为浊橙色钙质粗粉砂, 含大量钙质结核; 为一期快速洪泛沉积及其上覆泛滥平原沉积。

9)72.83~104.20 m:上部72.83~80.60 m为浊棕色— 棕色黏土, 含较多钙质结核, 且较硬(硬黏土); 下部80.6~104.2 m 以棕色黏土、黏土质粉砂为主, 发育钙质结核(图 2-g, 2-h), 偶见贝壳碎片, 局部有机质富集, 黏土质粉砂普遍发育水平层理; 系泛滥平原夹河间洼地沉积。

10)104.20~118.29 m:以棕色— 棕灰色黏土、粉砂质黏土、粉砂为主; 黏土层发育少量钙核, 局部有机质富集呈灰色; 粉砂层发育水平层理、交错层理, 偶见碳质纹层; 114.8~115.2 m见黑棕色含贝壳粉砂, 所含淡水双壳类丽蚌未定种Lamprotula sp.碎片占剖开岩心表面积的20%, 钙质结核占5%~10%(图 2-i, 2-j); 为泛滥平原夹湖沼— 湖泊沉积。

11)118.29~145.80 m:上部以浊棕色黏土、粉砂质黏土为主, 局部呈浊红棕色, 普遍含钙质结核, 局部富含有机质; 下部以黄棕色黏土质粉砂为主, 局部见水平层理; 为泛滥平原夹湖沼沉积。

12)145.80~154.05 m:灰绿色玄武质沉凝灰岩、玄武质凝灰角砾岩(图 3-d)。

13)154.05~222.16 m:上部154.05~172.50 m 以浊棕色— 棕色黏土、粉砂质黏土为主, 较硬, 普遍含钙质结核, 部分为钙质胶结, 局部轻微潜育化(图 2-k); 中部172.5~176.30 m为浊棕色— 棕色粉砂, 含少量碳质纹层, 构成水平层理, 局部发育槽状交错层理; 下部176.30~222.16 m以浊黄棕色黏土、粉砂质黏土为主, 含较多小钙质结核, 普遍沿裂隙发育网纹状潜育化(图 2-l), 见少量铁锰质斑点。该层为多期分支河道夹分支间湾沉积。

14)222.16~233.80 m:为灰黑色— 黑色玄武岩夹玄武质凝灰岩(图 3-e); 玄武岩呈块状构造, 致密, 细— 微细粒构造, 见气孔和杏仁构造, 地球化学结果显示拉斑玄武岩特征(另文发表)。

15)233.80~261.85 m:顶部为黏土质烘烤层, 与上覆玄武岩界面呈不规则接触; 中上部为浊黄棕色黏土, 致密块状, 含钙质结核(图 2-m); 底部为接近5 m厚的棕色粉砂, 发育平行层理及小型槽状交错层理, 为一期水动力较强的分支河道沉积, 上覆泛滥平原沉积。

16)261.85~278.90 m:上部以浊红棕色黏土、粉砂质黏土为主, 轻度钙质胶结, 质地较硬脆, 局部发生较强潜育化, 呈灰橄榄色; 下部为4.2 m 厚棕色粉砂, 发育水平层理, 氧化作用较强。此层为一期分支河道与间湾沉积。

17)278.90~331.55 m:上部278.9~290.2 m以棕色黏土质粉砂、粉砂为主, 局部见含淡水贝壳碎屑的薄砂层(图 2-n), 发育小型交错层理、槽状层理、斜层理; 中部290.20~310.95 m为棕色黏土, 中上部皆发生不均匀潜育化而呈现灰橄榄色色斑, 含钙质结核; 下部310.95~327.95 m 棕色黏土夹粉砂层, 局部为钙质胶结砂岩(图 2-o); 底部327.95~331.55 m 为黄棕色— 棕色粉砂, 发育板状交错层理, 其余偶见小型交错层理(图 2-p)。此层为一期大分支河道转为进积河流与间湾沉积。

18)331.55~386.70 m:上部331.55~377.80 m以棕色— 灰橄榄色黏土为主, 裂隙较发育, 裂隙面光滑, 偶见钙核, 含少量易碎的锰质结核, 局部发生潜育化呈灰橄榄色(图 2-q); 下部377.8~386.7 m 以浊黄棕色粉砂、钙质粉砂为主, 发育小— 中型槽状交错层理。此层为分支间湾及浅湖沉积。

19)386.7~407.1 m:上部386.70~399.45 m为浊红棕色黏土夹粉砂, 黏土硬、脆, 多见裂隙面及不均匀潜育化, 粉砂层发育水平层理; 下部399.45~407.10 m为亮红棕色(褐红色)粉砂(图 2-r), 发育板状交错层理和斜层理, 强烈氧化成褐色— 褐红色。此层为一期大分支河道沉积发育期。

20)407.10~445.80 m:上部407.10~419.55 m为暗红、红棕色黏土, 质地较硬、脆, 局部见较大裂隙面, 并沿裂隙面发生潜育化、钙化(图 2-s); 下部419.55~445.80 m 棕色— 暗红棕色粗粉砂, 见水平层理和低角度斜层理, 偶见碳质条纹; 其中434.2 m少量钙质胶结, 因轻度潜育化而呈黄绿色。此层同样为一期大分支河道转为分支间湾沉积。

21)445.8~481.1 m:上部 445.8~458.6 m 为浊棕色粉砂质黏土; 普遍含钙质结核, 轻度不均匀潜育化致局部呈淡黄绿色; 下部458.6~481.1 m 以红棕— 棕色黏土、黏土质粉砂为主, 钙质胶结半固结, 质地较硬, 易碎, 局部沿裂隙发生潜育化。此层为极浅水的泛滥平原沉积。

22)481.10~500.25 m:上部481.1~494.4 m 为灰色钙质黏土与暗红棕色黏土互层, 局部发育较强不规则潜育化(图 2-t); 下部494.40~500.25 m 为暗红棕色粉砂质黏土, 系遭强烈潴育化所致, 含少量钙质结核。此层为泛滥平原与河间洼地交互沉积。

3 实验与结果分析
3.1 14C测年

对CK3孔与其北1187 m 处X1孔关键层位进行了14C测年(表 1)。依据CK3孔9 m余和11 m 余2个液闪仪测试泥炭样品年龄, 确定其上覆地层为中全新世海侵层。

表1 CK3与X1孔14C测年数据 Table1 14C dating data from Boreholes CK3 and X1

X1孔34.0~34.3 m 为相当基底泥炭层位的绿灰色黏土, 据其2个样品AMS14C测年, 判断上覆地层为MIS3海侵层; 在X1孔33.6 m贝屑层中, 以单瓣无破损的滨海相双壳类缩缢蛏Sinonovacula constricta(Lamarck)壳体测年, 为早于43 500 a BP, 当属于再搬运材料。依据层位对比, 判断CK3孔28.70~32.51 m 含海相贝壳碎片的灰棕色黏土夹薄层灰黄棕色粉砂层, 同样为MIS3海侵层。

3.2 古地磁取样与实验

采取古地磁样品时将岩心从中间对半切开, 以2 cm× 2 cm× 2 cm 的无磁塑料盒内装置样品, 间距在0.2~1.0 m 之间; 个别层位因岩心扰动或缺失, 采样间距放宽超过1 m, 砂层未采取样品。全孔共采得古地磁样品500块; 其中包括14块火山岩柱状样及其间湖泊夹层样品3块。此外, 采取了粒度分析样品738块。

全部样品的古地磁测试在中国地质科学院地质力学研究所古地磁实验室完成, 所有的退磁及剩磁测试仪器均置于磁屏蔽屋内, 采用热退磁和交变场退磁相结合, 以避免热退过程中造成粗粒度样品的破碎散落。对关键层位的122块样品(包括火山样品)测量天然剩磁后, 按25~50 ℃ 的间隔从室温至580 ℃ 或680 ℃、在美制TD-48高容量热退磁仪中进行系统热退磁; 对378块粒度相对稍粗的样品采用交变退磁法, 使用美国GSD-2型交变退磁仪, 按3~10 mT 步长自5 mT 到60 mT 进行。剩磁测量在美制立式2G-755R超导磁力仪上进行, 所有退磁过程和剩磁测量均在零磁空间(< 300 nT)中进行。

粒度样品测试在南京大学地理与海洋学院地表过程实验室完成。测试仪器为英国产Mastersizer 2000粒度仪。粒度仪测量范围为0.2~2000 μ m, 重复测量的相对误差小于3%, 每个样品测试时间为3min左右。

3.3 古地磁结果分析

样品的剩磁组分均利用国际上通用的Enkin古地磁软件包进行主向量分析。代表性样品的剩磁矢量正交投影图如图4所示(由于钻探过程岩心已转动, 图中所示的磁偏角无实际意义, 只作分析数据参考)。测试的部分样品剩磁方向从第2步开始就稳定趋向原点; 而大部分样品的剩磁显示出2个分量, 在200~300 ℃ 或12~15 mT 下就可洗去, 系次生粘滞剩磁分量; 第2分量在300 ℃ 或15 mT 以后保持稳定, 并逐步趋向原点, 代表了原生特征剩磁的方向。退磁结果显示, 交变退磁的样品, 一般在20~50 mT 磁场之间获得稳定的特征剩磁。也有占样品总数6.4%的32块样品在退磁过程中方向不稳定, 未能分离出其特征剩磁方向; 绝大部分热退磁样品可在300~580 ℃ 获得稳定的特征剩磁, 有部分样品(约占热退磁样品数的1/4)直到680 ℃ 剩磁才衰减到零, 并获得稳定的特征剩磁, 表明沉积物中特征剩磁的载体有磁铁矿和赤铁矿。测试过程中发现, 强氧化条件下形成的强潴育化沉积物、多含钙质结核层位的样品, 测试数据质量一般不高; 淡灰色、灰绿色沉积物样品剩磁较弱, 且难以退净, 故而最终选择468块样品编绘磁倾角曲线。

图4 CK3孔代表性样品的退磁正交矢量投影图Fig.4 Orthogonal(Zijderveld)vector plots of representative specimens from Borehole CK3

CK3钻孔地层岩性柱、磁倾角随深度的变化曲线、测井曲线和获得的磁极性柱及与标准极性柱(Cande and Kent, 1995)的对比如图5所示, 所见古地磁极性亚时、极性带转换界线与标准极性柱有很好的可比性, 自上而下确定为4个极性带:

1)0~89.15 m:为以正极性为主的布容极性带, 其间夹有1个较薄的负极性漂移。

2)89.15~275.90 m:为以负极性为主的松山极性带, 其中2段显示正极性, 122.6~140.5 m判断为哈拉米洛(Jaramillo)亚带, 177.4~218.4 m为奥尔都维(Olduvai)亚带。

3)275.9~355.1 m:正向极性时判断为高斯极性带。298.0~324.7 m存在1段负极性带, 判断为凯纳(Kaena)亚时, 而在凯纳亚时下部存在连续2个样品呈现负极性, 样品数量过少, 考虑是遭受河流侵蚀、地层缺失所致, 暂定为马莫斯极性(Mammoth)亚时。

4)355.1 m以下:为吉尔伯特(Gilbert)负极性带, 其中456.7~480.4 m为科奇蒂(Cochiti)亚带。

图5 CK3孔测井曲线、平均粒径曲线、磁倾角(I0)曲线和古地磁极性柱Fig.5 Logging curves, curves of mean grain size and magnetic dip, and paleomagnetic column of Borehole CK3

3.4 测井曲线与平均粒径曲线的匹配

用测井曲线判断沉积相需要岩心研究配合(马正, 1982; 欧阳健等, 1999; Donselaar and Schmidt, 2005; Tä navsuu-Milkeviciene and Sarg, 2012)。以自然电位、视电阻率和自然伽马(γ )测井曲线分别配以各自的镜像曲线, 可以3类或2类曲线的共同变化趋势确定粒序和沉积旋回(王强和李从先, 2009), 继而以沉积物的粒度特征用于沉积环境类型的辨识和物质运动方式的判定(Russell, 1939; Doeglas, 1968; Visher, 1969)。在同一沉积环境中, 底质沉积物粒度分布的空间变化, 反映了多种动力搬运作用的综合(McLaren, 1981; Pedreros et al., 1996; 孙有斌等, 2001)。

据CK3孔的野外描述与岩心照片及详细粒度分析统计, 全孔厚度大于1 m 以上的砂层为19层, 与视电阻率、自然电位反映的砂层吻合率分别高达 90% 与85% 左右, 产生的误差基本是由于岩心缺失或钻探磨损2种因素产生的深度误差所致。火山岩层由于成分、结构、构造不同, 具有不同类型的测井响应特征(石强, 1996)。影响测井伽马值唯一的因素是放射性核素(U、Th和40K)的含量, 已经发现酸性凝灰岩的自然伽马曲线有高异常响应, 但并非是泥岩, 故而依据伽马值计算泥质含量会产生较大误差; 各类火山事件黏土岩的伽马值相差悬殊, 但其中泥质含量几乎相等(冯宝华, 2008)。上述现象在CK3孔尚未出现。该孔222.16~233.80 m 所见玄武岩岩性致密, 裂隙较少, 测井曲线上具有低自然伽马、高电阻率的特征; 145.80~154.05 m对应低伽马曲线层段, 自然电位、视电阻率曲线几乎没有反应; 40 m余与10余米2段火山岩层段的3类测井曲线匹配关系更是复杂, 对此留待进一步探讨。

平均粒径是累积曲线上与累积百分含量分别为16%、50%、84%相对应的3个粒径的平均值(Folk and Word, 1957), 包括了样品全部粒径的2/3以上的成分。图5所示CK3孔剖面中记述的相对粗颗粒层位, 平均粒径曲线与之基本完全对应, 较好地反映了第四纪松散沉积物的岩性; 同时, 平均粒径曲线变化亦较好地展示了钻孔地层沉积物粒序变化的趋势。

利用岩心与自然电位、视电阻率及自然伽马测井曲线3条测井曲线相结合, 深度标定, 加之平均粒径曲线的印证, 使钻孔沉积旋回划分更加准确。CK3孔反映出的沉积环境主要为分支河道、浅水湖沼及泛滥平原等。

4 讨论
4.1 海侵地层与火山堆积层年代

海侵层的基底泥炭及其测年是确定海侵发生年代的重要依据。渤海湾西岸第2海侵层底部基底泥炭的14C年龄大多为30~36 ka, 稍高地貌部位年龄更年轻些, 对应MIS3中、晚期(施雅风和于革, 2003; 王强等, 2008)。前述3.1小节已经总结CK3孔11.97 m 以上为中全新世海侵层, 28.70~32.51 m 为近40 ka 以来的MIS3海侵层。晚更新世末次间冰期(赵松龄等, 1978)相当MIS5时期(王强和金权, 1989; 王强和李从先, 2009), 晚更新世早期海侵层是以100~128 ka 的古地磁布莱克(Blake)亚时为底界(赵松龄等, 1978), 但是多年来并未得到更多的证实, 这一现象可能与滨海平原海侵发生时易遭受侵蚀有关。按照沉积特征, 判断56.12~59.74 m 为MIS5海侵层。目前看, 除河北南堡地区以外的海侵层埋深较大外, CK3孔3期海侵层埋深与渤海湾西岸的规律(王强等, 1986)相符。

根据钻孔顶部的14C测年、沉积特征, 判断CK3孔全新统底界在13.0 m; B/M界线、即中更新统底界为89.15 m, M/G界线、即第四系底界为275.9 m。根据近钻孔底部科奇蒂亚时的沉积速率外推, 大致得出底部年龄约为4.38 Ma。

根据磁性地层研究结果, CK3孔所见4期火山堆积中, 第1期火山层位(222.16~233.80 m)对应奥尔都维极性亚时下界, 据沉积速率推算, 其发生年代约为2.12 Ma; 第2期火山堆积层(145.80~154.05 m)位于哈拉米洛亚时下界, 其发生年代约为1.34 Ma。根据钻孔沉积特征可知, 第3期火山沉积(33.37~48.02 m)处于区域第2与第3海侵层之间, 第4期火山(13.14~16.18 m)处于第1与第2海侵层之间, 下伏泛滥平原沉积相当末次盛冰期地层, 依据前述X1孔相关层位14C测年, 以及火山隆丘东侧钻孔样品ESR测年数据37± 6.7 ka(尹功明等, 2013), 其发生年代分别对应于35~70 ka 和10~15 ka 之间。取230 m处玄武岩样品进行氩— 氩绝对年龄测试, 获得14 Ma 数据。鉴于该孔与区域钻孔磁性地层学研究有较高的可信度, 而被动大陆边缘多来自软流层的玄武岩, 有可能由于老物质的再熔融混入, 会使得绝对年龄数据出入较大, 故未采纳该测年数据。

4.2 CK3孔沉积速率变化

依据与Cande 和 Kent(1995)提出的标准古地磁极性柱对比确定的地层年代, 编绘了CK3孔沉积速率— 年龄图(图 6), 可见3.58~4.38 Ma 对应深度355.10~500.23 m, 此层段平均沉积速率约为186 m/Ma。从岩心沉积特征及测井曲线(图 5)配合分析可知, 此层段为全孔分支河道最发育处, 砂层沉积颗粒相对较粗, 为快速堆积期, 其间的细粒沉积物为分支间湾及浅水湖沼沉积。2.58~3.58 Ma 对应275.9~355.1 m, 平均沉积速率约79.2 m/Ma, 沉积环境以分支间湾、浅水湖沼为主, 沉积颗粒相对较细。2.58 Ma 以来, 沉积速率约为106.9 m/Ma, 主要为泛滥平原的分支河道与间湾亚相, 局部最多为浅水湖沼, 仅在晚期出现海侵沉积; 全孔所见4期火山喷发为快速堆积。

图6 CK3孔沉积速率— 年代图Fig.6 Sedimentation rate as a function of magnetostratigraphic ages of Borehole CK3

研究结果显示, 本孔在4.38~3.58 Ma 期间沉积速率相对较大; 3.58~2.58 Ma 沉积速率相对较小; 2.58 Ma 以来沉积速率介于前2个时段速率之间, 只在1.95 Ma 和1.07 Ma 前后相对较大; 与滨海地带周边钻孔对比(图 7), 可见埕宁隆起在高斯极性时处于相对较高位置; 此前的3.58~4.18 Ma 期间, 该地处于快速沉降期。

图7 海兴小山CK3孔与河北、山东、天津钻孔磁性地层对比Fig.7 Comparison of magnetostratigraphy of Borehole CK3 in Xiaoshan of Haixing with other boreholes in Hebei, Shandong and Tianjin

在CK3孔所在埕宁隆起以南的济阳坳陷中, 古地磁样品密度最高是山东省东营辛2-4孔(文中将其编号为DX2-4), 曾解释哥德堡(Gö thenburg)和布莱克亚时的界线分别位于井深28.0 m和161.0 m, 242.5 m为B/M界线, 274.5 m 为 M/G 界线, 452 m为Ga/Gi(Gauss/Gilbert)界线, 下伏吉尔伯特极性带(周墨清等, 1992)。

在古地磁研究早期, 互相参考各地钻孔地层极性柱划分是必然的, 可以对比的特征地层是上新世棕红色黏土层以及晚更新世以来的海侵层。依据胜利油田地质研究院古生物实验室分析的该孔微体生物标本, 确定61 m 以内地层见3期海侵层, 符合整个渤海湾西岸海侵层的总结(王强等, 1986)。按照这样的海侵层序对比, DX2-4孔60 m 处由4个样品组成负极性漂移应解释为布莱克亚时, 即晚更新世开始的层位; 由此, 该孔B/M界线相应亦需上移, 解释为在118 m(图 7), M/G和Ga/Gi界线仍维持274.5 m 和452 m 的原研究意见。

依据天津市塘沽市区北部G2孔磁性地层研究结果(肖国强等, 2014), 修订天津汉沽大神堂HG1孔的极性柱解释, B/M界线上提到141 m, M/G和Ga/Gi界线修订为300 m 和479 m(图 7)。

CK3孔科奇蒂亚时结束到高斯极性时开始这一时段为快速堆积期, 厚层棕红色粉砂沉积系暴露于大气下的氧化背景中(Yao et al., 2012, 2014), 在测井曲线上对应视电阻率高幅曲线, 与钻探现场编录记载的砂层深度几乎完全吻合。

4.3 特征颜色层段对比

在不可能逐孔进行古地磁测试研究的情况下, 井下地层对比必然寻求更直观的地质特征做时代判断, 沉积物颜色段或区域稳定的岩石地层段特征即是标志。

最早判断华北平原进入第四系是以黄褐色沉积地层为主(李应培, 1965)。在取得一些钻孔古地磁测试研究结论后, 提出奥尔都维亚带以下为棕红色沉积, 再向下甚至出现紫红色、猪肝色沉积(杨子赓等, 1979; 邵时雄等, 1983; 陈望和和倪明云, 1987), 石油系统也提出过这样的总结(大港油田地质志编辑委员会, 1991; 单怀广和张慧娟, 1990; 郝治国等, 1998)。对此现象往往考虑对比的是华北山区“ 三趾马红土” 、即上新世到中新世的保德— 静乐红土的特征颜色(王强和金权, 1989)。其后, 提出河北平原宏观“ 岩石” 界线年代, 即棕色泥岩段顶板为1.8 Ma, 下伏紫色泥岩段顶板为3.15 Ma, 且以前者为第四系下限(邵时雄等, 1984)。

CK3孔下部地层已经发现棕红色甚至紫红色的沉积物(图 8), 本质上应该说是强风化背景下沉积, 反映了前第四纪的湿热气候状况。但是磁性地层学研究的开展, 使得简单依据颜色对比地层时代的方法受到挑战。首先是冀中、冀南钻孔地层分别在165 m和125 m即出现这样的棕红色沉积; 其次, 保德阶已经确定终结在吉尔伯特极性时之前约5.35 Ma(岳乐平等, 2004); 再次之, 天津塘沽地区深1223 m 全取心深钻研究显示, 300 m 以下地层几乎没有多厚的棕红色黏土出现, 而是多见厚层灰色砂层(肖国强等, 2014); 显然, 沉降较快的低洼部位易出现积水环境, 即会出现厚层灰色砂层, 而古地貌稍高部位即会多氧化状态下沉积。

图8 CK3孔下部强风化地层岩心照片Fig.8 Core photographs of heavily weathered deposits at the lower part of Borehole CK3

4.4 对古地磁极性柱解释的讨论

在确认样品密度已经足够的情况下, 对CK3孔极性柱的解释亦有过其他考虑。如80 m上下6个样品显示的负极性漂移, 是否是个负极性亚时?类似现象已经在山西运城盐湖钻孔发现, 曾笼统地对应Cande 和Kent(1995)建立的极性柱中0.493~0.504 Ma 的C1n-1亚时(王强等, 2000), 对CK3孔所见则未作进一步分辨; 同样的情况见于曹妃甸Bg10孔113~115 m, 标注为e2。哈拉米洛和奥尔都维亚时之间报道过有5次极性漂移和亚时(Valet et al., 2008), 中国黄土地区(岳乐平和薛祥煦, 1996)及山西运城盐湖湖区P3钻孔(王强等, 2000)亦报道过有后哈拉米洛(Post-Jaramillo)— 蓝田(Lantian)亚时。CK3孔111~114 m 之间5个样品构成的正极性漂移, 作者存疑, 故暂以阴影标注, 未作确定。237~244 m 的正极性漂移, 是否相当留尼旺(Ré union)亚时?鉴于该亚时历时仅10 ka(Cande and Kent, 1995), 同样以阴影标注, 以示存疑。CK3孔中马莫斯亚时是依据2个样品确定, 系考虑采自薄层分支间湾沉积层, 可能受到分支河道冲刷、发生剥蚀、侵蚀而致。

显然, 早期没有热退磁仪时期的工作并非皆不可取, 如图7中HG1和Ba3孔, 虽然已经发表了二、三十年, 但由近年附近钻孔的工作成果看, 依然得到肯定。实际早期很多钻孔古地磁样品采集密度已经达到每米3~4个样, 与文中交代舍弃测试效果不好的样品数量不同的是, 当年的报道并未叙及该问题, 也造成现在使用其研究成果的疑惑。国内外很多极性亚时的报道需要更多的交流, 如Cande和Kent(1995)提出的极性年表布容极性时里就没有哥德堡、布莱克亚时, 而是仅有一个0.493~0.504 Ma 的C1n-1亚时, 但是珠江口万顷沙钻孔地层研究中确认有哥德堡和拉尚亚时(王建华等, 2009)。

Singer(2014)依据2.6 Ma以来全球10个负极性倒转和27个极性漂移的火山熔岩流 40Ar/39Ar测年数据, 以及海洋沉积中极性漂移的天文年代计算, 建立了不稳定极性年表(GITS, Geomagnetic Instability Time Scale), 依然肯定了Cande和Kent(1995)没有提及的诸如Cobb Mountain、Blake、Laschamp等中国学者使用过的极性亚时。但中国北方沿海平原钻孔地层少见玄武岩, 难以获得充分的绝对测年材料。显然松散沉积物磁性地层学工作尚有待进一步深入。同时, 亦应该看到, 陆相地层比海洋沉积存在着更多的地层间断(王强和金权, 1989), 在测试质量得到保证的情况下, 并非所有报道的极性漂移或亚时皆可能在一个钻孔完全得到证实, 而尚需更多的工作验证。

5 结论

河北海兴小山500.25 m深CK3孔的磁性地层研究发现, 该孔钻遇了布容正极性时、松山负极性时和高斯正极性时, 以及吉尔伯特负极性时的科奇蒂极性亚时, 钻孔底部地层年龄约为4.38 Ma。

根据钻孔顶部的14C测年、钻孔地层沉积特征及区域对比判断, 全新统底界在13.0 m; B/M界线、即中更新统底界为89.15 m, M/G界线、即第四系底界为275.9 m。第四纪时期, 黄骅坳陷呈现向北东方向倾斜。

CK3孔揭露了4期火山堆积, 根据磁性地层年龄结果换算, 第1期发生年代为2.12 Ma前后; 第2期发生在1.34 Ma 前后; 依据14C和ESR测年、地层关系, 第3期与第4期发生年代大致在35~70 ka 和10~15 ka, 与此前认识相比在第四纪火山活动年代上获得了较好的时限。

华北平原包括滨海平原地层, 皆需要更多高质量的古地磁研究, 方可提高第四纪古地理研究精度和水平。

致谢 古地磁测试得到了中国地质科学院地质力学研究所古地磁实验室的大力支持和帮助, 两位评审专家提出了中肯的批评、修改意见, 谨致谢忱!

参考文献
1 陈望和, 倪明云. 1987. 河北省第四纪地质[M]. 北京: 地质出版社, 1-186. [文内引用:3]
2 陈宇坤, 李振海, 邵永新, . 2008. 天津地区第四纪年代地层剖面研究[J]. 地震地质, 30(2): 383-398. [文内引用:1]
3 大港油田地质志编辑委员会. 1991. 中国石油地质志, 卷4, 大港油田[M]. 北京: 石油工业出版社, 74-113. [文内引用:3]
4 冯宝华. 2008. 火山沉积岩及其矿床研究中伽马测井曲线解释的误区[J]. 物探与化探, 32(2): 171-174. [文内引用:1]
5 冯金良. 1997. 小山风积物特征及其形成时代探讨[J]. 地理学与国土研究, 13(2): 53-56. [文内引用:1]
6 郝治国, 郑国光, 崔占堂, . 1998. 冀中油气区第三纪地层与微体古生物群[M]. 北京: 科学出版社, 6-66. [文内引用:1]
7 河北省地质矿产局. 1989. 河北省北京市天津市区域地质志[M]. 北京: 地质出版社, 590-616. [文内引用:3]
8 胡云壮, 张金起, 白耀楠, . 2014. 3. 45 Ma以来滦河冲积扇中部唐山TD1孔记录的区域构造和气候演化[J]. 古地理学报, 16(2): 249-262. [文内引用:1]
9 李华梅, 王俊达. 1983. 渤海湾北岸平原钻孔岩心的古地磁研究[J]. 地球化学, 12(3): 196-204. [文内引用:1]
10 李华梅, 王俊达, 安芷生. 1977. 河北衡水地区钻孔岩心的古地磁研究[J]. 地球化学, 6(3): 210-218. [文内引用:1]
11 李三忠, 索艳慧, 戴黎明, . 2010. 渤海湾盆地形成与华北克拉通的破坏[J]. 地学前缘, 17(4): 64-89. [文内引用:3]
12 李元芳, 高善明, 安凤桐. 1982. 滦河三角洲地区第四纪海相地层及其古地理意义的初步研究[J]. 海洋与湖沼, 13(5): 433-439. [文内引用:1]
13 李应培. 1965. 对华北平原第四系的初步认识[J]. 第四纪研究, 4(1): 173-180. [文内引用:1]
14 刘嘉麒. 1999. 中国火山[M]. 北京: 科学出版社, 1-166. [文内引用:1]
15 马正. 1982. 应用自然电位测井曲线解释沉积环境[J]. 石油与天然气地质, 3(1): 25-40. [文内引用:1]
16 欧阳健, 王贵文, 吴继余, . 1999. 测井地质分析与油气藏定量评价[A]. 北京: 石油工业出版社, 6-65. [文内引用:1]
17 单怀广, 张慧娟. 1990. 山东油气区早第三纪介形类[A]. 见: 《中国油气区地层古生物编辑委员会》编: 中国油气区地层古生物(二), (山东油气区专辑)[M]. 北京: 石油工业出版社, 1-56. [文内引用:1]
18 邵时雄, 安仲元, 韩书华. 1984. 河北平原新构造运动主要特征分析[J]. 海洋地质与第四纪地质, 4(4): 67-77. [文内引用:2]
19 邵时雄, 张玉芳, 韩书华. 1983. 河北平原第四纪火山堆积及火山活动分期特征[J]. 海洋地质与第四纪地质, 3(2): 87-94. [文内引用:1]
20 施林峰, 翟子梅, 王强, . 2009. 从天津CQJ4孔探讨中国东部海侵层的年代问题[J]. 地质论评, 55(3): 375-384. [文内引用:1]
21 石强. 1996. 火成岩测井解释方法研究[J]. 中国海上油气, 10(6): 402-406. [文内引用:1]
22 施雅风, 于革. 2003. 40-30 ka B. P. 中国暖湿气候和海侵的特征与成因探讨[J]. 第四纪研究, 23(1): 1-11. [文内引用:1]
23 孙有斌, 高抒, 鹿化煜. 2001. 前处理方法对北黄海沉积物粒度的影响[J]. 海洋与湖沼, 32(6): 665-671. [文内引用:1]
24 王建华, 曹玲珑, 王晓静, . 2009. 珠江三角洲万顷沙地区晚第四纪沉积相与古环境演变[J]. 海洋地质与第四纪地质, 29(6): 35-41. [文内引用:3]
25 汪品先, 闵秋宝, 卞云华, . 1981. 我国东部第四纪海侵地层的初步研究[J]. 地质学报, 55(1): 1-13. [文内引用:3]
26 王强, 金权. 1989. 中国东部第四纪间断地层学初步研究[A]. 见: 杨子赓, 林和茂主编. 中国近海及沿海地区第四纪进程与事件[M]. 北京: 海洋出版社, 132-142. [文内引用:2]
27 王强, 李从先. 2009. 中国东部沿海平原第四系层序类型[J]. 海洋地质与第四纪地质, 29(4): 39-51. [文内引用:3]
28 王强, 李彩光, 田国强, . 2000. 7. 1 Ma 以来运城盆地地表系统巨变与盐湖形成的构造学背景[J]. 中国科学(D辑): 地球科学, 30(4): 420-428. [文内引用:1]
29 王强, 李凤林, 李玉德, . 1986. 对渤海西、南岸平原第四纪海侵命名的讨论[J]. 海洋学报, 8(1): 72-82. [文内引用:2]
30 王强, 刘立军, 徐海振, , 2003. 华北平原第四系下限的再研究[J]. 地质调查与研究, 26(1): 52-60. [文内引用:1]
31 王强, 宋大伟, 王云生. 1992. 第四系[A]. 见: 邓永高主编. 中华人民共和国地质矿产部地质专报, 一, 区域地质, 第29号, 天津市区域地质志[M]. 北京: 地质出版社, 116-156. [文内引用:1]
32 王强, 袁桂邦, 张玉发, . 2008. MIS3阶段以来河北黄骅北部地区海侵与气候期对比[J]. 第四纪研究, 28(1): 79-95. [文内引用:2]
33 王淑芳, 王云生. 1988. 天津迎宾馆钻孔磁性地层的划分与对比[J]. 水文地质与工程地质, (1): 41-45. [文内引用:3]
34 武汉地质学院北京研究生院石油地质岩相组, 大港油田石油勘探开发研究院勘探室岩相组. 1987. 黄骅坳陷第三系沉积相及沉积环境[M]. 北京: 地质出版社, 1-180. [文内引用:2]
35 肖国强, 杨吉龙, 赵长荣, . 2014. 天津滨海地区G2孔磁性地层年代及其构造指示[J]. 地质通报(印刷中). [文内引用:1]
36 肖国桥, 郭正堂, 陈宇坤, . 2008. 渤海湾西岸BZ1钻孔的磁性地层学研究[J]. 第四纪研究, 28(5): 909-916. [文内引用:2]
37 胥勤勉, 袁桂邦, 张金起, . 2011. 渤海湾沿岸晚第四纪地层划分及地质意义[J]. 地质学报, 85(8): 1352-1367. [文内引用:2]
38 杨子赓, 李幼军, 丁秋玲, . 1979. 试论河北平原东部第四纪几个基本问题[J]. 地质学报, 53(4): 263-279. [文内引用:2]
39 姚政权, 郭正堂, 陈宇坤, . 2006. 渤海湾海陆交互相沉积的磁性地层学[J]. 海洋地质与第四纪地质, 26(1): 9-15. [文内引用:1]
40 尹功明, 赵波, 许建东, . 2013. 河北沧州小山火山的ESR年代学研究[J]. 岩石学报, 29(12): 4415-4420. [文内引用:1]
41 袁桂邦, 胥勤勉, 王艳, . 2014. 渤海湾北岸Bg10孔磁性地层研究及其地质意义[J]. 地质学报, 88(2): 1-13. [文内引用:1]
42 岳乐平, 薛祥煦. 1996. 中国黄土古地磁学[M]. 北京: 地质出版社, 78-88. [文内引用:1]
43 岳乐平, 邓涛, 张云翔, . 2004. 保德阶层型剖面磁性地层学研究[J]. 地层学杂志, 28(1): 48-63. [文内引用:2]
44 张宏才, 黄庆福, 赵松龄, . 1978. 中国北方沿海两个沉积岩心的古地磁研究[J]. 海洋与湖沼, 9(2): 183-192. [文内引用:2]
45 张志良. 1988. 河北平原3个钻孔岩心古地磁测量及其第四纪地层年龄的确定[J]. 天津地质学会志, 6(2): 91-99. [文内引用:1]
46 翟光明主编. 1990. 中国石油地质志, 卷16, 沿海大陆架及毗邻海域油气区[M]. 北京: 石油工业出版社, 75-122. [文内引用:3]
47 赵长荣, Hus J, 阎玉忠, . 2003. 渤海湾西岸湾顶晚更新世—全新世年代地层序列与地磁极漂移[J]. 地质调查与研究, 26(3): 183-192. [文内引用:1]
48 赵松龄, 杨光复, 苍树溪, . 1978. 关于渤海湾西岸海相地层与海岸线问题[J]. 海洋与湖沼, 9(1): 15-24. [文内引用:1]
49 周墨清, 葛宗诗, 郭青霞, . 1992. 东营市东辛2-4井磁性地层划分对比研究[J]. 海洋地质与第四纪地质, 12(3): 57-69. [文内引用:1]
50 Cand e S C, Kent D V. 1995. Revised calibration of the geomagnetic polarity timescale for the Late Cretaceous and Cenozoic[J]. Journal of Geophysical Research, 100: 6093-6095. [文内引用:1]
51 Doeglas D T. 1968. Grain-size indices, classifications and environment[J]. Sedimentology, 10: 83-100. [文内引用:1]
52 Donselaar M E, Schmidt J M. 2005. Integration of outcrop and borehole image logs for high-resolution facies interpretation: Example from a fluvial fan in the Ebro Basin, Spain[J]. Sedimentology, 52: 1021-1042. [文内引用:1]
53 Folk R L, Word W C. 1957. Brazos river bar: A study in the signification of grain size parameters[J]. Journal of Sedimentary Petrology, 27: 3-27. [文内引用:1]
54 McLaren P A. 1981. A interpretation of trends in grain size measurements[J]. Journal of Sedimentary Petrology, 51: 611-624. [文内引用:1]
55 Pedreros R, Howa H L, Michel D, et al. 1996. Application of grain size trend analysis for determination of sediment transport pathways in intertidal areas[J]. Marine Geology, 135: 35-49. [文内引用:1]
56 Russel R D. 1939. Effects of transportation of sedimentary particles[A]. In: Trask P D(ed). Recent Marine Sediments[C]. The Society of Economic Paleontologists and Mineralogists, Tulsa(Oklahoma), 32-47. [文内引用:2]
57 Singer B S. 2014. A Quaternary geomagnetic instability time scale[J]. Quaternary Geochronology, 21: 29-52. [文内引用:1]
58 Tänavsuu-Milkeviciene K, Sarg J F. 2012. Evolution of an organic-rich lake basin-stratigraphy, climate and tectonics: Piceance Creek basin, Eocene Green River Formation[J]. Sedimentology, 59: 1735-1768. [文内引用:2]
59 Valet J, Plenier G, Herrero B. 2008. Geomagnetic excursions reflect aborted polarity strate[J]. Earth and Planetary Science Letters, 274: 472-478. [文内引用:1]
60 Visher G S. 1969. Grain size distribution and depositional processes[J]. Journal of Sedimentary Petrology, 39: 1074-1106. [文内引用:1]
61 Yao Zhengquan, Guo Zhengtang, Xiao Guoqiao, et al. 2012. Sedimentary history of the western Bohai coastal plain since the late Pliocene: Implications on tectonic, climatic and sea-level changes[J]. Journal of Asian Earth Science, 54-55: 192-202. [文内引用:1]
62 Yao Zhengquan, Shi Xuefa, Liu Qingsong, et al. 2014. Paleomagnetic and astronomical dating of sediment core BH08 from the Bohai Sea, China: Implications for glacial-interglacial sedimentation[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 393: 90-101. [文内引用:1]