一種微古地貌恢復(fù)方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及油氣及煤層氣地震勘探與開發(fā)領(lǐng)域，特別是涉及一種微古地貌恢復(fù)方 法。
【背景技術(shù)】
[0002] 古地貌是控制含油氣盆地沉積相發(fā)育與分布的主要因素，在一定程度上控制著后 期油藏的儲(chǔ)蓋組合，精確的古地貌恢復(fù)對(duì)油氣勘探至關(guān)重要。古地貌恢復(fù)即利用地震測井 資料恢復(fù)出某個(gè)年代時(shí)的地形形態(tài)。古地貌的恢復(fù)有助于認(rèn)識(shí)古地理、古生物的分布以及 構(gòu)造演化特征等，也有助于揭示物源與沉積體系的配置關(guān)系，揭示有效砂體和優(yōu)質(zhì)儲(chǔ)層的 控制因素及分布規(guī)律，從而有效的指導(dǎo)油氣勘探。
[0003] 傳統(tǒng)的古地貌恢復(fù)方法有殘留厚度和補(bǔ)償厚度印模法、回剝和填平補(bǔ)齊法。近年 發(fā)展了沉積學(xué)分析方法，并將層序地層學(xué)理論應(yīng)用于古地貌恢復(fù)，形成了層拉平古地貌恢 復(fù)法，結(jié)合陸相地層的特征，發(fā)展了高分辨率層序地層學(xué)法。但是，目前對(duì)古地貌恢復(fù)的研 究大都停留在定性階段，無法進(jìn)行定量研究，沉積記錄資料的多少?zèng)Q定了恢復(fù)精度的高低。 微古地貌的恢復(fù)則需要有更多的沉積記錄資料以及更精確的恢復(fù)方法。

【發(fā)明內(nèi)容】

