4. 单井埋藏史曲线判读实战
埋藏史曲线是把一口井从沉积到现今的“深度—时间”轨迹画在一张图上，横轴为地质时间（Ma，向左变老），纵轴为埋深（m，向下增大）。判读的核心不是读点，而是读斜率、拐点、间断和层序关系。本节给出一口虚拟井（示例井 X1）的设定数据，按“分层回剥—剥蚀量估算—沉降阶段划分—事件标注”的顺序完整走一遍，所有数值均为教学假设，用于演示方法而非实测结果。
一、示例井的基础数据与参照约定
约定：深度以地表为 0 m、向下为正；地层年龄取底界年龄；剥蚀量指被剥蚀掉的地层厚度，用镜质体反射率、孔隙度—深度趋势或区域地层对比三种方法交叉约束。示例井 X1 的地层与关键界面假设如下：
| 地层/界面 | 底界年龄(Ma) | 现今底界深度(m) | 岩性 |
|---|---|---|---|
| 第四系—新近系 Q+N | 23 | 600 | 砂泥互层 |
| 古近系 E | 65 | 1800 | 泥岩为主 |
| 白垩系 K（不整合面） | 100 | 2400 | 砂岩夹泥岩 |
| 侏罗系 J | 160 | 3100 | 煤系泥岩 |
| 三叠系 T（基底顶） | 250 | 3600 | 致密砂岩 |
回剥的基本关系是：某界面的古埋深 = 现今埋深 − 上覆各地层在“该界面沉积之后”所经历的厚度恢复量。对每一层用孔隙度—深度关系 φ = φ₀·e^(−cz) 做去压实校正，把现今厚度换算为地质历史中的厚度。示例取 φ₀ = 0.55、c = 0.00045 m⁻¹，泥岩层压实效应明显，砂岩层近似刚性。这一步决定了曲线在深部的形态是否被“压扁”，是判读前必须完成的预处理。
二、剥蚀量的估算与不整合识别
X1 井在 K 底界（100 Ma）处存在明显不整合：上覆古近系底界 65 Ma 直接覆盖在白垩系 100 Ma 之上，缺失 35 Ma 的地层记录。曲线在此处表现为埋深突然变浅的“回跳”——这是剥蚀而非沉积的特征。估算剥蚀量可用三种手段互相校验：
- Ro 法：不整合面之下 Ro 值出现突变台阶，用 Ro 与埋深关系反推最大古埋深，再减去现今埋深即得剥蚀量。示例中不整合面下 Ro 由 0.85% 跳至 1.15%，反推剥蚀量约 500 m。
- 孔隙度趋势外推：把不整合面以下的孔隙度—深度趋势线向上延伸至地表，趋势线与现今地表之间的“缺失厚度”即剥蚀量，示例得 450～550 m。
- 地层对比法：邻井保留有该期地层，对比厚度差得 480 m。
三者取 500 m 作为 X1 井 K 不整合面的剥蚀量估计值，并在曲线上把 100 Ma 处的埋深点向上修正 500 m 后再继续向老地层回剥。若不修正，深部地层的古埋深会被系统性低估，热演化史也随之偏低。
三、沉降阶段划分与曲线形态判读
把 X1 井曲线按斜率变化切成四段，每段的判读依据如下：
- ① 250～160 Ma（三叠纪—侏罗纪）：快速沉降。曲线陡直，100 Ma 内埋深由 0 增至约 3100 m，平均沉降速率约 34 m/Ma，对应盆地断陷或快速充填期，泥岩压实快，是烃源岩成熟的关键埋深窗口。
- ② 160～100 Ma（白垩纪）：稳定沉降。斜率明显变缓，60 Ma 内仅增加约 700 m，速率约 12 m/Ma，属坳陷期缓慢沉降，地层厚度稳定，是储层发育的有利阶段。
- ③ 100～65 Ma：抬升剥蚀。曲线出现回跳，缺失 35 Ma 记录，剥蚀约 500 m，对应区域挤压抬升事件，需重点标注不整合。
- ④ 65～0 Ma（古近纪至今）：再次沉降。曲线重新变陡，65 Ma 内埋深由约 1900 m 增至 2400 m 以上，速率约 8～10 m/Ma，属坳陷再沉降，决定现今成熟度。
判读时要注意：斜率是“视沉降速率”，必须用去压实后的厚度计算，否则深部会因压实而显得沉降变慢，造成阶段误判。此外，曲线上的每一个拐点都应回到地层表核对年龄，避免把地层界面误读为构造事件。
四、关键地质事件标注与成果表达
完成上述三步后，在曲线上标注四类事件：地层界面（T、J、K、E 底界）、不整合与剥蚀量（K 面，500 m）、沉降阶段界线（160、100、65 Ma）、成熟度窗口（Ro = 0.5% 与 1.0% 对应的古埋深区间）。示例中 Ro = 0.5% 约出现在 120 Ma、埋深 2800 m 处，Ro = 1.0% 约出现在 90 Ma、埋深 3300 m 处，据此可判断主力烃源岩在早白垩世末进入生烃窗，抬升剥蚀期生烃一度停滞，古近纪再沉降后二次生烃。
按上述流程，一口虚拟井的埋藏史曲线就能从“一条折线”转化为包含埋深恢复、剥蚀量、沉降阶段和事件时序的完整判读成果，为后续热演化史与生烃史分析提供可靠输入。