[发明专利]一种识别沉积地层中的米兰科维奇周期的方法及系统

权利要求书

1.一种识别沉积地层中的米兰科维奇周期的方法，其特征在于，包括以下步骤：

S10：获得可以反映泥质含量的测井曲线；

S20：对步骤S10中的测井曲线作连续小波变换，得到不同小波尺度对应的小波分解曲线；

S30：统计步骤S20中各小波尺度对应的小波分解曲线中包含的周期个数，作小波尺度与小波分解曲线周期个数对应关系图；

S40：在步骤S30得到的小波尺度与小波分解曲线周期个数对应关系图中，搜索同时符合米兰科维奇不同轨道参数周期之间的周期比值和频率比值的点组；

S50：对步骤S40得到的点组对应的小波分解曲线按照米兰科维奇周期的时间意义得到相应的地层年龄；

S60：引入地质约束条件，将步骤S50中得到的地层年龄与所述地质约束条件对比，确定米兰科维奇周期。

2.权利要求1所述的识别沉积地层中的米兰科维奇周期的方法，其特征在于，步骤S10中采用的测井曲线为自然伽马曲线。

3.根据权利要求1所述的识别沉积地层中的米兰科维奇周期的方法，其特征在于，步骤S20中所述连续小波变换为Morlet连续小波变换，其中Morlet连续小波定义式如下：

其中：Ψ_o为Morlet小波，t为变量。

4.根据权利要求3所述的识别沉积地层中的米兰科维奇周期的方法，其特征在于，步骤S20中所述的连续小波变换，根据下式进行：

其中：

WT_f为f的连续小波变换，

Ψ为小波函数，

t为变量，

a为尺度因子，

b为平移因子。

5.根据权利要求1所述的识别沉积地层中的米兰科维奇周期的方法，其特征在于，步骤S30包括以下子步骤：

S31：对步骤S20中的小波分解曲线作光滑去噪处理，剔除曲线中的扰动噪音；

S32：应用极值点法计算步骤S31中的小波分解曲线中包含的周期个数；

S33：以小波尺度为横坐标，周期个数为纵坐标，将步骤S32中得到的不同小波尺度对应分解曲线包含的周期个数投在图中，得到小波尺度与小波分解曲线周期个数对应关系曲线。

6.根据权利要求1所述的识别沉积地层中的米兰科维奇周期的方法，其特征在于，步骤S40中，在所述小波尺度与小波分解曲线周期个数对应关系图上，以小波尺度比值与所述周期比率对应，周期个数比值与所述频率比率对应。

7.根据权利要求6所述的识别沉积地层中的米兰科维奇周期的方法，其特征在于，所述步骤S40中，所述搜索过程设置合理的误差范围，在误差范围内，均认为搜索到的点组是有效的。

8.根据权利要求1所述的识别沉积地层中的米兰科维奇周期的方法，其特征在于，步骤S50中的地层年龄获得方法为以所述周期对应的周期个数和所述周期对应的时间相乘得到地层年龄。

9.根据权利要求1所述的识别沉积地层中的米兰科维奇周期的方法，其特征在于，步骤S60中的所述的地质约束条件为地层年代格架或井底地层地质年代范围。

10.一种识别沉积地层中的米兰科维奇周期的系统，其特征在于，包括：

曲线获得模块，所述曲线获得模块用于获得可以反映泥质含量的测井曲线；

小波变换模块，所述小波变换模块用于对曲线获得模块中得到的测井曲线作连续小波变换，得到不同小波尺度对应的小波分解曲线；

周期统计模块，所述周期统计模块用于对小波变换模块得到的各小波尺度对应的小波分解曲线中包含的周期个数，作小波尺度与小波分解曲线周期个数对应关系图；

点组搜索模块，所述点组搜索模块用于在周期统计模块得到的小波尺度与小波分解曲线周期个数对应关系图中，搜寻同时符合米兰科维奇不同轨道参数周期之间的周期比值和频率比值的点组；

地层年龄计算模块，所述地层年龄计算模块用于对点组搜索模块得到的点组对应的小波分解曲线按照米兰科维奇周期的时间意义进行定年计算相应的地层年龄；和

米兰科维奇周期确定模块，所述米兰科维奇周期确定模块用于引入地质约束条件，将地层年龄计算模块中得到的地层年龄与所述地质约束条件对比，确定米兰科维奇周期。