Method and storage medium for quantitative reconstruction of paleowater depth based on milankovitch cycles

A storage medium which is computer readable and stores thereon a method for quantitative reconstruction of paleowater depth based on Milankovitch cycles, wherein the method comprises:

• selecting a calibration well;
• selecting a rock sample based on the calibration well and acquiring lithofacies sensitive well logging data of the rock sample;
• performing macro-and microelement analysis and assay on the rock sample, and calculating a single-point paleowater depth by utilizing an abundance of the macro-and microelement;
• denoising the lithofacies sensitive well logging data to obtain denoised well logging data;
• performing bandpass filtering analysis and multitaper method (MTM) spectral analysis in a time domain and a depth domain on the denoised well logging data, and determining whether the denoised well logging data accords with a Milankovitch astronomical signal;
• if the denoised well logging data accords with the Milankovitch astronomical signal, performing Evolutive Harmonic Analysis (Eha) & Evolutive Power Spectral Analysis and calculating Evolutive Average Spectral Misfit (eAsm) by Monte Carlo simulation, thereby obtaining a depth-domain spectrum of the denoised well logging data;
• based on Milankovitch astronomical cycle theory, tracking a minimum of a null hypothesis significance level and obtaining a sedimentation rate in the depth-domain spectrum;
• establishing an equal-depth correspondence between the obtained sedimentation rate and the single-point paleowater depth, and establishing and fitting a sedimentation rate-paleowater depth equation;
• verifying the sedimentation rate-paleowater depth equation to obtain a verified sedimentation rate-paleowater depth equation;
• and
• calculating a complete sequence of paleowater depths of the calibration well by the verified sedimentation rate-paleowater depth equation to reconstruct a spatial and temporal distribution of paleowater depths in a target area.