[发明专利]一种基于去压实作用批量化计算盆地埋藏史的方法

权利要求书

1.一种基于去压实作用批量化计算盆地埋藏史的方法，其特征在于，该方法具体包括以下步骤：

步骤(1)将多口井的分层和岩性资料按指定表格格式整理到一个表格中；

步骤(2)将步骤(1)整理好的表格按单井进行提取，再分别对每一口单井按单一岩性进行去压实，回剥至地层分界处，最后将所有单井的埋藏史数据整合到一起得到多井埋藏史数据；

步骤(3)选择步骤(2)中多井埋藏史数据中的具体某一地质时期，输出该时期各地层底界在多井中的结果；

步骤(4)对步骤(3)中指定时期、不同地层底界数值按各井的坐标进行插值，绘制不同地层底界的古构造图；

步骤(5)选取步骤(2)中多井数据中具体一口井的埋藏史数据，绘制单井埋藏史图。

2.根据权利要求1所述的一种基于去压实作用批量化计算盆地埋藏史的方法，其特征在于：所述的步骤(1)具体包括以下步骤：

步骤(1.1)将钻井获得的单井编录数据整理到一个数据文档，并且需要将井号放在Excel的第一列，地层代号放在第二列，岩性分层的底界深度放在第三列，岩性放在第四列；

步骤(1.2)将岩性统一归为泥岩、粉砂岩、砂岩和砾岩4大类，并在Excel表中将泥岩类修改为数字1，粉砂岩类修改为数字2，砂岩修改为数字3，砾岩修改为数字4；

步骤(1.3)修改并统一不同井的地层代号名称；

步骤(1.4)筛选无效井，删除单井中只有一套地层分层的井。

3.根据权利要求2所述的一种基于去压实作用批量化计算盆地埋藏史的方法，其特征在于：所述的步骤(2)具体包括以下步骤：

步骤(2.1)按岩性搜集研究区岩石孔隙度和深度数据，计算拟合研究区泥岩、粉砂岩、砂岩和砾岩的压实系数和初始孔隙度；

步骤(2.2)在获得各岩性的压实系数和初始空隙度的基础上，按各岩性层进行去压实校正，恢复到没有压实作用时各岩性层的厚度；

步骤(2.3)根据去压实后的单井数据，首先剥掉第一层岩性单元，再回到步骤(2.2)中去掉第二层的压实作用，再进行剥离，从上至下逐层去压实、剥离每一层岩性单元后，当某一岩性底界对应一岩性地层单元底界时，记录该地层单元的底界深度，为该时期古埋深。

4.根据权利要求3所述的一种基于去压实作用批量化计算盆地埋藏史的方法，其特征在于：所述的步骤(2.1)中采用公式(1)对泥岩、粉砂岩、砂岩、砾岩的孔隙度和深度进行指数关系拟合，

其中，为深度为z时的孔隙度，为初始孔隙度，C为压实系数。

5.根据权利要求4所述的一种基于去压实作用批量化计算盆地埋藏史的方法，其特征在于：所述的步骤(2.2)中，采用单一岩性进行去压实恢复，单一岩性的骨架物质总量为：

其中h_s为深度z₁到z₂时的孔隙度。

6.根据权利要求5所述的一种基于去压实作用批量化计算盆地埋藏史的方法，其特征在于：所述的步骤(2.2)中，将式(1)带入式(2)得到：

上式(3)中，z₂为要求解的岩性底界，z₁为上一套岩性的底界深度z₂，以此类推最上部的第一层岩性单元z₁为井口值0。

7.根据权利要求6所述的一种基于去压实作用批量化计算盆地埋藏史的方法，其特征在于：所述的步骤(2.2)中式(3)是z₂＝f(z₂)的隐函数，通过迭代计算求解出去压实后的每一套岩性底界深度z₂，迭代次数要求不超过100次，相对容差小于10^-5。