CN113065236A - 一种连续沉积物古气候代用指标他生旋回变化的分析方法 - Google Patents

一种连续沉积物古气候代用指标他生旋回变化的分析方法 Download PDF

Info

Publication number
CN113065236A
CN113065236A CN202110294450.6A CN202110294450A CN113065236A CN 113065236 A CN113065236 A CN 113065236A CN 202110294450 A CN202110294450 A CN 202110294450A CN 113065236 A CN113065236 A CN 113065236A
Authority
CN
China
Prior art keywords
paleoclimate
convolution
change
index
depth
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110294450.6A
Other languages
English (en)
Other versions
CN113065236B (zh
Inventor
张新荣
方石
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jilin University
Original Assignee
Jilin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jilin University filed Critical Jilin University
Priority to CN202110294450.6A priority Critical patent/CN113065236B/zh
Publication of CN113065236A publication Critical patent/CN113065236A/zh
Application granted granted Critical
Publication of CN113065236B publication Critical patent/CN113065236B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Complex Calculations (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)

Abstract

本发明公开了一种连续沉积物古气候代用指标他生旋回变化的分析方法,该方法有效描述他生旋回信息特征,进而分析沉积剖面形成过程中经历的古气候背景,为深入分析古气候区域对比中存在的差异,促进全球变化研究提供基础技术支撑。本发明清楚地分析剖面上古气候相对变化的过程中的幅度和程度;相对简单、明了地将多个他生旋回分离开,对分析古气候等外界因素的影响和控制机制有很重要的作用。本发明解决现在技术所存在的因古气候代用指标的综合信息中他生旋回和自生旋回的叠加,强烈干扰古气候信息的有效性,影响古气候重建结果分析的准确性和全球对比研究进行的问题。

Description

一种连续沉积物古气候代用指标他生旋回变化的分析方法
技术领域
本发明涉及地质学、沉积学及古气候全球变化技术领域,特别涉及,一种连续沉积物古气候代用指标他生旋回变化的分析方法。
背景技术
全球气候变化是国际环境外交和国家可持续发展的重大问题之一。地球表层的沉积物,如黄土、古土壤、深海沉积、冰芯、湖泊、泥炭、石笋等蕴含着丰富的气候、环境等全球变化信息,在重建过去全球变化历史、预测未来全球变化趋势中起着至关重要的作用。
古气候是过去全球变化的一个分支。古气候信息的一个主要来源是沉积物中的古气候代用指标,如:岩石矿物指标、生物化石指标(孢粉、植硅体、有孔虫、介形虫、硅藻等)、地球化学指标(常量、微量元素、同位素等)、地球物理指标(磁化率、粒度、色度等)、腐殖化度等。这些代用指标记录了其形成时的古气候、古环境信息,与其他组分一起沉积下来保存在沉积地层中。但地层中保留的是古气候代用指标的综合信息,既含有其受气候等外界因素干扰形成的他生旋回信息,也含有其自生旋回信息。其中,他生旋回对衡量当地的古气候特征,重建古气候过程至关重要,分析他生旋回的变化过程,在古气候重建和对比研究中非常关键。
针对上述问题,设计一种分析方法,解决因古气候代用指标的综合信息中他生旋回与自生旋回的叠加,强烈干扰古气候信息的有效性,影响古气候重建结果分析的准确性和全球对比研究进行的问题。
发明内容
针对上述缺陷,本发明解决的技术问题在于,提供一种连续沉积物古气候代用指标他生旋回变化的分析方法,以解决现在技术所存在的因古气候代用指标的综合信息中他生旋回和自生旋回的叠加,强烈干扰古气候信息的有效性,影响古气候重建结果分析的准确性和全球对比研究进行的问题。
本发明提供了一种连续沉积物古气候代用指标他生旋回变化的分析方法,包括:
步骤1、确定用于分析的古气候代用指标,获取沉积剖面上不同深度采样点沉积物的古气候代用指标参数ai及其对应深度值di
步骤2、选择典型深度的沉积物进行测年,获取其绝对年龄校正值,并通过线性插值法获取剖面上不同深度采样点沉积物的年龄,估计沉积旋回的时间尺度;
步骤3、根据古气候代用指标参数ai及其对应深度值di,获取沉积物古气候代用指标在剖面上的旋回cj并确定典型旋回c0及参数代表值a0和旋回深度跨度b0
步骤4、基于典型旋回c0及古气候代用指标参数a0和深度跨度b0,获取各旋回cj的古气候代用指标参数的相对变化幅度
Figure BDA0002983792710000021
和各旋回cj的深度跨度的相对变化幅度ψj
步骤5、获取各旋回cj的古气候代用指标的古年均温变率
Figure BDA0002983792710000022
Figure BDA0002983792710000023
Figure BDA0002983792710000024
为横坐标,以ψj=(dj-d0)/d0为纵坐标,绘制各旋回c0的古气候代用指标变率幅度图,即
Figure BDA0002983792710000025
图,其中aj为各旋回cj的参数代表值,d0为旋回c0的深度跨度,dj为各旋回的深度跨度;
步骤6、以深度为纵坐标,各旋回的古气候代用指标的古年均温变率kj为横坐标,绘制kj在剖面上的变化曲线P,并得到各旋回cj的古气候代用指标变化散点图和各旋回cj的古气候代用指标变率曲线图;
步骤7、分析各旋回cj的古气候代用指标散点图、变率曲线图上各旋回的相似性、相异性,获取剖面沉积形成过程中他生旋回系统的变化过程,得到剖面上蕴含的古气候演变过程。
优选地,所述步骤3的具体步骤包括:
步骤3.1、根据古气候代用指标参数ai及其对应深度值di,获取剖面上沉积物古气候代用指标参数ai随深度di的变化曲线C;
步骤3.2、划分沉积物古气候代用指标在剖面上的旋回cj及各旋回的曲线形态类型;
步骤3.3、通过分析各旋回曲线形态特征,确定剖面上沉积物古气候代用指标参数变化曲线的典型旋回c0及典型旋回c0的参数代表值a0和旋回深度跨度b0
优选地,所述步骤4的具体步骤包括:
步骤4.1、获取各旋回cj的古气候代用指标参数ai的相对变化幅度
Figure BDA0002983792710000036
并得到
Figure BDA0002983792710000037
在剖面上的变化曲线A;
步骤4.2、获取各旋回cj的深度跨度bj的相对变化幅度ψj,并得到ψj在剖面上的变化曲线B。
优选地,所述步骤4.1的具体步骤包括:
步骤4.1.1、基于典型旋回c0及古气候代用指标参数a0和深度跨度b0,确定古气候代用指标参数ai在各旋回cj上的变化幅度
Figure BDA0002983792710000031
其中各旋回cj的参数代表值aj选择旋回cj的底部、中部或顶部样品的古气候参数值;
步骤4.1.2、确定沉积物古气候代用指标变化幅度
Figure BDA0002983792710000032
其中
Figure BDA0002983792710000033
步骤4.1.3、以深度为纵坐标,各旋回古气候参数的相对变化值
Figure BDA0002983792710000034
为横坐标,绘制各旋回cj的古气候代用指标参数相对变化幅度
Figure BDA0002983792710000035
在剖面上的变化曲线A。
优选地,所述步骤4.2的具体步骤包括:
步骤4.2.1、基于典型旋回c0及古气候代用指标参数a0和深度跨度b0,确定各旋回深度跨度bj,其中bj为旋回底部深度与顶部深度的差值,即bj=dj底-dj顶
步骤4.2.2、基于各旋回深度跨度bj,计算剖面上沉积物古气候代用指标各旋回cj的深度跨度bj相对典型旋回c0的深度跨度的变化值ψj,得到ψj=(dj-d0)/d0,其中d0为旋回c0的深度跨度,dj为各旋回的深度跨度;
步骤4.2.3、以深度为纵坐标,各旋回cj的深度跨度bj相对典型旋回c0的变化值ψj为横坐标,绘制各旋回深度跨度相对变化值在剖面上的变化曲线B。
由上述方案可知,本发明提供的一种连续沉积物古气候代用指标他生旋回变化的分析方法,是一种研究泥炭等连续沉积物沉积过程中受气候等外界因素干扰而产生的他生旋回变化的技术方法,该方法有效描述他生旋回信息特征,进而分析沉积剖面形成过程中经历的古气候背景,为深入分析古气候区域对比中存在的差异,促进全球变化研究提供基础技术支撑。本发明清楚地分析剖面上古气候相对变化的过程中的幅度和程度;相对简单、明了地将多个他生旋回分离开,对分析古气候等外界因素的影响因素和控制机制有很重要的作用,有助于快速进行古气候重建结果的解释,促进全球对比工作的有效开展。本发明作用效果显著,适于广泛推广。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种连续沉积物古气候代用指标他生旋回变化的分析方法的过程框图一;
图2为本发明实施例提供的一种连续沉积物古气候代用指标他生旋回变化的分析方法的过程框图二;
图3为本发明实施例提供的一种连续沉积物古气候代用指标他生旋回变化的分析方法的吉祥泥炭剖面植硅体反演的古年均温度变化曲线C;
图4为本发明实施例提供的一种连续沉积物古气候代用指标他生旋回变化的分析方法的吉祥泥炭剖面植硅体反演的古年均温度变化旋回划分;
图5为本发明实施例提供的一种连续沉积物古气候代用指标他生旋回变化的分析方法的吉祥泥炭剖面植硅体古年均温各旋回相对变化幅度图;
图6为本发明实施例提供的一种连续沉积物古气候代用指标他生旋回变化的分析方法的吉祥泥炭剖面植硅体古年均温各旋回时间跨度变化幅度图;
图7为本发明实施例提供的一种连续沉积物古气候代用指标他生旋回变化的分析方法的泥炭植硅体古年均温各旋回相对变化与时间跨度相对变化散点图;
图8为本发明实施例提供的一种连续沉积物古气候代用指标他生旋回变化的分析方法的吉祥泥炭剖面植硅体古年均温各旋回变率曲线图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
请一并参阅图1至图8,现对本发明提供的一种连续沉积物古气候代用指标他生旋回变化的分析方法的一种具体实施方式进行说明。该种连续沉积物古气候代用指标他生旋回变化的分析方法的技术方案是通过古气候代用指标的曲线形态分析来实现的,具体步骤包括:
S1、确定用于分析的古气候代用指标,获取沉积剖面上不同深度采样点沉积物的古气候代用指标参数ai及其对应深度值di(cm),其中古气候代用指标可以是岩石矿物、孢粉、植硅体、有孔虫、变形虫、硅藻、常量元素、微量元素、碳同位素、氧同位素、氮同位素、磁化率、粒度、色度、腐殖化度等的任意一种;
S2、选择典型深度的沉积物进行测年,获取其绝对年龄校正值,并通过线性插值法获取剖面上不同深度采样点沉积物的年龄,估计沉积旋回的时间尺度,其中典型深度指沉积剖面上沉积物颜色、岩性变化处所对应深度;
S3、根据古气候代用指标参数ai及其对应深度值di,获取沉积物古气候代用指标在剖面上的旋回cj并确定典型旋回c0及参数代表值a0和旋回深度跨度b0
S3具体的实现步骤可以为:
S3.1、根据古气候代用指标参数ai及其对应深度值di,获取剖面上沉积物古气候代用指标参数ai随深度di的变化曲线C;
S3.2、划分沉积物古气候代用指标在剖面上的旋回cj及各旋回的曲线形态类型;
S3.3、通过分析各旋回曲线形态特征,确定剖面上沉积物古气候代用指标参数变化曲线的典型旋回c0及典型旋回c0的参数代表值a0和旋回深度跨度b0
S4、基于典型旋回c0及古气候代用指标参数a0和深度跨度b0,获取各旋回cj的古气候代用指标参数的相对变化幅度
Figure BDA0002983792710000051
和各旋回cj的深度跨度的相对变化幅度ψj
S4具体的实现步骤可以为:
S4.1、获取各旋回cj的古气候代用指标参数ai的相对变化幅度
Figure BDA0002983792710000061
并得到
Figure BDA0002983792710000062
在剖面上的变化曲线A;
S4.1具体的实现步骤可以为:
S4.1.1、基于典型旋回c0及古气候代用指标参数a0和深度跨度b0,确定古气候代用指标参数ai在各旋回cj上的变化幅度
Figure BDA0002983792710000063
其中各旋回cj的参数代表值aj可以选择旋回cj的底部、中部或顶部样品的古气候参数值;
S4.1.2、确定沉积物古气候代用指标变化幅度
Figure BDA0002983792710000064
其中
Figure BDA0002983792710000065
S4.1.3、以深度为纵坐标,各旋回古气候参数的相对变化值
Figure BDA0002983792710000066
为横坐标,绘制各旋回cj的古气候代用指标参数相对变化幅度
Figure BDA0002983792710000067
在剖面上的变化曲线A。
S4.2、获取各旋回cj的深度跨度bj的相对变化幅度ψj,并得到ψj在剖面上的变化曲线B。
S4.2具体的实现步骤可以为:
S4.2.1、基于典型旋回c0及古气候代用指标参数a0和深度跨度b0,确定各旋回深度跨度bj,其中bj为旋回底部深度与顶部深度的差值,即bj=dj底-dj顶
S4.2.2、基于各旋回深度跨度bj,计算剖面上沉积物古气候代用指标各旋回cj的深度跨度bj相对典型旋回c0的深度跨度的变化ψj,得到ψj=(dj-d0)/d0,其中d0为旋回c0的深度跨度,dj为各旋回的深度跨度;
S4.2.3、以深度为纵坐标,各旋回cj的深度跨度bj相对典型旋回c0的变化值ψj为横坐标,绘制各旋回深度跨度相对变化值在剖面上的变化曲线B。
S5、获取各旋回cj的古气候代用指标的古年均温变率
Figure BDA0002983792710000068
Figure BDA0002983792710000069
/a0为横坐标,以ψj=(dj-d0)/d0为纵坐标,绘制各旋回c0的古气候代用指标变率幅度图,即
Figure BDA00029837927100000610
图,其中aj为各旋回cj的参数代表值,d0为旋回c0的深度跨度,dj为各旋回的深度跨度;
S6、以深度为纵坐标,各旋回的古气候代用指标的古年均温变率kj为横坐标,绘制kj在剖面上的变化曲线P,并得到各旋回cj的古气候代用指标变化散点图和各旋回cj的古气候代用指标变率曲线图;
S7、分析各旋回cj的古气候代用指标散点图、变率曲线图上各旋回的相似性、相异性,获取剖面沉积形成过程中他生旋回系统的变化过程,得到剖面上蕴含的古气候演变过程。
本发明提供一种分析沉积物中的古气候代用指标在气候等外界因素干扰下形成的他生旋回信息的方法,有效描述他生旋回信息特征,进而分析沉积剖面形成过程中经历的古气候背景,为深入分析古气候区域对比中存在的差异,促进全球变化研究提供基础技术支撑。
古气候变化涵盖年、百年、千年、万年、十万年、百万年、千万年、亿年等周期性旋回尺度。不同时间尺度的变化相互交织,进一步使古气候变化的方式和变率表现出很大的复杂性。第四纪以来地球各尺度的气候变化主要表现为逐渐变冷、快速回暖的不对称性。此处称具有这种变化特征的古气候参数曲线为钟形曲线。分析这种钟形曲线的他生旋回特征,得到古气候代用指标在沉积地层中的他生旋回的变化过程,可以在一定程度上得到相对有效的区域古气候信息,保证某一地区古气候重建结果解释的相对准确,促进古气候全球对比研究的科学进行。
与现有技术相比,该种连续沉积物古气候代用指标他生旋回变化的分析方法可以清楚地分析剖面上古气候相对变化的过程中的幅度和程度;可以相对简单、明了地将多个他生旋回分离开,对分析古气候等外界因素的影响因素和控制机制有很重要的作用,有助于快速进行古气候重建结果的解释,促进全球对比工作的有效开展。本发明解决现在技术所存在的因古气候代用指标的综合信息中他生旋回和自生旋回的叠加,强烈干扰古气候信息的有效性,影响古气候重建结果分析的准确性和全球对比研究进行的问题,作用效果显著,适于广泛推广。
示例性的:对长白山地区吉祥泥炭连续沉积进行分析。AMS14C测年结果显示该剖面形成于3500yrB.P.以来,属于晚全新世地层。使用该种连续沉积物古气候代用指标他生旋回变化的分析方法的具体分析步骤如下:
1、获取泥炭沉积剖面上不同深度采样点泥炭植硅体反演的古年均温ai及其对应深度值di(cm);
2、选择典型深度的泥炭沉积物进行AMS14C测年,获取其绝对年龄分别为2400±40yrB.P.,1790±40yrB.P.,及1230±35yrB.P.,通过线性插值法计算得到剖面底部泥炭的年龄约3500yrB.P.,据此确定该晚全新世剖面古气候存在百年尺度周期旋回;
3、获取剖面上泥炭植硅体古年均温ai在剖面的变化曲线C;
4、划分泥炭植硅体古年均温ai在剖面上的周期性旋回cj,即c0、c1、c2、c3、c4、c5、c6,且曲线形态类型以钟形为主;
5、分析各旋回曲线形态特征,确定曲线C上底部第一个旋回作为典型旋回c0
6、确定泥炭植硅体古温度参数a0在典型旋回c0上的变化幅度
Figure BDA0002983792710000081
其中a0选择旋回c0的底部348cm处的温度,即a0=a348=4.1℃;
7、确定沉积物古气候代用指标变化幅度
Figure BDA0002983792710000082
根据
Figure BDA0002983792710000083
算得
Figure BDA0002983792710000084
Figure BDA0002983792710000085
Figure BDA0002983792710000086
8、以深度为纵坐标,古气候参数相对变化值
Figure BDA0002983792710000087
为横坐标,绘制各旋回cj的古气候代用指标参数相对变化幅度在剖面上的变化曲线A,请参阅图5;
9、确定泥炭植硅体古温度参数a0在典型旋回c0上的深度跨度b0,其中b0为旋回底部深度与顶部深度的差值,即b0=d0底-d0顶=d348-d278=348-278=70cm;
10、确定泥炭植硅体古温度参数变化曲线各旋回的深度跨度的变化幅度ψj,其中ψj=(bj-b0)/b0,根据公式算得ψ1=0,ψ2=-0.79,ψ3=-0.42,ψ4=-0.68,ψ5=-0.84,ψ6=-0.63;
11、以深度为纵坐标,各旋回cj的深度跨度bj相对变化幅度值ψj为横坐标,绘制各旋回深度跨度相对变化值在剖面上的变化曲线B,请参阅图6;
12、计算各旋回cj的植硅体古年均温变率kj
Figure BDA0002983792710000088
k0=0,k1=0;k2=10.8,k3=17.3,k4=9.4,k5=1.9,k6=1.5;
13、以
Figure BDA0002983792710000089
为横坐标,以ψj=(dj-d0)/d0为纵坐标,绘制各旋回c0的古气候代用指标变率幅度图,
Figure BDA00029837927100000810
图,请参阅图7;
14、以深度为纵坐标,各旋回的植硅体古年均温变率kj为横坐标,绘制kj在剖面上的变化曲线P,请参阅图8;
15、分析各旋回cj的古气候代用指标散点图、变率图上各旋回的相似性、相异性,分析剖面沉积形成过程中他生旋回系统的变化过程,分析剖面上蕴含的古气候演变过程。
最终观察得到,连续沉积物泥炭中的植硅体古年均温气候变化曲线涵盖了7个逐渐变冷,快速变暖的周期性旋回c0、c1、c2、c3、c4、c5、c6,曲线形态类型以钟形为主。c0旋回为底部含泥炭的沙砾层,属于剖面早期形成阶段,选择它作为整个剖面的对比依据。
由上述数据可见,剖面上古年均温第二个气候旋回c1开始时的年均温度较c0旋回偏高,反映了相对较高温度的第二个旋回的开始,之后的第三个,第四个和第五个旋回相比c0而言,初始时的古年均温较低,但落幅不大;第五个旋回c6和第六个旋回c7初始时的古年均温比前五个旋回都低,显示为整个剖面上最冷的2个气候周期。旋回深度跨度变化幅度以c0、c1相对稳定且大,c2迅速减小,c3、c4、c5、c6则都明显小于c0和c1;这表明古气候系统逐渐变得不稳定或稳定时长相对较小。古年均温变化的
Figure BDA0002983792710000091
图所指示的是不同旋回受沉积系统外界因素干扰差异的表现,由于3500年内本区的地貌地质背景变化不大,因此可将其控制因素理解为气候控制因素。图7和图8显示,早期c0和c1两个旋回的初始年均温和深度变化幅度相似,归为一组,其气候背景均为相对稳定的暖湿特征;中期c2、c3和c4旋回的初始年均温和深度变化幅度相似,归为一组,其气候背景较早期冷凉,且不稳定;剖面上部c5和c6则为一组,其气候背景继续变得较凉而不稳定。所以,
Figure BDA0002983792710000092
图可以显示相似性质的他生旋回,对分析古气候变化的影响因素和机制有重要作用。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其它实施例的不同之处,各个实施例之间相同或相似部分互相参见即可。本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (5)

1.一种连续沉积物古气候代用指标他生旋回变化的分析方法,其特征在于,包括:
步骤1、确定用于分析的古气候代用指标,获取沉积剖面上不同深度采样点沉积物的古气候代用指标参数ai及其对应深度值di
步骤2、选择典型深度的沉积物进行测年,获取其绝对年龄校正值,并通过线性插值法获取剖面上不同深度采样点沉积物的年龄,估计沉积旋回的时间尺度;
步骤3、根据古气候代用指标参数ai及其对应深度值di,获取沉积物古气候代用指标在剖面上的旋回cj并确定典型旋回c0及参数代表值a0和旋回深度跨度b0
步骤4、基于典型旋回c0及古气候代用指标参数a0和深度跨度b0,获取各旋回cj的古气候代用指标参数的相对变化幅度
Figure FDA0002983792700000011
和各旋回cj的深度跨度的相对变化幅度ψj
步骤5、获取各旋回cj的古气候代用指标的古年均温变率
Figure FDA0002983792700000012
Figure FDA0002983792700000013
Figure FDA0002983792700000014
为横坐标,以ψj=(dj-d0)/d0为纵坐标,绘制各旋回c0的古气候代用指标变率幅度图,即
Figure FDA0002983792700000015
图,其中aj为各旋回cj的参数代表值,d0为旋回c0的深度跨度,dj为各旋回的深度跨度;
步骤6、以深度为纵坐标,各旋回的古气候代用指标的古年均温变率kj为横坐标,绘制kj在剖面上的变化曲线P,并得到各旋回cj的古气候代用指标变化散点图和各旋回cj的古气候代用指标变率曲线图;
步骤7、分析各旋回cj的古气候代用指标散点图、变率曲线图上各旋回的相似性、相异性,获取剖面沉积形成过程中他生旋回系统的变化过程,得到剖面上蕴含的古气候演变过程。
2.根据权利要求1所述的一种连续沉积物古气候代用指标他生旋回变化的分析方法,其特征在于,所述步骤3的具体步骤包括:
步骤3.1、根据古气候代用指标参数ai及其对应深度值di,获取剖面上沉积物古气候代用指标参数ai随深度di的变化曲线C;
步骤3.2、划分沉积物古气候代用指标在剖面上的旋回cj及各旋回的曲线形态类型;
步骤3.3、通过分析各旋回曲线形态特征,确定剖面上沉积物古气候代用指标参数变化曲线的典型旋回c0及典型旋回c0的参数代表值a0和旋回深度跨度b0
3.根据权利要求2所述的一种连续沉积物古气候代用指标他生旋回变化的分析方法,其特征在于,所述步骤4的具体步骤包括:
步骤4.1、获取各旋回cj的古气候代用指标参数ai的相对变化幅度
Figure FDA0002983792700000021
并得到
Figure FDA0002983792700000022
在剖面上的变化曲线A;
步骤4.2、获取各旋回cj的深度跨度bj的相对变化幅度ψj,并得到ψj在剖面上的变化曲线B。
4.根据权利要求3所述的一种连续沉积物古气候代用指标他生旋回变化的分析方法,其特征在于,所述步骤4.1的具体步骤包括:
步骤4.1.1、基于典型旋回c0及古气候代用指标参数a0和深度跨度b0,确定古气候代用指标参数ai在各旋回cj上的变化幅度
Figure FDA0002983792700000023
其中各旋回cj的参数代表值aj选择旋回cj的底部、中部或顶部样品的古气候参数值;
步骤4.1.2、确定沉积物古气候代用指标变化幅度
Figure FDA0002983792700000024
其中
Figure FDA0002983792700000025
步骤4.1.3、以深度为纵坐标,各旋回古气候参数的相对变化值
Figure FDA0002983792700000026
为横坐标,绘制各旋回cj的古气候代用指标参数相对变化幅度
Figure FDA0002983792700000027
在剖面上的变化曲线A。
5.根据权利要求4所述的一种连续沉积物古气候代用指标他生旋回变化的分析方法,其特征在于,所述步骤4.2的具体步骤包括:
步骤4.2.1、基于典型旋回c0及古气候代用指标参数a0和深度跨度b0,确定各旋回深度跨度bj,其中bj为旋回底部深度与顶部深度的差值,即bj=dj底-dj顶
步骤4.2.2、基于各旋回深度跨度bj,计算剖面上沉积物古气候代用指标各旋回cj的深度跨度bj相对典型旋回c0的深度跨度的变化值ψj,得到ψj=(dj-d0)/d0,其中d0为旋回c0的深度跨度,dj为各旋回的深度跨度;
步骤4.2.3、以深度为纵坐标,各旋回cj的深度跨度bj相对典型旋回c0的变化值ψj为横坐标,绘制各旋回深度跨度相对变化值在剖面上的变化曲线B。
CN202110294450.6A 2021-03-19 2021-03-19 一种连续沉积物古气候代用指标他生旋回变化的分析方法 Active CN113065236B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110294450.6A CN113065236B (zh) 2021-03-19 2021-03-19 一种连续沉积物古气候代用指标他生旋回变化的分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110294450.6A CN113065236B (zh) 2021-03-19 2021-03-19 一种连续沉积物古气候代用指标他生旋回变化的分析方法

Publications (2)

Publication Number Publication Date
CN113065236A true CN113065236A (zh) 2021-07-02
CN113065236B CN113065236B (zh) 2022-03-11

Family

ID=76562383

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110294450.6A Active CN113065236B (zh) 2021-03-19 2021-03-19 一种连续沉积物古气候代用指标他生旋回变化的分析方法

Country Status (1)

Country Link
CN (1) CN113065236B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2048026A1 (en) * 1991-07-26 1993-01-27 Yavor Yossifov Shopov Apparatus for reading of luminescence records in minerals
CN103345001A (zh) * 2013-06-25 2013-10-09 中国地质大学(北京) 一种测定古湖泊水深的方法
CN103377196A (zh) * 2012-04-13 2013-10-30 中国石油化工股份有限公司 一种基于地理信息系统的岩相古地理资料处理方法
CN104749394A (zh) * 2015-04-17 2015-07-01 中国地质大学(北京) 基于沿岸沙坝厚度定量方法的古风力测量方法及装置
CN104765966A (zh) * 2015-04-17 2015-07-08 中国地质大学(北京) 一种古湖泊中沿岸沙坝原始厚度的测量方法及装置
CN104898181A (zh) * 2015-06-02 2015-09-09 中国地质大学(北京) 古时期的最大理论波高和累计频率波高确定方法
CN106019402A (zh) * 2016-05-12 2016-10-12 山东科技大学 综合判别与验证巨厚煤层内部古气候旋回的方法
CN106204591A (zh) * 2016-07-11 2016-12-07 青岛海洋地质研究所 一种海域沉积物岩心中古环境信息载体的取样方法
CN107153215A (zh) * 2017-04-19 2017-09-12 中国石油大学(华东) 一种基于干旱气候约束下的沉积模式构建方法
CN109685402A (zh) * 2019-02-21 2019-04-26 中国石油化工股份有限公司 一种湖泊沉积环境的综合量化恢复方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2048026A1 (en) * 1991-07-26 1993-01-27 Yavor Yossifov Shopov Apparatus for reading of luminescence records in minerals
CN103377196A (zh) * 2012-04-13 2013-10-30 中国石油化工股份有限公司 一种基于地理信息系统的岩相古地理资料处理方法
CN103345001A (zh) * 2013-06-25 2013-10-09 中国地质大学(北京) 一种测定古湖泊水深的方法
CN104749394A (zh) * 2015-04-17 2015-07-01 中国地质大学(北京) 基于沿岸沙坝厚度定量方法的古风力测量方法及装置
CN104765966A (zh) * 2015-04-17 2015-07-08 中国地质大学(北京) 一种古湖泊中沿岸沙坝原始厚度的测量方法及装置
CN104898181A (zh) * 2015-06-02 2015-09-09 中国地质大学(北京) 古时期的最大理论波高和累计频率波高确定方法
CN106019402A (zh) * 2016-05-12 2016-10-12 山东科技大学 综合判别与验证巨厚煤层内部古气候旋回的方法
CN106204591A (zh) * 2016-07-11 2016-12-07 青岛海洋地质研究所 一种海域沉积物岩心中古环境信息载体的取样方法
CN107153215A (zh) * 2017-04-19 2017-09-12 中国石油大学(华东) 一种基于干旱气候约束下的沉积模式构建方法
CN109685402A (zh) * 2019-02-21 2019-04-26 中国石油化工股份有限公司 一种湖泊沉积环境的综合量化恢复方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
姜文超 等: "沉积学发展及油气勘探领域热点问题概述", 《地质与资源》 *
方石 等: "松辽盆地合隆地区泉头组河流沉积体系判别及其形成条件分析", 《吉林大学学报(地球科学版)》 *
王治祥: "中新世以来轨道尺度的古气候变化在青藏高原东北缘湖盆记录中的沉积响应", 《中国优秀博硕士学位论文全文数据库(博士)基础科学辑》 *

Also Published As

Publication number Publication date
CN113065236B (zh) 2022-03-11

Similar Documents

Publication Publication Date Title
Littke et al. Reflectance of dispersed vitrinite in Palaeozoic rocks with and without cleavage: Implications for burial and thermal history modeling in the Devonian of Rursee area, northern Rhenish Massif, Germany
Huang et al. Astronomical forcing of Middle Permian terrestrial climate recorded in a large paleolake in northwestern China
Wright et al. Sea-level fluctuations driven by changes in global ocean basin volume following supercontinent break-up
Fitzpatrick et al. Physical properties of the WAIS Divide ice core
Müller et al. A tectonic-rules-based mantle reference frame since 1 billion years ago–implications for supercontinent cycles and plate–mantle system evolution
CN108846540A (zh) 致密砂岩气田的采收率标定方法及装置
CN109212627B (zh) 一种基于细分层系的薄储层描述及增储规模确定方法
CA2896072A1 (fr) Procede d'exploitation d'un gisement petrolier a partir d'une technique de positionnement des puits a forer
CN111472765A (zh) 目标井的地层划分方法和装置
CN113420936A (zh) 一种基于动态时间规划优化算法的地层对比方法
CN113065236B (zh) 一种连续沉积物古气候代用指标他生旋回变化的分析方法
Maurer et al. Statistical time-series analysis and sedimentological tuning of bedding rhythms in a Triassic basinal succession (Southern Alps, Italy)
Carlucci et al. Trilobite biofacies along an Ordovician (Sandbian) carbonate buildup to basin gradient, southwestern Virginia
CN111983671B (zh) 基于微古地貌恢复的浅水湖盆储层预测方法和装置
CN110389381B (zh) 一种基于地震属性的砂体储层预测方法及装置
Wang et al. Late Triassic sedimentary records reveal the hydrological response to climate forcing and the history of the chaotic Solar System
CN110879920B (zh) 一种获取沉积盆地坳陷区古地表热流的方法、装置及系统
Ping et al. The application of RASC/CASC methods to quantitative biostratigraphic correlation of Neogene in northern South China Sea
Yu et al. A high-resolution timescale for the Miocene Shanwang diatomaceous shale lagerstätte (China): development of Wavelet Scale Series Analysis for cyclostratigraphy
Tang et al. Astronomical forcing of favorable sections of lacustrine tight reservoirs in the lower Shangganchaigou formation of the Gasi Area, Western Qaidam Basin
CN113138414A (zh) 深层砂岩储层孔隙度的预测方法及装置
Shen et al. Quantitative palaeogeographical reconstruction of the North China Block during the Carboniferous and Permian transition: Implications for coal accumulation and source rock development
CN110133719B (zh) 基于90°相移数据体、种子点追踪的砂砾岩期次划分方法
Batista-Rodríguez et al. Geometry of ophiolites in Eastern Cuba from 3D inversion of aeromagnetic data, constrained by surface geology
CN114740537A (zh) 一种页岩层系等时对比方法及电子设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant