[发明专利]一种基于超压抑制有机质成熟度的镜质体反射率预测模型

权利要求书

1.一种基于超压抑制有机质成熟度的镜质体反射率预测模型，其特征在于按照以下步骤进行：

(1)构建预测镜质体反射率数学模型；

(2)求取模型中引入的超压对频率因子和活化能的影响因子；

(3)验证改进后的镜质体反射率数学模型的有效性；

所述步骤(1)的方法如下：

考虑超压对镜质体热演化过程中频率因子或活化能的影响，基于EASY％R_o模型引入了频率因子(λ₁)和活化能(λ₂)两个参量，对镜质体反射率预测模型进行了改进，得到改进后的镜质体反射率预测模型A-E_a-R_o：

d_ω/d_t＝-kω (1)

k＝A·exp(-E_a/RT) (2)

F＝1-ω/ω₀ (3)

％R_o＝exp(-1.6+0.37F) (4)

d_ω′/d_t＝-k′ω′ (5)

k′＝A′·exp(-E′_a/RT) (6)

F′＝1-ω′/ω₀ (7)

R′_o＝exp(-1.6+0.37F′) (8)

A′＝A·exp(-Δp/λ₁) (9)

E′_a＝E_a+λ₂·Δp (10)

联合式(1)-(10)，整理得出本次新模型公式

式中，ω和ω′为EASY％R_o、A-E_a-R_o模型中未反应物含量；ω₀为总反应物的初始浓度；k和k′为EASY％R_o、A-E_a-R_o模型中反应速率；A和A′为EASY％R_o、A-E_a-R_o模型中频率因子A＝1.0×10¹³/s；E_a和E_a'为EASY％R_o、A-E_a-R_o模型反应活化能；R为通用气体常数；T为绝对温度；λ₁和λ₂为频率因子系数和活化能系数；ΔP为超压，孔隙流体压力与静水压力之差，MPa；％R_o和R_o'为EASY％R_o、A-E_a-R_o模型预测值。

2.按照权利要求1所述一种基于超压抑制有机质成熟度的镜质体反射率预测模型，其特征在于：所述步骤(2)的方法如下：

①准备数据，包括同一深度下对应的超压、温度和EASY％R_o值；

②优选λ₁、λ₂；

利用非线性约束规划求解方法，给λ₁、λ₂任意赋值，输入每组数据，求取每个深度下对应的A-E_a-R_o预测值R_o'及每个预测值的绝对误差，然后求取所有绝对误差的和；调整λ₁、λ₂赋值，使用牛顿迭代方法，重复上述步骤，直至绝对误差和最小，输出λ₁、λ₂值，得到优选出来的值；

％R_on为第n个EASY％R_o值，R_o'为A-E_a-R_o预测值，R_o(m)为实测R_o值。