Method and system for optimally selecting carbon storage site based on multi-frequency band seismic data and equipment

A system for optimally selecting a carbon storage site based on multi-frequency band seismic data, comprising:

• well logging instruments configured to generate original seismic data and well logging data of a known well site;
• a data obtaining module configured for obtaining the original seismic data and well logging data of the known well site, then obtaining near-well geological interpretation results and three-dimensional post-stack seismic data, and obtaining the depth data of a maker layer of a target horizon and the isochronous three-dimensional distribution of the marker layer based on the three-dimensional post-stack seismic data;
• a pretreatment module configured for eliminating the abnormal values and performing standardized pretreatment to obtain standardized well logging data based on the well logging data;
• a spread spectrum simulation module configured for performing seismic wavelet spread spectrum simulation to obtain spread spectrum simulated wavelets based on the three-dimensional post-stack seismic data;
• an isochronous stratigraphic framework model obtaining module configured for building an isochronous stratigraphic framework model of a target horizon, and calculating the geometric structure and spatial distribution of the fault-karst based on the spread spectrum simulated wavelets;
• a Bimf component obtaining module configured for performing multi-scale decomposition to obtain multi-scale Bimf components of the seismic data based on the three-dimensional post-stack seismic data and the spread spectrum simulated wavelets;
• a time-depth conversion relationship obtaining module configured for performing well-to-seismic calibration and ignoring small-scale Bimf components in the multi-scale Bimf components of the seismic data to obtain the time-depth conversion relationship based on the spread spectrum simulated wavelets and standardized well logging data;
• a wave impedance inversion data volume obtaining module configured for performing waveform-indicated seismic wave impedance inversion to obtain a wave impedance inversion data volume based on the well logging data, the frequency-boosted seismic data and the time-depth conversion relationship, and on the basis of stratigraphic trend information provided by the isochronous stratigraphic framework model;
• a stable stratum wave impedance data volume obtaining module configured for based on the three-dimensional post-stack seismic data, calculating a three-dimensional variance attribute volume, delineating the a sedimentary stable area, selecting the virtual-well site based on the sedimentary stable area, and obtaining the stable stratum wave impedance data volume through the virtual well cross-well wave impedance interpolation;
• an abnormal wave impedance data volume obtaining module configured for calculating the difference between the stable stratum wave impedance data volume and the wave impedance inversion data volume to obtain the abnormal wave impedance data volume;
• a fault-karst wave impedance data volume obtaining module configured for retaining the abnormal wave impedance data in the spatial geometric contour of the fault-karst by removing areas lower than the average value in the three-dimensional variance attribute volume to obtain the fault-karst wave impedance data volume including the geometric structure and internal wave impedance characteristics of the fault-karst;
• a model interpretation module configured for comparing the near-well geological interpretation result with the fault-karst wave impedance data volume, delineating the characteristic value interval of a cave reservoir, the characteristic value interval of a fractured zone and the characteristic value interval of a surrounding rock, and obtaining a fault-karst reservoir bed interpretation model;
• and
• according to a fault-karst reservoir bed interpretation model, the migration and reservoir space is calculated after carbon dioxide injection, and a trap structure with the preset size threshold value is found, so that the migration and reservoir space where the injected carbon dioxide is stored in a high-density free phase under the conditions of being blocked by the tight carbonate rocks at the side of the fault zone and sealed by a top marl cover coat is selected as the position of the carbon storage box.