[0004] 本發(fā)明的目的是要提供一種恢復(fù)精度較高、可定量研究的微古地貌的恢復(fù)方法。
[0005] 本申請(qǐng)的發(fā)明人發(fā)現(xiàn)，對(duì)于當(dāng)前的一目標(biāo)地形，其上的沉積體分布與該目標(biāo)地形 對(duì)應(yīng)的古地貌之間具有對(duì)應(yīng)關(guān)系。換句話說，盡管在一個(gè)古地貌上形成沉積體之后會(huì)因地 質(zhì)構(gòu)造運(yùn)動(dòng)而使得地形發(fā)生變化直至演化成當(dāng)前的目標(biāo)地形，但是其上的沉積體依然會(huì)保 留有其在古地貌時(shí)的沉積體基本相似的分布或者說形態(tài)。本申請(qǐng)的發(fā)明人進(jìn)一步認(rèn)識(shí)到， 如果將微古地貌對(duì)應(yīng)于顱骨，將沉積體對(duì)應(yīng)于軟組織，那么微古地貌與沉積體之間的聯(lián)系 可以類比于顱骨與軟組織之間的聯(lián)系，因此可以將類似于顱面復(fù)原的技術(shù)應(yīng)用于微古地貌 恢復(fù)。
[0006] 基于發(fā)明人的上述發(fā)現(xiàn)和認(rèn)識(shí)，本發(fā)明提供了一種微古地貌恢復(fù)方法，用于基于 目標(biāo)地形恢復(fù)出與所述目標(biāo)地形對(duì)應(yīng)的微古地貌，包括：
[0007] 測量所述目標(biāo)地形上的當(dāng)前沉積體分布；
[0008] 在一個(gè)或多個(gè)已知沉積體分布中確定與所述當(dāng)前沉積體分布相似的一個(gè)已知沉 積體分布作為參考沉積體分布；其中，每一所述已知沉積體分布對(duì)應(yīng)于一已知古地貌，將所 述參考沉積體分布對(duì)應(yīng)的已知古地貌作為參考古地貌；
[0009] 確定將所述參考沉積體分布變形至所述當(dāng)前沉積體分布的第一變形關(guān)系；和 [0010] 利用所述第一變形關(guān)系將所述參考古地貌進(jìn)行變形，并將變形得到的結(jié)果作為與 所述目標(biāo)地形對(duì)應(yīng)的所述微古地貌。
[0011] 可選地，該方法還包括：
[0012] 計(jì)算所述當(dāng)前沉積體分布與所述多個(gè)已知沉積體分布中任一已知沉積體分布的 相似度；
[0013] 其中，將具有最高相似度的已知沉積體分布作為所述參考沉積體分布。
[0014] 可選地，所述相似度由所述當(dāng)前沉積體分布與所述已知沉積體分布之間的第二變 形關(guān)系來表征；
[0015] 其中，所述第二變形關(guān)系所對(duì)應(yīng)的變形量越小則所述相似度越高。
[0016] 可選地，所述相似度是基于所述當(dāng)前沉積體分布與所述已知沉積體分布的多個(gè)對(duì) 應(yīng)的特征點(diǎn)的數(shù)據(jù)計(jì)算得到的。
[0017] 可選地，該方法還包括：
[0018] 采用特征提取算法確定所述當(dāng)前沉積體分布和/或所述已知沉積體分布的特征 點(diǎn)；可選地，所述特征提取算法包括在顱面恢復(fù)中使用的自動(dòng)或半自動(dòng)特征點(diǎn)標(biāo)定方法。
[0019] 可選地，所述第二變形關(guān)系包括所述當(dāng)前沉積體分布與所述已知沉積體分布中的 一個(gè)變形到另一個(gè)的平移矩陣和/或旋轉(zhuǎn)矩陣。
[0020] 可選地，將具有最小的模的所述平移矩陣和/或所述旋轉(zhuǎn)矩陣所對(duì)應(yīng)的所述已知 沉積體分布作為所述參考沉積體分布。
[0021] 可選地，將所述當(dāng)前沉積體分布與所述參考沉積體分布之間的所述第二變形關(guān)系 作為所述第一變形關(guān)系。
[0022] 可選地，所述相似度是采用配準(zhǔn)算法計(jì)算得到的；可選地，所述配準(zhǔn)算法包括三維 點(diǎn)云配準(zhǔn)算法；可選地，所述三維點(diǎn)云配準(zhǔn)算法包括主成分分析方法。
[0023] 可選地，該方法還包括：
[0024] 建立一數(shù)據(jù)庫，所述數(shù)據(jù)庫中存儲(chǔ)有多種已知古地貌以及各自對(duì)應(yīng)的已知沉積體 分布的數(shù)據(jù)。
[0025] 本發(fā)明從顱骨復(fù)原的角度，將類似于顱面復(fù)原的技術(shù)應(yīng)用于微古地貌的恢復(fù)，通 過計(jì)算參考沉積體分布與當(dāng)前沉積體分布的變形關(guān)系并將該變形關(guān)系應(yīng)用于參考古地貌 使其進(jìn)行變形，從而得到與目標(biāo)地形對(duì)應(yīng)的微古地貌。本發(fā)明提供的微古地貌的恢復(fù)方法， 使古地貌恢復(fù)精度提高至一個(gè)新的高度，實(shí)現(xiàn)了其恢復(fù)的定量化。并且，本發(fā)明提供的微古 地貌的恢復(fù)方法可以對(duì)多處不同的地形進(jìn)行恢復(fù)，具有很強(qiáng)的適用性。
[0026] 根據(jù)下文結(jié)合附圖對(duì)本發(fā)明具體實(shí)施例的詳細(xì)描述，本領(lǐng)域技術(shù)人員將會(huì)更加明 了本發(fā)明的上述以及其他目的、優(yōu)點(diǎn)和特征。
【附圖說明】
[0027] 后文將參照附圖以示例性而非限制性的方式詳細(xì)描述本發(fā)明的一些具體實(shí)施例。 附圖中相同的附圖標(biāo)記標(biāo)示了相同或類似的部件或部分。本領(lǐng)域技術(shù)人員應(yīng)該理解，這些 附圖未必是按比例繪制的。附圖中：
[0028]圖1是根據(jù)本發(fā)明微古地貌恢復(fù)方法的一個(gè)實(shí)施例的流程圖；
[0029] 圖2是將參考沉積體分布Ps通過第一變形關(guān)系T/變形到當(dāng)前沉積體分布Ρε的結(jié) 構(gòu)示意圖；
[0030] 圖3是將參考古地貌Gs通過第一變形關(guān)系Τ/變形就得到了目標(biāo)地形所對(duì)應(yīng)的微 古地貌的結(jié)構(gòu)示意圖。
【具體實(shí)施方式】
[0031] 圖1是根據(jù)本發(fā)明的微古地貌恢復(fù)方法的一個(gè)實(shí)施例的流程圖，其用于基于目標(biāo) 地形恢復(fù)出與所述目標(biāo)地形對(duì)應(yīng)的微古地貌。該方法可以包括步驟S101，測量目標(biāo)地形上 的當(dāng)前沉積體分布Ρε。在步驟S101中，可利用地震、測井及鉆井資料來獲取所需的當(dāng)前沉 積體分布ρε，這是本領(lǐng)域技術(shù)人員所熟知的。
[0032] 在步驟S103中，可以在一個(gè)或多個(gè)已知沉積體分WPe(E= 1、2、…、η，η為大于 等于1的正整數(shù)）中確定與當(dāng)前沉積體分布Ρε相似的一個(gè)已知沉積體分布作為參考沉積 體分布1\，其中Ps為Ρ:至Ρ"中的一個(gè)。每一已知沉積體分布ΡΕ可以對(duì)應(yīng)于一個(gè)已知古地 貌GE。對(duì)應(yīng)地，可以將該參考沉積體分布Ps對(duì)應(yīng)的已知古地貌作為參考古地貌Gs。這些已 知的沉積體分布PE和古地貌GE可以從現(xiàn)有的經(jīng)驗(yàn)知識(shí)中收集獲得。例如，可以將通過傳統(tǒng) 的各種古地貌恢復(fù)方法恢復(fù)的古地貌及其對(duì)應(yīng)的沉積體分布作為這里的已知的沉積體分 布PjP古地貌Ge。
[0033] 在一個(gè)實(shí)施例中，可以建立一數(shù)據(jù)庫DB，并將上述的各個(gè)已知古地貌心以及各自 對(duì)應(yīng)的已知沉積體分布匕的數(shù)據(jù)存儲(chǔ)在該數(shù)據(jù)庫DB中。下面的表1示意性地示出了該數(shù) 據(jù)庫DB。
[0034] 表 1
[0035]
[0036] 在一個(gè)實(shí)施例中，^以在步驟SIO^中計(jì)算當(dāng)前沉'體分布Ρε與任一已知沉積體 分布？1；的相似度，并將具有最高相似度的已知沉積體分布ΡΕ作為參考沉積體分布Ps。在一 個(gè)實(shí)施例中，該相似度可以由當(dāng)前沉積體分布匕與各個(gè)已知沉積體分布ρ￡之間的第二變形 關(guān)系Τ2Ε來表征。第二變形關(guān)系Τ2Ε所對(duì)應(yīng)的變形量越小則表明相似度越高。
[0037] 為了方便且準(zhǔn)確地計(jì)算該相似度或第二變形關(guān)系Τ2Ε，可以采用特征提