一种微古地貌恢复方法
【技术领域】
[0001] 本发明涉及油气及煤层气地震勘探与开发领域，特别是涉及一种微古地貌恢复方 法。
【背景技术】
[0002] 古地貌是控制含油气盆地沉积相发育与分布的主要因素，在一定程度上控制着后 期油藏的储盖组合，精确的古地貌恢复对油气勘探至关重要。古地貌恢复即利用地震测井 资料恢复出某个年代时的地形形态。古地貌的恢复有助于认识古地理、古生物的分布以及 构造演化特征等，也有助于揭示物源与沉积体系的配置关系，揭示有效砂体和优质储层的 控制因素及分布规律，从而有效的指导油气勘探。
[0003] 传统的古地貌恢复方法有残留厚度和补偿厚度印模法、回剥和填平补齐法。近年 发展了沉积学分析方法，并将层序地层学理论应用于古地貌恢复，形成了层拉平古地貌恢 复法，结合陆相地层的特征，发展了高分辨率层序地层学法。但是，目前对古地貌恢复的研 究大都停留在定性阶段，无法进行定量研究，沉积记录资料的多少决定了恢复精度的高低。 微古地貌的恢复则需要有更多的沉积记录资料以及更精确的恢复方法。
【发明内容】
[0004] 本发明的目的是要提供一种恢复精度较高、可定量研究的微古地貌的恢复方法。
[0005] 本申请的发明人发现，对于当前的一目标地形，其上的沉积体分布与该目标地形 对应的古地貌之间具有对应关系。换句话说，尽管在一个古地貌上形成沉积体之后会因地 质构造运动而使得地形发生变化直至演化成当前的目标地形，但是其上的沉积体依然会保 留有其在古地貌时的沉积体基本相似的分布或者说形态。本申请的发明人进一步认识到， 如果将微古地貌对应于颅骨，将沉积体对应于软组织，那么微古地貌与沉积体之间的联系 可以类比于颅骨与软组织之间的联系，因此可以将类似于颅面复原的技术应用于微古地貌 恢复。
[0006] 基于发明人的上述发现和认识，本发明提供了一种微古地貌恢复方法，用于基于 目标地形恢复出与所述目标地形对应的微古地貌，包括：
[0007] 测量所述目标地形上的当前沉积体分布；
[0008] 在一个或多个已知沉积体分布中确定与所述当前沉积体分布相似的一个已知沉 积体分布作为参考沉积体分布；其中，每一所述已知沉积体分布对应于一已知古地貌，将所 述参考沉积体分布对应的已知古地貌作为参考古地貌；
[0009] 确定将所述参考沉积体分布变形至所述当前沉积体分布的第一变形关系；和 [0010] 利用所述第一变形关系将所述参考古地貌进行变形，并将变形得到的结果作为与 所述目标地形对应的所述微古地貌。
[0011] 可选地，该方法还包括：
[0012] 计算所述当前沉积体分布与所述多个已知沉积体分布中任一已知沉积体分布的 相似度；
[0013] 其中，将具有最高相似度的已知沉积体分布作为所述参考沉积体分布。
[0014] 可选地，所述相似度由所述当前沉积体分布与所述已知沉积体分布之间的第二变 形关系来表征；
[0015] 其中，所述第二变形关系所对应的变形量越小则所述相似度越高。
[0016] 可选地，所述相似度是基于所述当前沉积体分布与所述已知沉积体分布的多个对 应的特征点的数据计算得到的。
[0017] 可选地，该方法还包括：
[0018] 采用特征提取算法确定所述当前沉积体分布和/或所述已知沉积体分布的特征 点；可选地，所述特征提取算法包括在颅面恢复中使用的自动或半自动特征点标定方法。
[0019] 可选地，所述第二变形关系包括所述当前沉积体分布与所述已知沉积体分布中的 一个变形到另一个的平移矩阵和/或旋转矩阵。
[0020] 可选地，将具有最小的模的所述平移矩阵和/或所述旋转矩阵所对应的所述已知 沉积体分布作为所述参考沉积体分布。
[0021] 可选地，将所述当前沉积体分布与所述参考沉积体分布之间的所述第二变形关系 作为所述第一变形关系。
[0022] 可选地，所述相似度是采用配准算法计算得到的；可选地，所述配准算法包括三维 点云配准算法；可选地，所述三维点云配准算法包括主成分分析方法。
[0023] 可选地，该方法还包括：
[0024] 建立一数据库，所述数据库中存储有多种已知古地貌以及各自对应的已知沉积体 分布的数据。
[0025] 本发明从颅骨复原的角度，将类似于颅面复原的技术应用于微古地貌的恢复，通 过计算参考沉积体分布与当前沉积体分布的变形关系并将该变形关系应用于参考古地貌 使其进行变形，从而得到与目标地形对应的微古地貌。本发明提供的微古地貌的恢复方法， 使古地貌恢复精度提高至一个新的高度，实现了其恢复的定量化。并且，本发明提供的微古 地貌的恢复方法可以对多处不同的地形进行恢复，具有很强的适用性。
[0026] 根据下文结合附图对本发明具体实施例的详细描述，本领域技术人员将会更加明 了本发明的上述以及其他目的、优点和特征。
【附图说明】
[0027] 后文将参照附图以示例性而非限制性的方式详细描述本发明的一些具体实施例。 附图中相同的附图标记标示了相同或类似的部件或部分。本领域技术人员应该理解，这些 附图未必是按比例绘制的。附图中：
[0028]图1是根据本发明微古地貌恢复方法的一个实施例的流程图；
[0029] 图2是将参考沉积体分布Ps通过第一变形关系T/变形到当前沉积体分布Ρε的结 构示意图；
[0030] 图3是将参考古地貌Gs通过第一变形关系Τ/变形就得到了目标地形所对应的微 古地貌的结构示意图。
【具体实施方式】
[0031] 图1是根据本发明的微古地貌恢复方法的一个实施例的流程图，其用于基于目标 地形恢复出与所述目标地形对应的微古地貌。该方法可以包括步骤S101，测量目标地形上 的当前沉积体分布Ρε。在步骤S101中，可利用地震、测井及钻井资料来获取所需的当前沉 积体分布ρε，这是本领域技术人员所熟知的。
[0032] 在步骤S103中，可以在一个或多个已知沉积体分WPe(E= 1、2、…、η，η为大于 等于1的正整数）中确定与当前沉积体分布Ρε相似的一个已知沉积体分布作为参考沉积 体分布1\，其中Ps为Ρ:至Ρ"中的一个。每一已知沉积体分布ΡΕ可以对应于一个已知古地 貌GE。对应地，可以将该参考沉积体分布Ps对应的已知古地貌作为参考古地貌Gs。这些已 知的沉积体分布PE和古地貌GE可以从现有的经验知识中收集获得。例如，可以将通过传统 的各种古地貌恢复方法恢复的古地貌及其对应的沉积体分布作为这里的已知的沉积体分 布PjP古地貌Ge。
[0033] 在一个实施例中，可以建立一数据库DB，并将上述的各个已知古地貌心以及各自 对应的已知沉积体分布匕的数据存储在该数据库DB中。下面的表1示意性地示出了该数 据库DB。
[0034] 表 1
[0035]
[0036] 在一个实施例中，^以在步骤SIO^中计算当前沉'体分布Ρε与任一已知沉积体 分布？1；的相似度，并将具有最高相似度的已知沉积体分布ΡΕ作为参考沉积体分布Ps。在一 个实施例中，该相似度可以由当前沉积体分布匕与各个已知沉积体分布ρ￡之间的第二变形 关系Τ2Ε来表征。第二变形关系Τ2Ε所对应的变形量越小则表明相似度越高。
[0037] 为了方便且准确地计算该相似度或第二变形关系Τ2Ε，可以采用特征提