[發明專利]一種富含有機質泥頁巖孔隙度的演化過程的恢復方法

有效

• 申請號: 201710515666.4
• 申請日: 2017-06-29
• 公開（公告）號: CN109211745B
• 公開（公告）日: 2020-10-20
• 發明（設計）人: 王保華;陸建林;李浩;宋振響;左宗鑫;王苗;宋在超
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發研究院
• 主分類號: G01N15/08
• 分類號: G01N15/08
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 吳大建;張杰

1.一種富含有機質泥頁巖孔隙度的演化過程的恢復方法，包括：

在目標工區的目的層內采樣得到多個樣品，并獲取分別對應于每個樣品的的一組實測數據，以及所述目的層的應力場分布數據；

根據所有樣品的實測數據獲得有機碳恢復系數關于成熟度的第一關系表達式，以及生烴孔隙度與構造孔隙度的和關于生烴消耗的有機碳的第二關系表達式，包括：

基于如下表達式建立所述第一關系表達式：

f_恢復(Ro(x,y))＝A₁(Ro(x,y))³+A₂(Ro(x,y))²+A₃Ro(x,y)+A₄

其中，f_恢復(Ro(x,y))表示目的層內對應于地質點(x,y)處的所述有機碳恢復系數，Ro(x,y)表示目的層內對應于地質點(x,y)處的所述成熟度，A₁、A₂、A₃、A₄分別表示通過統計回歸確定的系數，

基于如下表達式建立所述第二關系表達式：

Φ_生+構(ΔC)＝B₁ΔC+B₂

其中，Φ_生+構(ΔC)表示生烴孔隙度與構造孔隙度的和，ΔC表示生烴消耗的有機碳，B₁、B₂分別表示通過統計回歸確定的系數；

根據所述第一關系表達式與所述第二關系表達式，建立所述生烴孔隙度關于所述成熟度的第一數學模型，包括：

基于如下表達式建立所述第一數學模型：

Φ_生(x,y)＝Φ_生+構((f_恢復(Ro(x,y))-1)·Toc_殘(x,y))-Φ_構(x,y)

其中，Φ_生+構(·)表示所述生烴孔隙度與構造孔隙度的和關于生烴消耗的有機碳的第二關系表達式，f_恢復(Ro(x,y))表示目的層內對應于地質點(x,y)處的所述有機碳恢復系數，Ro(x,y)表示目的層內對應于地質點(x,y)處的所述成熟度，Toc_殘(x,y)表示目的層內對應于地質點(x,y)處的所述殘余有機碳，Φ_生(x,y)表示目的層內對應于地質點(x,y)處的所述生烴孔隙度，Φ_構(x,y)表示目的層內對應于地質點(x,y)處的所述構造孔隙度；

根據所述第二關系表達式及所述應力場分布數據，建立所述構造孔隙度關于構造應力的第二數學模型，包括：

基于如下表達式建立所述第二數學模型：

其中，Φ_構(x,y)表示目的層內對應于地質點(x,y)處的所述構造孔隙度，P(x,y)表示目的層內對應于地質點(x,y)處的所述構造應力，Φ_構0表示所述未發生生烴演化時的構造孔隙度的值，P₀表示與所述未發生生烴演化時的構造孔隙度的值相對應的構造應力值，P₁與P₂分別表示對應于所述構造應力的取值范圍的最小值與最大值，Φ_構1與Φ_構2分別表示對應于所述構造孔隙度的取值范圍的最小值與最大值；

根據所述第一數學模型、所述第二數學模型以及泥頁巖壓實模型，確定所述目的層內每個地質點在其地質演化的各個歷史階段中的孔隙度的數值。

2.根據權利要求1所述的恢復方法，其特征在于，所述實測數據包括埋深、成熟度、殘余有機碳以及總孔隙度；所述應力場分布數據包括根據所述目的層的構造埋深、地層厚度、巖性及實測的應力數據，并采用薄板彎曲模型模擬計算得到的數據。