CN114966856B - 基于多频带地震资料的碳封存场址优选方法、系统和设备 - Google Patents

基于多频带地震资料的碳封存场址优选方法、系统和设备 Download PDF

Info

Publication number
CN114966856B
CN114966856B CN202210919676.5A CN202210919676A CN114966856B CN 114966856 B CN114966856 B CN 114966856B CN 202210919676 A CN202210919676 A CN 202210919676A CN 114966856 B CN114966856 B CN 114966856B
Authority
CN
China
Prior art keywords
data
seismic
wave impedance
model
matrix
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210919676.5A
Other languages
English (en)
Other versions
CN114966856A (zh
Inventor
田飞
张江云
张旺
郑文浩
单小彩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Institute of Geology and Geophysics of CAS
Original Assignee
Institute of Geology and Geophysics of CAS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Geology and Geophysics of CAS filed Critical Institute of Geology and Geophysics of CAS
Priority to CN202210919676.5A priority Critical patent/CN114966856B/zh
Publication of CN114966856A publication Critical patent/CN114966856A/zh
Priority to US17/902,842 priority patent/US11852771B1/en
Application granted granted Critical
Publication of CN114966856B publication Critical patent/CN114966856B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48—Processing data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/51—Migration
    • G01V2210/514—Post-stack
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/70—Combining sequestration of CO2 and exploitation of hydrocarbons by injecting CO2 or carbonated water in oil wells

Abstract

本发明属于环境监测领域，具体涉及了一种基于多频带地震资料的碳封存场址优选方法、系统和设备，旨在解决断裂带附近发育的断溶体分布复杂，难以准确定位识别断溶体储层的空间几何结构的问题。本发明包括：基于所述三维叠后地震资料，进行地震子波扩频模拟，获得扩频模拟子波；构建目标层位的等时地层格架模型，并计算断溶体几何结构和空间分布；进而进行波形指示反演获得波阻抗反演数据体，通过虚拟井井间波阻抗插值，获取稳定地层波阻抗数据体；计算稳定地层波阻抗数据体与波阻抗反演数据体之差获得波阻抗异常数据体，进而获得断溶体储层解释模型，进而确定碳封存箱的位置。本发明提高了断溶体储层解释模型的精确性，辅助碳封存箱选址。

Description

基于多频带地震资料的碳封存场址优选方法、系统和设备
技术领域
本发明属于二氧化碳的捕获、利用与封存领域，具体涉及了一种基于多频带地震资料的碳封存场址优选方法、系统和设备。
背景技术
为实现全球控制温升的目标，碳捕获、利用与封存（CCUS）技术可将化石能源使用过程产生的二氧化碳封存于地下，从而降低其在大气中的含量，对于减缓全球气候变暖有着十分重要的意义。目前，将二氧化碳封存的方法主要有地质封存、海洋封存、化学封存等。
其中，地质封存是通过把二氧化碳注入到油气藏、深部咸水层、玄武岩、不可采煤层等，使二氧化碳在超临界状态下封存于地下，该方法可在一定条件下实现对二氧化碳的永久性封存。
油气行业封存二氧化碳具有天然优势，利用油气藏封存二氧化碳不仅可以降低大气中二氧化碳的浓度，减缓温室效应，同时还可以通过向正在生产的油气井中注入二氧化碳达到驱油并增加原油产量的目的，具有很好的经济效益。
理论研究和实践表明，当前可进行二氧化碳地质封存的地质体主要有油气藏、深部咸水层、玄武岩、不可采煤层等，其中，由于油气藏自身良好的封闭性，与其他地质体相比，二氧化碳封存于其中的泄漏风险最小，并且由于油气藏中已部署了生产井和注入井，对二氧化碳进行封存也更为方便，经济效益也更好。
碳酸盐岩沉积盆地中，多次构造运动形成纵横交错的走滑断裂带，走滑剪切过程中应力集中，伴生发育了分布密集的裂缝及破碎区域。经大气溶蚀作用与热液溶蚀作用，断层两侧形成了大量形态各异、规模各异的溶洞、裂缝与溶蚀孔洞。在上覆泥灰岩、泥岩等盖层封堵以及侧向致密灰岩遮挡下，形成由不规则缝洞体组成的断溶体圈闭，封闭性良好，是油气资源重要的储集空间同时也为二氧化碳封存提供了天然场所，即碳封存箱。
而断裂带附近发育的碳封存箱分布复杂，几何结构多样，如何准确定位识别碳封存箱的空间几何结构，分析二氧化碳储集空间与封存性能评价是CCUS的主要攻关方向。
为解决这一技术问题，本发明采用测井与地震资料结合的地球物理手段，依据碳封存箱与围岩的地震波阻抗差异性，可实现碳封存箱的定位与结构表征，为二氧化碳储集空间与封存性能评价提供技术支撑。
发明内容
为了解决现有技术中的上述问题，即断裂带附近发育的断溶体分布复杂，几何结构多样，难以准确定位识别碳封存箱的空间几何结构，进而辅助二氧化碳封存的问题，本发明提供了一种基于多频带地震资料的碳封存场址优选方法，所述方法包括：
获取已知井位的原始地震资料和测井资料，进而获取井旁地质解释结果和三维叠后地震资料；基于所述三维叠后地震资料，获取目标层位标志层深度数据和标志层的等时三维展布；
基于所述测井资料，进行剔除异常值和标准化的预处理获得标准化测井数据；
基于所述三维叠后地震资料，进行地震子波扩频模拟，获得扩频模拟子波；
基于所述扩频模拟子波，构建目标层位的等时地层格架模型，并计算断溶体几何结构和空间分布；
基于所述三维叠后地震资料和扩频模拟子波，进行多尺度分解获得地震资料多尺度Bimf分量；
基于所述扩频模拟子波和标准化测井数据，进行井震标定并忽略地震资料多尺度Bimf分量中的小尺度Bimf分量，获得时深转换关系；
基于所述测井资料、提频后的地震资料和时深转换关系，在所述等时地层格架模型提供的地层趋势信息基础上，进行波形指示地震波阻抗反演获得波阻抗反演数据体；
基于所述三维叠后地震资料，计算三维方差属性体，并划定沉积稳定区，基于所述沉积稳定区选择虚拟井井位，并通过虚拟井井间波阻抗插值，获取稳定地层波阻抗数据体；
计算所述稳定地层波阻抗数据体与波阻抗反演数据体之差，获得波阻抗异常数据体，突出了断裂带对稳定沉积地层波阻抗的影响；
通过去除三维方差属性体中低于平均值的区域，保留断溶体的空间几何轮廓内的波阻抗异常数据，获得包含断溶体几何结构与内部波阻抗特征的断溶体波阻抗数据体；
比对所述井旁地质解释结果和断溶体波阻抗数据体，划定洞穴储层特征值区间、裂缝破碎带特征值区间和围岩特征值区间，获得断溶体储层解释模型；
根据所述断溶体储层解释模型，计算二氧化碳注入后的运移储集空间，寻找具有预设尺寸阈值的圈闭构造，使得注入后的二氧化碳在断裂带旁侧致密碳酸盐岩封堵以及顶部泥灰岩盖层密闭条件下以高密度游离相存储的运移储集空间，选定为碳封存箱位置。
在一些优选的实施方式中，所述地震资料多尺度Bimf分量，其获得方法为：
设置初始的叠前道集矩阵：
Figure 928636DEST_PATH_IMAGE001
其中，
Figure 791419DEST_PATH_IMAGE002
为矩阵行坐标，
Figure 436027DEST_PATH_IMAGE003
为矩阵列坐标，
Figure 626837DEST_PATH_IMAGE004
为叠后三维地震资料数据的值，
Figure 257669DEST_PATH_IMAGE005
为矩阵元素值；
设置变量
Figure 850325DEST_PATH_IMAGE006
为分解的Bimf层数，
Figure 793398DEST_PATH_IMAGE006
的初始值为1，定义初始化矩阵变量
Figure 471504DEST_PATH_IMAGE007
，获得变量矩阵
Figure 906027DEST_PATH_IMAGE008
为：
Figure 353189DEST_PATH_IMAGE009
计算关于
Figure 729813DEST_PATH_IMAGE010
的局部最大值矩阵和局部最小值矩阵：
在变量矩阵
Figure 160794DEST_PATH_IMAGE011
中以预设尺寸的观测矩阵进行观测，提取观测矩阵中的元素值，并获取观测矩阵的最大值；
滑动所述观测矩阵，直至观测矩阵中心遍历所有矩阵变量h，将测得的观测矩阵的最大值赋值于对应观测矩阵的中心元素位置，获得局部最大值矩阵
Figure 133429DEST_PATH_IMAGE012
；
通过设置观测矩阵获取观测矩阵中的最小值，进而获得局部最小值矩阵
Figure 700677DEST_PATH_IMAGE013
；
通过样条插值的方法获取最大值包络面和最小值包络面；
插值公式为：
Figure 982622DEST_PATH_IMAGE014
其中，A、B和C表示插值过度矩阵：
Figure 369741DEST_PATH_IMAGE015
Figure 270701DEST_PATH_IMAGE016
Figure 833401DEST_PATH_IMAGE017
其中，
Figure 896035DEST_PATH_IMAGE018
和
Figure 160663DEST_PATH_IMAGE019
为预设的大于0且小于1的值，
Figure 599734DEST_PATH_IMAGE020
为扩频模拟子波；
求得
Figure 875995DEST_PATH_IMAGE021
为最大值包络面，将过度公式B中的局部最大值矩阵
Figure 250476DEST_PATH_IMAGE022
替换为局部最小值矩阵
Figure 612187DEST_PATH_IMAGE023
，获得
Figure 976653DEST_PATH_IMAGE024
为最小值包络面；
计算最大值包络面和最小值包络面的均值，获得中间过度量：
Figure 638579DEST_PATH_IMAGE025
将变量矩阵
Figure 183961DEST_PATH_IMAGE026
减去中间过度量，获得地震资料对应的第j层Bimf分量：
Figure 298547DEST_PATH_IMAGE027
获取余量
Figure 204055DEST_PATH_IMAGE028
作为计算下一层Bimf分量的输入量：
Figure 720487DEST_PATH_IMAGE029
Figure 436771DEST_PATH_IMAGE030
表示第l层的余量，
Figure 507495DEST_PATH_IMAGE031
表示第
Figure 92060DEST_PATH_IMAGE032
层的余量，第
Figure 322053DEST_PATH_IMAGE032
层的余量作为第l层计算Bimf分量的输入量；
迭代计算Bimf分量，直至满足终止条件：
Figure 333871DEST_PATH_IMAGE033
获得地震资料多尺度Bimf分量；其中，r为预设的终止阈值。
在一些优选的实施方式中，所述时深转换关系，其获取方式为：
基于每个已知井位的测井资料中的声波时差曲线和密度曲线做乘积运算获取波阻抗曲线，进而计算反射系数曲线；
以目标层段地震主频为依据构建雷克子波，将雷克子波与反射系数曲线褶积计算后，得到合成地震记录；
将每个钻井井位井眼处的标志层深度数据与标志层的等时三维展布对应，计算合成地震记录与井旁地震道扩频模拟子波的相关性，当波形相关性高于第一相关阈值时，初步井震标定完成，获得测井深度与地震反射波双程旅行时之间的初步时深转化关系：
Figure 157471DEST_PATH_IMAGE034
其中，
Figure 155514DEST_PATH_IMAGE035
表示声波测井标志层深度对应的地震资料的双程旅行时；
Figure 646538DEST_PATH_IMAGE036
为声波时差；
Figure 688312DEST_PATH_IMAGE037
为测井曲线数据采样间隔；
Figure 999208DEST_PATH_IMAGE038
为地震波双程旅行时；
将所述地震资料多尺度Bimf分量以尺度由大到小的顺序，逐步添加Bimf分量获得更新后的有效地震资料；
每添加一个Bimf分量获得更新后的地震资料计算一次更新后的合成地震记录与有效地震资料的第二相关性；
随着有效地震资料逐步添加Bimf分量，第二相关性起初呈现上升趋势，当第二相关性出现下降趋势时，取第二相关性峰值时的更新后的合成地震记录与更新后的地震资料，计算第二时深转化关系；
将所述第二时深转化关系作为最终测得的时深转化关系。
在一些优选的实施方式中，所述波阻抗反演数据体，其获取方式包括：
本发明的另一方面，提出了一种基于多频带地震资料的碳封存场址优选系统，包括：
资料获取模块，配置为获取已知井位的原始地震资料和测井资料，进而获取井旁地质解释结果和三维叠后地震资料；基于所述三维叠后地震资料，获取目标层位标志层深度数据和标志层的等时三维展布；
预处理模块，配置为基于所述测井资料，进行剔除异常值和标准化的预处理获得标准化测井数据；
扩频模拟模块，配置为基于所述三维叠后地震资料，进行地震子波扩频模拟，获得扩频模拟子波；
等时地层格架模型获取模块，配置为基于所述扩频模拟子波，构建目标层位的等时地层格架模型，并计算断溶体几何结构和空间分布；
Bimf分量获取模块，配置为基于所述三维叠后地震资料和扩频模拟子波，进行多尺度分解获得地震资料多尺度Bimf分量；
时深转换关系获取模块，配置为基于所述扩频模拟子波和标准化测井数据，进行井震标定并忽略地震资料多尺度Bimf分量中的小尺度Bimf分量，获得时深转换关系；
波阻抗反演数据体获取模块，配置为基于所述测井资料、提频后的地震资料和时深转换关系，在所述等时地层格架模型提供的地层趋势信息基础上，进行波形指示地震波阻抗反演获得波阻抗反演数据体；
未定地层波阻抗数据体获取模块，配置为基于所述三维叠后地震资料，计算三维方差属性体，并划定沉积稳定区，基于所述沉积稳定区选择虚拟井井位，并通过虚拟井井间波阻抗插值，获取稳定地层波阻抗数据体；
波阻抗异常数据体获取模块，配置为计算所述稳定地层波阻抗数据体与波阻抗反演数据体之差，获得波阻抗异常数据体；
断溶体波阻抗数据体获取模块，配置为通过去除三维方差属性体中低于平均值的区域，保留断溶体的空间几何轮廓内的波阻抗异常数据，获得包含断溶体几何结构与内部波阻抗特征的断溶体波阻抗数据体；
模型解释模块，配置为比对所述井旁地质解释结果和断溶体波阻抗数据体，划定洞穴储层特征值区间、裂缝破碎带特征值区间和围岩特征值区间，获得断溶体储层解释模型；
根据所述断溶体储层解释模型，计算二氧化碳注入后的运移储集空间，寻找具有预设尺寸阈值的圈闭构造，使得注入后的二氧化碳在断裂带旁侧致密碳酸盐岩封堵以及顶部泥灰岩盖层密闭条件下以高密度游离相存储的运移储集空间，选定为碳封存箱位置。
本发明的第三方面，提出了一种电子设备，包括：至少一个处理器；以及与至少一个所述处理器通信连接的存储器；其中，所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的基于多频带地震资料的碳封存场址优选方法。
本发明的第四方面，提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现上述的基于多频带地震资料的碳封存场址优选方法。
本发明的有益效果：
（1）本发明一方面通过Bimf地震资料多尺度分解与井震标定交互方式提升地震资料有效频带，提升了地震资料对于地下地质体细节表征能力。（2）另一方面通过分析稳定地层波阻抗数据体与断裂带波阻抗反演结果突出断溶体与围岩波阻抗差异性，增强了断溶体结构的细节特征。分析沉积稳定区波阻抗与断层，提高了断溶体储层解释模型的精确性，解决了断裂带附近发育的断溶体分布复杂，几何结构多样，难以准确定位识别断溶体储层的空间几何结构，进而辅助二氧化碳封存的问题。（3）在后期进行碳封存工作中，不仅需要考虑储集空间的封存能力，即碳封存箱空间的量化识别。还要从安全角度出发，考察地质体的圈闭构造，即盖层的岩性、厚度是否具有封闭性。该发明根据波阻抗的差异性有效识别并定位了碳封存箱位置，刻画碳封存箱几何特征及其发育环境。有助于后期对碳封存储量评估的工作。（4）由于断溶体内部的不均一性，二氧化碳储存空间难以定量评价。本发明所获得的断溶体储层解释模型不仅具有断溶体几何结构特征，而且还具有丰富的内部波阻抗结构信息，对于评价断溶体内部有效二氧化碳储集空间具有重要意义。以此结论作为碳封存箱选址依据可靠性较强。
附图说明
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述，本申请的其它特征、目的和优点将会变得更明显：
图1是本发明实施例中基于多频带地震资料的碳封存场址优选方法的流程示意图；
图2是本发明实施例中原始地震资料示意图；
图3是本发明实施例中等时地层格架模型的示意图；
图4是本发明实施例中以r取值为0.2时地震资料多尺度Bimf分量示意图；
图5为本发明实施例中虚拟井井位网格示意图；
图6为本发明实施例中稳定地层波阻抗数据体的效果示意图；
图7为本发明实施例中波阻抗异常数据体效果示意图；
图8为本发明实施例中二氧化碳封存选址结果与封存箱内部几何结构示意图。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅用于解释相关发明，而非对该发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与有关发明相关的部分。
需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
本发明提供一种基于多频带地震资料的碳封存场址优选方法，本方法包括：
本发明的一种基于多频带地震资料的碳封存场址优选方法，通过将所述背景地层波阻抗模型与波阻抗反演数据体进行做差运算的方式，断裂带附近发育的碳封存箱分布复杂，几何结构多样，难以准确定位识别碳封存箱的空间几何结构，进而辅助二氧化碳封存的问题。
为了更清晰地对本发明基于多频带地震资料的碳封存场址优选方法进行说明，下面结合图1对本发明实施例中各步骤展开详述。
本发明第一实施例的基于多频带地震资料的碳封存场址优选方法，各步骤详细描述如下：
获取已知井位的原始地震资料和测井资料，进而获取井旁地质解释结果和三维叠后地震资料；基于所述三维叠后地震资料，获取目标层位标志层深度数据和标志层的等时三维展布；原始地震资料如图2所示；
过观察地震波形数据确定鹰山组顶界标志层：与上覆巴楚组泥岩层相比，碳酸盐岩内幕区波形反射不规则，无一定方向，振幅可强可弱，同相轴可长可短连续性差；且具有非系统性同相轴反射终止和分叉现象。根据二者分界面强轴确定鹰山组顶面，作为标志层。
所述测井资料为借助常规测井仪器，在工区井眼处探测深度范围为5500~5750m的九条常规测井曲线数据；测井曲线数据采样间隔为0.01m。借助成像测井、钻井、录井、岩心等其他信息，得到个别深度段的确定的岩性、物性信息，进而划定围岩、洞穴型储层、裂缝型储层、孔洞型储层的各项参数阈值范围，根据各项参数阈值范围获得井旁地质解释结果。
本实施例运用三维地震勘探方法，借助地震波激发源与地震信号检波器，得到工区面积约27km2的三维叠后地震资料数据，信号记录双程旅行时为4s，采样点时间间隔为1ms，探测深度超过6000m。
本实施例用于5m以上的目标体探测，对地震资料主频要求较高，应在50~60Hz范围内。当缝洞体发育层段内地震资料振幅数据体主频低于50HZ时，则需要采用基于地震子波模拟的一维经验模态分解算法进行地震数据的拓频降噪处理，得到高分辨率、高信噪比的扩频地震资料。
基于所述测井资料，进行剔除异常值和标准化的预处理获得标准化测井数据；
所述异常值即为离群点，基于所得九种常规测井曲线数据，将九种常规测井曲线数据中每种曲线的所有井眼处数据绘制统计直方图，合理调整区间阈值，去除与中位数偏离最大的前百分之5的数据，保留框选区间内的数据点，即为剔除异常值的测井资料。
所述标准化处理，基于剔除异常值的测井资料，叠合工区内所有已知井位的单种剔除异常值的测井数据，绘制曲线柱状图，通过整合阈值获取标准化测井数据。
剔除异常值和标准化的步骤，可以消除井间仪器带来的影响。
基于所述三维叠后地震资料，进行地震子波扩频模拟，获得扩频模拟子波；地震子波模拟是一种在确保处理后的地震数据具有较高保真度高的前提下，通过拓宽有效频带从而提高地震信号分辨率的数据处理方法。
在本实施例中，所述扩频模拟子波，其获得方法为：
将原始地震资料的地震记录褶积模型在频率域的表示为：
Figure 659996DEST_PATH_IMAGE039
其中，
Figure 880893DEST_PATH_IMAGE040
表示傅氏变换后的地震记录频谱，
Figure 234514DEST_PATH_IMAGE041
表示傅氏变换后的地震子波频谱，
Figure 160269DEST_PATH_IMAGE042
表示傅氏变换后的反射系数频谱，
Figure 624749DEST_PATH_IMAGE043
表示角频率；
将所述地震记录褶积模型在频率域的表示转化为地震记录线性系统频域的表示：
Figure 700152DEST_PATH_IMAGE044
其中，
Figure 959095DEST_PATH_IMAGE045
为地震记录线性系统频域表示，
Figure 979004DEST_PATH_IMAGE046
为地震子波线性系统频域表示，
Figure 637387DEST_PATH_IMAGE047
为反射系数线性系统频域表示；
反傅氏变换为：
Figure 691931DEST_PATH_IMAGE048
其中，
Figure 121775DEST_PATH_IMAGE049
表示地震记录频谱的复赛谱序列，
Figure 769925DEST_PATH_IMAGE050
表示地震子波频谱的复赛谱序列，
Figure 841786DEST_PATH_IMAGE051
表示反射系数频谱的复赛谱序列，
Figure 609891DEST_PATH_IMAGE052
表示地震波到时；
通过低通滤波器将复赛谱中的子波复赛谱序列和反射系数复赛谱序列进行分离，提取子波振幅；子波振幅得到提取利用了子波与反射系数序列平滑程度的差异在复赛谱中易于区分的特性，子波能量出现于原点附近，而反射系数序列则远离原点；
基于所述子波振幅，通过最小二乘法模拟地震子波振幅谱：
Figure 741795DEST_PATH_IMAGE053
其中，
Figure 470717DEST_PATH_IMAGE054
表示大于等于0的常数，
Figure 221635DEST_PATH_IMAGE055
表示地震子波频谱，由
Figure 985192DEST_PATH_IMAGE056
经傅氏变换而来，
Figure 412631DEST_PATH_IMAGE057
和
Figure 894428DEST_PATH_IMAGE058
为待求的关于频率
Figure 917879DEST_PATH_IMAGE059
的多项式，用于拟合地震子波振幅谱；
基于所述模拟地震子波振幅谱，获得子波最大相位分量和最小相位分量；
设子波
Figure 801521DEST_PATH_IMAGE060
的最大相位分量为
Figure 396932DEST_PATH_IMAGE061
、最小相位分量为
Figure 366025DEST_PATH_IMAGE062
，则子波
Figure 317800DEST_PATH_IMAGE063
为：
Figure 931315DEST_PATH_IMAGE064
振幅谱的复赛谱中表示为：
Figure 575923DEST_PATH_IMAGE065
其中，振幅谱的复赛谱
Figure 891367DEST_PATH_IMAGE066
在复赛谱的正、负轴上对称显示，
Figure 646833DEST_PATH_IMAGE067
为地震子波最小相位分量
Figure 239489DEST_PATH_IMAGE068
所对应的最小相位函数的复赛谱，
Figure 664785DEST_PATH_IMAGE069
为地震子波最大相位分量 QUOTE
Figure 608470DEST_PATH_IMAGE070
Figure 26682DEST_PATH_IMAGE070
所对应的最小相位函数的复赛谱；
基于所述振幅谱的复赛谱确定一组具有相同振幅谱的混合相位子波集合，调整俞氏子波参数，在保证地震子波主频完整的前提下，保持低频、扩展高频提升主频，提升有效带宽至预设的带宽阈值，获得扩频模拟子波。通常设置有效贷款为0~60Hz。
基于所述扩频模拟子波，构建目标层位的等时地层格架模型，并计算碳封存箱几何结构和空间分布；等时地层格架模型如图3所示；
基于所述三维叠后地震资料和扩频模拟子波，进行多尺度分解获得地震资料多尺度Bimf分量；本实施例通过将扩频模拟子波的逐级筛选，实现对地震信号的多尺度分解，能够快速提取不同尺度的特征信号，分离噪声和有用信号，实现地震信号内部趋势信息的提取。进而提升了地震资料对于地下地质体细节表征能力。本发明的目的是要得到高分辨率的结果，高分辨率意味着需要高频带的地震数据，但是现在有效频带的分布范围为20-35Hz，在高频带中存在大量的噪声，影响对断溶体的观测，通过忽略小尺度的Bimf分量，可避免高分辨率中噪声的影响，提高准确率。
在本实施例中，所述地震资料多尺度Bimf分量，其获得方法为：
将三维叠后地震资料沿“道集”方向切片，获得若干个叠前道集矩阵；
Figure 5002DEST_PATH_IMAGE071
。
通过EMD分解的方式获取Bimf分量：
设置初始的叠前道集矩阵：
Figure 601200DEST_PATH_IMAGE001
其中，
Figure 422394DEST_PATH_IMAGE002
为矩阵行坐标，
Figure 254084DEST_PATH_IMAGE003
为矩阵列坐标，
Figure 86911DEST_PATH_IMAGE004
为叠后三维地震资料数据的值，
Figure 119589DEST_PATH_IMAGE005
为矩阵元素值；
设置变量
Figure 772287DEST_PATH_IMAGE006
为分解的Bimf层数，
Figure 800811DEST_PATH_IMAGE006
的初始值为1，定义初始化矩阵变量
Figure 222565DEST_PATH_IMAGE007
，获得变量矩阵
Figure 691723DEST_PATH_IMAGE008
为：
Figure 831718DEST_PATH_IMAGE009
计算关于
Figure 395423DEST_PATH_IMAGE010
的局部最大值矩阵和局部最小值矩阵：
在变量矩阵
Figure 671683DEST_PATH_IMAGE010
中以预设尺寸的观测矩阵进行观测，提取观测矩阵中的元素值，并获取观测矩阵的最大值；所述预设尺寸可选用
Figure 46164DEST_PATH_IMAGE072
的小矩阵进行观测；
滑动所述观测矩阵，直至观测矩阵中心遍历所有矩阵变量h，将测得的观测矩阵的最大值赋值于对应观测矩阵的中心元素位置，获得局部最大值矩阵
Figure 673455DEST_PATH_IMAGE012
；
通过设置观测矩阵获取观测矩阵中的最小值，进而获得局部最小值矩阵
Figure 40851DEST_PATH_IMAGE013
；
通过样条插值的方法获取最大值包络面和最小值包络面；
插值公式为：
Figure 437197DEST_PATH_IMAGE014
其中，A、B和C表示插值过度矩阵：
Figure 513738DEST_PATH_IMAGE073
Figure 752958DEST_PATH_IMAGE074
Figure 799411DEST_PATH_IMAGE017
其中，
Figure 925630DEST_PATH_IMAGE018
和
Figure 500968DEST_PATH_IMAGE019
为预设的大于0且小于1的值，
Figure 102851DEST_PATH_IMAGE020
为扩频模拟子波；
求得
Figure 832557DEST_PATH_IMAGE021
为最大值包络面，将过度公式B中的局部最大值矩阵
Figure 937917DEST_PATH_IMAGE022
替换为局部最小值矩阵
Figure 825101DEST_PATH_IMAGE023
，获得
Figure 648701DEST_PATH_IMAGE024
为最小值包络面；
计算最大值包络面和最小值包络面的均值，获得中间过度量：
Figure 771378DEST_PATH_IMAGE025
将变量矩阵
Figure 121456DEST_PATH_IMAGE026
减去中间过度量，获得地震资料对应的第j层Bimf分量：
Figure 304176DEST_PATH_IMAGE027
获取余量
Figure 490438DEST_PATH_IMAGE028
作为计算下一层Bimf分量的输入量：
Figure 416806DEST_PATH_IMAGE029
Figure 886970DEST_PATH_IMAGE030
表示第l层的余量，
Figure 240591DEST_PATH_IMAGE031
表示第
Figure 773203DEST_PATH_IMAGE032
层的余量，第
Figure 378628DEST_PATH_IMAGE032
层的余量作为第l层Bimf分量的输入量；
迭代计算Bimf分量，直至满足终止条件：
Figure 968878DEST_PATH_IMAGE033
获得地震资料多尺度Bimf分量；其中，r为预设的终止阈值；r通常取值为0.2~0.3，本实施例以0.2为例，如图4所示，能够保证Bimf的数目和质量，保证其能够更好的反映波形的细节。
小尺度Bimf分量包含了大量的噪声和信号的细节信息以及边缘信息，而中大尺度Bimf分量则表征了信号的内部结构特征和信号的趋势特征。
通过忽略地震资料小尺度Bimf分量，发现合成地震记录与井旁地震到相关系数先增加后降低，选取相关系数最高处对应的削减后的地震资料进行后续运算。
经过子波扩频模拟后，拓展了地震数据的有效频带，高频部分得到了合理加强。在地震波形上表现为同相轴数目增多，更容易反映地震波反射信息的细节变化，并且在振幅、相位、频率方面改善了同一反射波组波形的一致性。在缝洞体地震响应上，“串珠状”反射特征尤为明显，且串珠内部形态的细节能够清晰地显示，代表了不同结构特征与充填物组合的复杂缝洞型储层地震反射，有助于后期精细地质解释。
基于所述扩频模拟子波和标准化测井数据，进行井震标定并忽略地震资料多尺度Bimf分量中的小尺度Bimf分量，获得时深转换关系；
基于所述测井资料、提频后的地震资料和时深转换关系，在所述等时地层格架模型提供的地层趋势信息基础上，进行波形指示地震波阻抗反演获得波阻抗反演数据体；
在本实施例中，所述时深转换关系，其获取方式为：
基于每个已知井位的测井资料中的声波时差曲线和密度曲线做乘积运算获取波阻抗曲线，进而计算反射系数曲线；
以目标层段地震主频为依据构建雷克子波，优选采用25Hz的雷克子波，将雷克子波与反射系数曲线褶积计算后，得到合成地震记录；
将每个钻井井位井眼处的标志层深度数据与标志层的等时三维展布对应，计算合成地震记录与井旁地震道扩频模拟子波的相关性，当波形相关性高于第一相关阈值时，第一相关阈值优选大于85%，初步井震标定完成，获得测井深度与地震反射波双程旅行时之间的初步时深转化关系：
Figure 227821DEST_PATH_IMAGE034
其中，
Figure 513309DEST_PATH_IMAGE035
表示声波测井标志层深度对应的地震资料的双程旅行时；
Figure 656846DEST_PATH_IMAGE036
为声波时差；
Figure 711389DEST_PATH_IMAGE037
为测井曲线数据采样间隔；
Figure 534377DEST_PATH_IMAGE038
为地震波双程旅行时；
将所述地震资料多尺度Bimf分量以尺度由大到小的顺序，逐步添加Bimf分量获得更新后的有效地震资料；
每添加一个Bimf分量获得更新后的有效地震资料计算一次更新后的合成地震记录与有效地震资料的第二相关性；
随着有效地震资料逐步添加Bimf分量，第二相关性起初呈现上升趋势，当第二相关性出现下降趋势时，取第二相关性峰值时的更新后的合成地震记录与更新后的地震，计算第二时深转化关系；
将所述第二时深转化关系作为最终测得的时深转化关系。
在本实施例中，所述波阻抗反演数据体，其获取方式包括：
基于所述扩频模拟子波，计算待判别地震道波形与已知井的合成地震记录的波形相关性，根据波形相关性最高的井对应的波阻抗曲线建立初始模型；
利用白噪声满足高斯分布的规律，将测井资料中的波阻抗数据表示为：
Figure 307160DEST_PATH_IMAGE075
其中，
Figure 988809DEST_PATH_IMAGE076
表示测井波阻抗曲线，
Figure 163438DEST_PATH_IMAGE077
表示待求解的地下地层实际波阻抗值，
Figure 764184DEST_PATH_IMAGE078
表示随机噪声；
根据中心极限定理，
Figure 148897DEST_PATH_IMAGE079
也满足高斯分布，确定初始目标函数为：
Figure 899816DEST_PATH_IMAGE080
其中，
Figure 928952DEST_PATH_IMAGE081
表示与后验信息有关的函数，
Figure 356391DEST_PATH_IMAGE082
表示基于最优样本数选取样本井对样本井的波阻抗曲线，进行匹配滤波后，求得后验概率统计分布密度，进而计算得到的波阻抗期望值，
Figure 838188DEST_PATH_IMAGE083
表示白噪声的协方差；
基于所述初始目标函数，通过最大后验估计，在目标函数中引入先验信息，获得稳定的目标函数为：
Figure 127218DEST_PATH_IMAGE084
其中，
Figure 10860DEST_PATH_IMAGE085
表示待模拟的特征参数，
Figure 609201DEST_PATH_IMAGE086
表示与地质和测井资料先验信息有关的函数，
Figure 578294DEST_PATH_IMAGE087
表示用于协调
Figure 405435DEST_PATH_IMAGE081
和
Figure 143584DEST_PATH_IMAGE086
之间的相互影响的平滑参数；
以所述稳定的目标函数作为初始模型的输入，通过马尔科夫链蒙特卡罗方法MCMC和Metropolis-Hastings抽样准则对后验概率分布抽样，不断优化初始模型的参数，选取目标函数取最大值时的解作为随机实现，取多次随机实现的均值作为期望值输出，将所述期望值输出作为波阻抗反演数据体。
基于所述三维叠后地震资料，计算三维方差属性体，并划定沉积稳定区，基于所述沉积稳定区选择虚拟井井位，并通过虚拟井井间波阻抗插值，获取稳定地层波阻抗数据体；虚拟经井位划分如图5所示；
在本实施例中，所述划定沉积稳定区，具体为：
设三维叠后地震资料中的各采样点数据为
Figure 788192DEST_PATH_IMAGE088
，其中
Figure 366285DEST_PATH_IMAGE089
表示地震测网线号，
Figure 121752DEST_PATH_IMAGE090
代表地震测网道号，
Figure 589773DEST_PATH_IMAGE091
代表1ms采样的地震记录采样点序号；
计算预设采样区域的采样点数据均方差：
Figure 405283DEST_PATH_IMAGE092
将采样区域进行纵向和横向平移，遍历计算所有采样区域的数据均方差，获得三维方差属性体；
对所述三维方差属性体进行切片，获取方差属性数据在平面上的分布特征，将方差属性值低于总体能量均值的区域作为沉积稳定区。
在本实施例中，所述稳定地层波阻抗数据体，其具体获得方式为：
基于所述沉积稳定区，划定预设尺寸的网格，将每个网格节点视为虚拟井井位；本实施例将沉积稳定区结施图建立50*50m的网格，将每个网格节点视为虚拟井井位。
设未知地下波阻抗模型参数m为
Figure 208022DEST_PATH_IMAGE093
维空间向量
Figure 767180DEST_PATH_IMAGE094
；
通过一次地震采集过程获得了观测数据d为k维数据向量
Figure 355287DEST_PATH_IMAGE095
；
通过非线性函数核G建立未知地下波阻抗模型参数和k维数据向量的联系，即获得正演模型：
Figure 810539DEST_PATH_IMAGE096
其中，
Figure 241520DEST_PATH_IMAGE097
表示与地下波阻抗模型参数m无关的随机噪声，服从高斯分布；
基于正演模型构建反演目标函数：
Figure 463423DEST_PATH_IMAGE098
其中，F为观测数据d与预测数据G(m)之间的均方误差；
观测数据d和待反演模型参数m之间存在严重的非线性，因此需要对反演目标函数进行线性化求解，解的精度依赖于先验信息的可靠性。
对所述反演目标函数进行线性化求解；
对所述正演模型进行泰勒展开并略去二次以上的高阶项，获得预测数据简略表达：
Figure 30671DEST_PATH_IMAGE099
其中，
Figure 797770DEST_PATH_IMAGE100
表示依据先验信息建立的初始模型，A为Jacobian矩阵，A的元素为一阶偏微分
Figure 716047DEST_PATH_IMAGE101
；
设
Figure 351428DEST_PATH_IMAGE102
，
Figure 897816DEST_PATH_IMAGE103
，
Figure 226029DEST_PATH_IMAGE104
表示随机噪声，则正演模型的迭代方程为：
Figure 241389DEST_PATH_IMAGE105
其中，
Figure 946040DEST_PATH_IMAGE106
表示迭代了
Figure 615443DEST_PATH_IMAGE107
次后的正演模型，
Figure 114558DEST_PATH_IMAGE108
表示迭代了k+1次后的正演模型，
Figure 617215DEST_PATH_IMAGE109
为加噪预测数据；
显然，该线性化反演方法只有在初始模型较为接近真实模型时才会得到精确的结果，且从数学计算上来讲，该线性方程组往往趋于病态，稳定性较差。当先验信息不足时，可尝试对待反演模型参数进行量化约束以控制解空间大小，以期得到稳定精确的反演结果，这就是宽带约束反演方法的基本思想，其本质是利用最大似然分布获取最优解。
通过一次地震采集过程获得叠前深度域偏移地震剖面
Figure 859977DEST_PATH_IMAGE110
，
假设反演深度域模型参数有先验概率分布
Figure 115378DEST_PATH_IMAGE111
，根据贝叶斯公式有概率分布为：
Figure 785394DEST_PATH_IMAGE112
从概率论的角度来说，反演的目的是获得最大的后验概率密度P(M=m|D=d)，从地球物理的角度来说，就是在已知深度偏移的情况下，存在多种深度域模型参数能通过正演过程形成深度域偏移剖面，解空间很大，反演的目的则为寻找出最接近于真实地球模型的解
加入随机噪音的概率分布为：
Figure 775346DEST_PATH_IMAGE113
其中，
Figure 556221DEST_PATH_IMAGE114
为噪音协方差矩阵，
Figure 807073DEST_PATH_IMAGE115
为噪音协方差矩阵求逆运算，det为取对应矩阵的行列式；
n为噪音且服从高斯分布均值为0，
Figure 772624DEST_PATH_IMAGE116
为一等值的对角矩阵且对角线元素为数据中心包含噪音的方差
Figure 374507DEST_PATH_IMAGE117
，则加入随机噪音的概率分布变形为：
Figure 568859DEST_PATH_IMAGE118
对于反演深度域模型参数先验概率分布
Figure 939797DEST_PATH_IMAGE119
，设置初始模型
Figure 341829DEST_PATH_IMAGE120
，且有
Figure 165428DEST_PATH_IMAGE121
，等价的深度域模型概率分布为：
Figure 429050DEST_PATH_IMAGE122
；
设等价的深度域模型概率分布服从高斯分布，有加入随机噪声的深度域模型概率分布为：
Figure 654495DEST_PATH_IMAGE123
假设深度域模型参数均值与初始模型
Figure 693340DEST_PATH_IMAGE124
等值，且互不影响，有简化的深度域模型概率分布：
Figure 269815DEST_PATH_IMAGE125
其中，
Figure 71549DEST_PATH_IMAGE126
表示真实模型相对于所给定初始模型的方差值；
加入随机噪音的概率分布变形中的分母
Figure 885921DEST_PATH_IMAGE127
依赖于数据采集和处理，设为常数，通过后验概率最大化等价于分子最大化，使后验概率分布式最大：
Figure 505121DEST_PATH_IMAGE128
使后验概率分布式最大等价式为：
Figure 427947DEST_PATH_IMAGE129
令关于模型参数m的偏导数等于零，且
Figure 892426DEST_PATH_IMAGE131
，
Figure 702250DEST_PATH_IMAGE132
，获得矩阵方程：
Figure 492352DEST_PATH_IMAGE133
Figure 902473DEST_PATH_IMAGE134
其中，
Figure 905064DEST_PATH_IMAGE135
，I为单位矩阵，
Figure 959608DEST_PATH_IMAGE136
为波阻抗迭代矩阵，整理得到随机逆反演的基本公式：
Figure 795977DEST_PATH_IMAGE137
Figure 303182DEST_PATH_IMAGE138
随机逆反演的基本公式的迭代公式为：
Figure 234097DEST_PATH_IMAGE139
；
其中，
Figure 674306DEST_PATH_IMAGE140
为由参数模型所形成的深度域合成地震记录；
通过不断迭代更新未知地下波阻抗模型的参数，直至
Figure 415997DEST_PATH_IMAGE141
出现正负值震荡，此时的
Figure 676077DEST_PATH_IMAGE142
为低频波阻抗反演数据；
基于虚拟井井位平面坐标与地震数据平面坐标的相关关系，进而确定虚拟井与所述低频波阻抗反演数据的一一对应关系，进而将虚拟井井旁低频波阻抗反演数据对虚拟井赋值，获得虚拟井井位的波阻抗数据；
基于波阻抗反演数据体提取虚拟井井位波阻抗数据，根据地层格架限定的时窗范围针对所有虚拟井波阻抗数据的插值计算，设定计算区域为地质格架内工区的全部范围，得到稳定地层波阻抗数据；稳定地层波阻抗数据效果如图6所示；
计算所述稳定地层波阻抗数据体与波阻抗反演数据体之差，获得波阻抗异常数据体；波阻抗异常数据体效果如图7所示；
通过去除三维方差属性体中低于平均值的区域，保留断溶体的空间几何轮廓内的波阻抗异常数据，获得包含断溶体几何结构与内部波阻抗特征的断溶体波阻抗数据体；
比对所述井旁地质解释结果和断溶体波阻抗数据体，划定洞穴储层特征值区间、裂缝破碎带特征值区间和围岩特征值区间，获得断溶体储层解释模型；
根据所述断溶体储层解释模型，评价盖层的岩性、厚度选定碳封存箱位置。
输入测井解释成果比如在钻井不通深度段的岩性解释，包括裂缝储层、孔洞储层、过度带、围岩；根据井位、测井数据与地震资料的时深转化关系将钻井轨迹、测井解释成果投影在断溶体波阻抗数据剖面图中。
比对测井解释结果与断溶体结构-特征值模型的特征值能量数据，划定特征值大于0.82的区域为裂缝储层，划定特征值处于0.63~0.82的区域为孔洞储层，划定特征值处于0.31~0.63的区域为过渡带，划定特征值小于0.31的区域为围岩，以此作为最终的断溶体储层解释模型；
根据所述断溶体储层解释模型，评价盖层的岩性、厚度选定碳封存箱位置。二氧化碳封存选址结果与封存箱内部几何结构如图8所示。
上述实施例中虽然将各个步骤按照上述先后次序的方式进行了描述，但是本领域技术人员可以理解，为了实现本实施例的效果，不同的步骤之间不必按照这样的次序执行，其可以同时(并行)执行或以颠倒的次序执行，这些简单的变化都在本发明的保护范围之内。
本发明第二实施例的基于多频带地震资料的碳封存场址优选系统，所述系统包括：
资料获取模块，配置为获取已知井位的原始地震资料和测井资料，进而获取井旁地质解释结果和三维叠后地震资料；基于所述三维叠后地震资料，获取目标层位标志层深度数据和标志层的等时三维展布；
预处理模块，配置为基于所述测井资料，进行剔除异常值和标准化的预处理获得标准化测井数据；
扩频模拟模块，配置为基于所述三维叠后地震资料，进行地震子波扩频模拟，获得扩频模拟子波；
等时地层格架模型获取模块，配置为基于所述扩频模拟子波，构建目标层位的等时地层格架模型，并计算断溶体几何结构和空间分布；
Bimf分量获取模块，配置为基于所述三维叠后地震资料和扩频模拟子波，进行多尺度分解获得地震资料多尺度Bimf分量；
时深转换关系获取模块，配置为基于所述扩频模拟子波和标准化测井数据，进行井震标定并忽略地震资料多尺度Bimf分量中的小尺度Bimf分量，获得时深转换关系；
波阻抗反演数据体获取模块，配置为基于所述测井资料、提频后的地震资料和时深转换关系，在所述等时地层格架模型提供的地层趋势信息基础上，进行波形指示地震波阻抗反演获得波阻抗反演数据体；
未定地层波阻抗数据体获取模块，配置为基于所述三维叠后地震资料，计算三维方差属性体，并划定沉积稳定区，基于所述沉积稳定区选择虚拟井井位，并通过虚拟井井间波阻抗插值，获取稳定地层波阻抗数据体；
波阻抗异常数据体获取模块，配置为计算所述稳定地层波阻抗数据体与波阻抗反演数据体之差，获得波阻抗异常数据体；
断溶体波阻抗数据体获取模块，配置为通过去除三维方差属性体中低于平均值的区域，保留断溶体的空间几何轮廓内的波阻抗异常数据，获得包含断溶体几何结构与内部波阻抗特征的断溶体波阻抗数据体；
模型解释模块，配置为比对所述井旁地质解释结果和断溶体波阻抗数据体，划定洞穴储层特征值区间、裂缝破碎带特征值区间和围岩特征值区间，获得断溶体储层解释模型；
根据所述断溶体储层解释模型，计算二氧化碳注入后的运移储集空间，寻找具有预设尺寸阈值的圈闭构造，使得注入后的二氧化碳在断裂带旁侧致密碳酸盐岩封堵以及顶部泥灰岩盖层密闭条件下以高密度游离相存储的运移储集空间，选定为碳封存箱位置。所属技术领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的系统的具体工作过程及有关说明，可以参考前述方法实施例中的对应过程，在此不再赘述。
需要说明的是，上述实施例提供的基于多频带地震资料的碳封存场址优选系统，仅以上述各功能模块的划分进行举例说明，在实际应用中，可以根据需要而将上述功能分配由不同的功能模块来完成，即将本发明实施例中的模块或者步骤再分解或者组合，例如，上述实施例的模块可以合并为一个模块，也可以进一步拆分成多个子模块，以完成以上描述的全部或者部分功能。对于本发明实施例中涉及的模块、步骤的名称，仅仅是为了区分各个模块或者步骤，不视为对本发明的不当限定。
本发明第三实施例的一种电子设备，包括：至少一个处理器；以及与至少一个所述处理器通信连接的存储器；其中，所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现上述的基于多频带地震资料的碳封存场址优选方法。
本发明第四实施例的一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现上述的基于多频带地震资料的碳封存场址优选方法。
所属技术领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的存储装置、处理装置的具体工作过程及有关说明，可以参考前述方法实施例中的对应过程，在此不再赘述。
术语“第一”、“第二”等是用于区别类似的对象，而不是用于描述或表示特定的顺序或先后次序。
术语“包括”或者任何其它类似用语旨在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备/装置不仅包括那些要素，而且还包括没有明确列出的其它要素，或者还包括这些过程、方法、物品或者设备/装置所固有的要素。
至此，已经结合附图所示的优选实施方式描述了本发明的技术方案，但是，本领域技术人员容易理解的是，本发明的保护范围显然不局限于这些具体实施方式。在不偏离本发明的原理的前提下，本领域技术人员可以对相关技术特征做出等同的更改或替换，这些更改或替换之后的技术方案都将落入本发明的保护范围之内。

Claims (10)

1.一种基于多频带地震资料的碳封存场址优选方法，其特征在于，所述方法包括：
获取已知井位的原始地震资料和测井资料，进而获取井旁地质解释结果和三维叠后地震资料；基于所述三维叠后地震资料，获取目标层位标志层深度数据和标志层的等时三维展布；
基于所述测井资料，进行剔除异常值和标准化的预处理获得标准化测井数据；
基于所述三维叠后地震资料，进行地震子波扩频模拟，获得扩频模拟子波；
基于所述扩频模拟子波，构建目标层位的等时地层格架模型，并计算断溶体几何结构和空间分布；
基于所述三维叠后地震资料和扩频模拟子波，进行多尺度分解获得地震资料多尺度Bimf分量；
基于所述扩频模拟子波和标准化测井数据，进行井震标定并忽略地震资料多尺度Bimf分量中的小尺度Bimf分量，获得时深转换关系；
基于所述测井资料、提频后的地震资料和时深转换关系，在所述等时地层格架模型提供的地层趋势信息基础上，进行波形指示地震波阻抗反演获得波阻抗反演数据体；
基于所述三维叠后地震资料，计算三维方差属性体，并划定沉积稳定区，基于所述沉积稳定区选择虚拟井井位，并通过虚拟井井间波阻抗插值，获取稳定地层波阻抗数据体；
计算所述稳定地层波阻抗数据体与波阻抗反演数据体之差，获得波阻抗异常数据体；
通过去除三维方差属性体中低于平均值的区域，保留断溶体的空间几何轮廓内的波阻抗异常数据，获得包含断溶体几何结构与内部波阻抗特征的断溶体波阻抗数据体；
比对所述井旁地质解释结果和断溶体波阻抗数据体，划定洞穴储层特征值区间、裂缝破碎带特征值区间和围岩特征值区间，获得断溶体储层解释模型；
根据所述断溶体储层解释模型，计算二氧化碳注入后的运移储集空间，寻找具有预设尺寸阈值的圈闭构造，使得注入后的二氧化碳在断裂带旁侧致密碳酸盐岩封堵以及顶部泥灰岩盖层密闭条件下以高密度游离相存储的运移储集空间，选定为碳封存箱位置。
2.根据权利要求1所述的基于多频带地震资料的碳封存场址优选方法，其特征在于，所述扩频模拟子波，其获得方法为：
将原始地震资料的地震记录褶积模型在频率域的表示为：
Figure 978309DEST_PATH_IMAGE001
其中，
Figure 269613DEST_PATH_IMAGE002
表示傅氏变换后的地震记录频谱，
Figure 271067DEST_PATH_IMAGE003
表示傅氏变换后的地震子波频谱，
Figure 938809DEST_PATH_IMAGE004
表示傅氏变换后的反射系数频谱，
Figure 342108DEST_PATH_IMAGE005
表示角频率；
将所述地震记录褶积模型在频率域的表示转化为地震记录线性系统频域的表示：
Figure 56511DEST_PATH_IMAGE006
其中，
Figure 545261DEST_PATH_IMAGE007
为地震记录线性系统频域表示，
Figure 751115DEST_PATH_IMAGE008
为地震子波线性系统频域表示，
Figure 8921DEST_PATH_IMAGE009
为反射系数线性系统频域表示；
反傅氏变换为：
Figure 907606DEST_PATH_IMAGE010
其中，
Figure 132920DEST_PATH_IMAGE011
表示地震记录频谱的复赛谱序列，
Figure 142465DEST_PATH_IMAGE012
表示地震子波频谱的复赛谱序列，
Figure 254777DEST_PATH_IMAGE013
表示反射系数频谱的复赛谱序列，t表示地震波到时；
通过低通滤波器将复赛谱中的子波复赛谱序列和反射系数复赛谱序列进行分离，提取子波振幅；
基于所述子波振幅，通过最小二乘法模拟地震子波振幅谱：
Figure 58785DEST_PATH_IMAGE014
其中，
Figure 256548DEST_PATH_IMAGE015
表示大于等于0的常数，
Figure 319051DEST_PATH_IMAGE016
表示地震子波频谱，由
Figure 551449DEST_PATH_IMAGE017
经傅氏变换而来，
Figure 791938DEST_PATH_IMAGE018
和
Figure 476997DEST_PATH_IMAGE019
为待求的关于频率
Figure 828344DEST_PATH_IMAGE020
的多项式，用于拟合地震子波振幅谱；
基于所述模拟地震子波振幅谱，获得子波最大相位分量和最小相位分量；
设子波
Figure 915249DEST_PATH_IMAGE021
的最大相位分量为
Figure 575906DEST_PATH_IMAGE022
、最小相位分量为
Figure 748261DEST_PATH_IMAGE023
，则子波
Figure 903299DEST_PATH_IMAGE024
为：
Figure 844710DEST_PATH_IMAGE025
振幅谱的复赛谱中表示为：
Figure 692581DEST_PATH_IMAGE026
其中，振幅谱的复赛谱
Figure 352232DEST_PATH_IMAGE027
在复赛谱的正、负轴上对称显示，
Figure 291719DEST_PATH_IMAGE028
为地震子波最小相位分量
Figure 87637DEST_PATH_IMAGE029
所对应的最大相位函数的复赛谱，
Figure 106409DEST_PATH_IMAGE030
为地震子波最大相位分量
Figure 987777DEST_PATH_IMAGE031
所对应的最小相位函数的复赛谱；
基于所述振幅谱的复赛谱确定一组具有相同振幅谱的混合相位子波集合，调整俞氏子波参数，在保证地震子波主频完整的前提下，保持低频、扩展高频提升主频，提升有效带宽至预设的带宽阈值，获得扩频模拟子波。
3.根据权利要求1所述的基于多频带地震资料的碳封存场址优选方法，其特征在于，所述地震资料多尺度Bimf分量，其获得方法为：
设置初始的叠前道集矩阵：
Figure 750197DEST_PATH_IMAGE032
其中，
Figure 666200DEST_PATH_IMAGE033
为矩阵行坐标，
Figure 573982DEST_PATH_IMAGE034
为矩阵列坐标，
Figure 942647DEST_PATH_IMAGE035
为叠后三维地震资料数据的值，
Figure 243178DEST_PATH_IMAGE036
为矩阵元素值；
设置变量
Figure 13688DEST_PATH_IMAGE037
为分解的Bimf层数，
Figure 108683DEST_PATH_IMAGE037
的初始值为1，定义初始化矩阵变量
Figure 479490DEST_PATH_IMAGE038
，获得变量矩阵
Figure 318133DEST_PATH_IMAGE039
为：
Figure 943150DEST_PATH_IMAGE040
计算关于
Figure 474625DEST_PATH_IMAGE041
的局部最大值矩阵和局部最小值矩阵：
在变量矩阵
Figure 552302DEST_PATH_IMAGE042
中以预设尺寸的观测矩阵进行观测，提取观测矩阵中的元素值，并获取观测矩阵的最大值；
滑动所述观测矩阵，直至观测矩阵中心遍历所有矩阵变量h，将测得的观测矩阵的最大值赋值于对应观测矩阵的中心元素位置，获得局部最大值矩阵
Figure 460216DEST_PATH_IMAGE043
；
通过设置观测矩阵获取观测矩阵中的最小值，进而获得局部最小值矩阵
Figure 454585DEST_PATH_IMAGE044
；
通过样条插值的方法获取最大值包络面和最小值包络面；
插值公式为：
Figure 891383DEST_PATH_IMAGE045
其中，A、B和C表示插值过度矩阵：
Figure 721936DEST_PATH_IMAGE046
Figure 167960DEST_PATH_IMAGE047
Figure 501990DEST_PATH_IMAGE048
其中，
Figure 375268DEST_PATH_IMAGE049
和
Figure 679735DEST_PATH_IMAGE050
为预设的大于0且小于1的值，
Figure 663871DEST_PATH_IMAGE051
为扩频模拟子波；
求得
Figure 383566DEST_PATH_IMAGE052
为最大值包络面，将过度公式B中的局部最大值矩阵
Figure 427745DEST_PATH_IMAGE053
替换为局部最小值矩阵
Figure 967311DEST_PATH_IMAGE054
，获得
Figure 755138DEST_PATH_IMAGE055
为最小值包络面；
计算最大值包络面和最小值包络面的均值，获得中间过度量：
Figure 313027DEST_PATH_IMAGE056
将变量矩阵
Figure 528108DEST_PATH_IMAGE057
减去中间过度量，获得地震资料对应的第l层Bimf分量：
Figure 289391DEST_PATH_IMAGE058
获取余量
Figure 615330DEST_PATH_IMAGE059
作为计算下一层Bimf分量的输入量：
Figure 44037DEST_PATH_IMAGE060
Figure 413707DEST_PATH_IMAGE061
表示第l层的余量，
Figure 193445DEST_PATH_IMAGE062
表示第
Figure 323075DEST_PATH_IMAGE063
层的余量，第
Figure 340709DEST_PATH_IMAGE063
层的余量作为计算第l层Bimf分量的输入量；
迭代计算Bimf分量，直至满足终止条件：
Figure 163172DEST_PATH_IMAGE064
获得地震资料多尺度Bimf分量；其中，r为预设的终止阈值。
4.根据权利要求1所述的基于多频带地震资料的碳封存场址优选方法，其特征在于，所述时深转换关系，其获取方式为：
基于每个已知井位的测井资料中的声波时差曲线和密度曲线做乘积运算获取波阻抗曲线，进而计算反射系数曲线；
以目标层段地震主频为依据构建雷克子波，将雷克子波与反射系数曲线褶积计算后，得到合成地震记录；
将每个钻井井位井眼处的标志层深度数据与标志层的等时三维展布对应，计算合成地震记录与井旁地震道扩频模拟子波的相关性，当波形相关性高于预设的第一相关阈值时，初步井震标定完成，获得测井深度与地震反射波双程旅行时之间的初步时深转化关系：
Figure 148314DEST_PATH_IMAGE065
其中，
Figure 81635DEST_PATH_IMAGE066
表示声波测井标志层深度对应的地震资料的双程旅行时；
Figure 688197DEST_PATH_IMAGE067
为声波时差；
Figure 415981DEST_PATH_IMAGE068
为测井曲线数据采样间隔；
Figure 904732DEST_PATH_IMAGE069
为地震波双程旅行时；
将所述地震资料多尺度Bimf分量以尺度由大到小的顺序，逐步添加Bimf分量获得更新后的有效地震资料；
每添加一个Bimf分量获得更新后的有效地震资料计算一次更新后的合成地震记录与有效地震资料的第二相关性；
随着有效地震资料逐步添加Bimf分量，第二相关性起初呈现上升趋势，当第二相关性出现下降趋势时，取第二相关性峰值时的更新后的合成地震记录与更新后的地震资料，计算第二时深转化关系；
将所述第二时深转化关系作为最终测得的时深转化关系。
5.根据权利要求1所述的基于多频带地震资料的碳封存场址优选方法，其特征在于，所述波阻抗反演数据体，其获取方式包括：
基于所述扩频模拟子波，计算待判别地震道波形与已知井的合成地震记录的波形相关性，根据波形相关性最高的井对应的波阻抗曲线建立初始模型；
利用白噪声满足高斯分布的规律，将测井资料中的波阻抗数据表示为：
Figure 380360DEST_PATH_IMAGE070
其中，
Figure 638166DEST_PATH_IMAGE071
表示测井波阻抗曲线，
Figure 802431DEST_PATH_IMAGE072
表示待求解的地下地层实际波阻抗值，
Figure 512899DEST_PATH_IMAGE073
表示随机噪声；
根据中心极限定理，
Figure 522443DEST_PATH_IMAGE074
也满足高斯分布，确定初始目标函数为：
Figure 900335DEST_PATH_IMAGE075
其中，
Figure 953610DEST_PATH_IMAGE076
表示与后验信息有关的函数，
Figure 416952DEST_PATH_IMAGE077
表示基于最优样本数选取样本井对样本井的波阻抗曲线，进行匹配滤波后，求得后验概率统计分布密度，进而计算得到的波阻抗期望值，
Figure 964608DEST_PATH_IMAGE078
表示白噪声的协方差；
基于所述初始目标函数，通过最大后验估计，在目标函数中引入先验信息，获得稳定的目标函数为：
Figure 197007DEST_PATH_IMAGE079
其中，
Figure 437495DEST_PATH_IMAGE080
表示待模拟的特征参数，
Figure 371822DEST_PATH_IMAGE081
表示与地质和测井资料先验信息有关的函数，
Figure 988748DEST_PATH_IMAGE082
表示用于协调
Figure 75653DEST_PATH_IMAGE076
和
Figure 487043DEST_PATH_IMAGE081
之间的相互影响的平滑参数；
以所述稳定的目标函数作为初始模型的输入，通过马尔科夫链蒙特卡罗方法MCMC和Metropolis-Hastings抽样准则对后验概率分布抽样，不断优化初始模型的参数，选取目标函数取最大值时的解作为随机实现，取多次随机实现的均值作为期望值输出，将所述期望值输出作为波阻抗反演数据体。
6.根据权利要求1所述的基于多频带地震资料的碳封存场址优选方法，其特征在于，所述划定沉积稳定区，具体为：
设三维叠后地震资料中的各采样点数据为
Figure 393819DEST_PATH_IMAGE083
，其中
Figure 548856DEST_PATH_IMAGE084
表示地震测网线号，
Figure 5114DEST_PATH_IMAGE085
代表地震测网道号，
Figure 321826DEST_PATH_IMAGE086
代表1ms采样的地震记录采样点序号；
计算预设采样区域的采样点数据均方差
Figure 981478DEST_PATH_IMAGE087
：
Figure 940207DEST_PATH_IMAGE088
将采样区域进行纵向和横向平移，遍历计算所有采样区域的数据均方差，获得三维方差属性体；
对所述三维方差属性体进行切片，获取方差属性数据在平面上的分布特征，将方差属性值低于总体能量均值的区域作为沉积稳定区。
7.根据权利要求1所述的基于多频带地震资料的碳封存场址优选方法，其特征在于，所述稳定地层波阻抗数据体，其具体获得方式为：
基于所述沉积稳定区，划定预设尺寸的网格，将每个网格节点视为虚拟井井位；
设未知地下波阻抗模型参数m为
Figure 736124DEST_PATH_IMAGE089
维空间向量
Figure 754896DEST_PATH_IMAGE090
；
通过一次地震采集过程获得了观测数据d为k维数据向量
Figure 888461DEST_PATH_IMAGE091
；
通过非线性函数核G建立未知地下波阻抗模型参数和k维数据向量的联系，即获得正演模型：
Figure 385302DEST_PATH_IMAGE092
其中，
Figure 35726DEST_PATH_IMAGE093
表示与地下波阻抗模型参数m无关的随机噪声，服从高斯分布；
基于正演模型构建反演目标函数：
Figure 225399DEST_PATH_IMAGE094
其中，F为观测数据d与预测数据G(m)之间的均方误差；
对所述反演目标函数进行线性化求解；
对所述正演模型进行泰勒展开并略去二次以上的高阶项，获得预测数据简略表达：
Figure 859643DEST_PATH_IMAGE095
其中，
Figure 143862DEST_PATH_IMAGE096
表示依据先验信息建立的初始模型，A为Jacobian矩阵，A的元素为一阶偏微分
Figure 914372DEST_PATH_IMAGE097
；
设
Figure 274946DEST_PATH_IMAGE098
，
Figure 130907DEST_PATH_IMAGE099
，
Figure 969550DEST_PATH_IMAGE100
表示随机噪声，则正演模型的迭代方程为：
Figure 860146DEST_PATH_IMAGE101
其中，
Figure 375309DEST_PATH_IMAGE102
表示迭代了
Figure 718566DEST_PATH_IMAGE103
次后的正演模型，
Figure 626479DEST_PATH_IMAGE104
表示迭代了t+1次后的正演模型，
Figure 106002DEST_PATH_IMAGE105
为加噪预测数据；
通过一次地震采集过程获得叠前深度域偏移地震剖面
Figure 808379DEST_PATH_IMAGE106
，
假设反演深度域模型参数有先验概率分布
Figure 638932DEST_PATH_IMAGE107
，根据贝叶斯公式有概率分布为：
Figure 803066DEST_PATH_IMAGE108
加入随机噪音的概率分布为：
Figure 402674DEST_PATH_IMAGE109
其中，
Figure 275952DEST_PATH_IMAGE110
为噪音协方差矩阵，
Figure 593801DEST_PATH_IMAGE111
为噪音协方差矩阵求逆运算，det为取对应矩阵的行列式；
n为噪音且服从高斯分布均值为0，
Figure 577938DEST_PATH_IMAGE112
为一等值的对角矩阵且对角线元素为数据中心包含噪音的方差
Figure 32053DEST_PATH_IMAGE113
，则加入随机噪音的概率分布变形为：
Figure 322570DEST_PATH_IMAGE114
对于反演深度域模型参数先验概率分布
Figure 127715DEST_PATH_IMAGE115
，设置初始模型
Figure 649963DEST_PATH_IMAGE116
，且有
Figure 224164DEST_PATH_IMAGE117
，等价的深度域模型概率分布为：
Figure 439245DEST_PATH_IMAGE118
；
设等价的深度域模型概率分布服从高斯分布，有加入随机噪声的深度域模型概率分布为：
Figure 466106DEST_PATH_IMAGE119
假设深度域模型参数均值与初始模型
Figure 41313DEST_PATH_IMAGE120
等值，且互不影响，有简化的深度域模型概率分布：
Figure 470020DEST_PATH_IMAGE121
其中，
Figure 856002DEST_PATH_IMAGE122
表示真实模型相对于所给定初始模型的方差值；
加入随机噪音的概率分布变形中的分母
Figure 635740DEST_PATH_IMAGE123
依赖于数据采集和处理，设为常数，通过后验概率最大化等价于分子最大化，使后验概率分布式最大：
Figure 499790DEST_PATH_IMAGE124
使后验概率分布式最大等价式为：
Figure 517425DEST_PATH_IMAGE125
令关于模型参数m的偏导数等于零，且
Figure 589155DEST_PATH_IMAGE126
，
Figure 590609DEST_PATH_IMAGE127
，获得矩阵方程：
Figure 258351DEST_PATH_IMAGE128
Figure 396071DEST_PATH_IMAGE129
其中，
Figure 389435DEST_PATH_IMAGE130
，
Figure 347027DEST_PATH_IMAGE131
为单位矩阵，
Figure 333306DEST_PATH_IMAGE132
为波阻抗迭代矩阵，整理得到随机逆反演的基本公式：
Figure 325533DEST_PATH_IMAGE133
Figure 224219DEST_PATH_IMAGE134
随机逆反演的基本公式的迭代公式为：
Figure DEST_PATH_IMAGE135
；
其中，
Figure 934686DEST_PATH_IMAGE136
为由参数模型所形成的深度域合成地震记录；
通过不断迭代更新未知地下波阻抗模型的参数，直至
Figure DEST_PATH_IMAGE137
出现正负值震荡，此时的
Figure 930848DEST_PATH_IMAGE138
为低频波阻抗反演数据；
基于虚拟井井位平面坐标与地震数据平面坐标的相关关系，进而确定虚拟井与所述低频波阻抗反演数据的一一对应关系，进而将虚拟井井旁低频波阻抗反演数据对虚拟井赋值，获得虚拟井井位的波阻抗数据；
基于波阻抗反演数据体提取虚拟井井位波阻抗数据，根据地层格架限 定的时窗范围针对所有虚拟井波阻抗数据的插值计算，设定计算区域为地质格架内工区的全部范围，得到稳定地层波阻抗数据体。
8.一种基于多频带地震资料的碳封存场址优选系统，其特征在于，所述系统包括：
资料获取模块，配置为获取已知井位的原始地震资料和测井资料，进而获取井旁地质解释结果和三维叠后地震资料；基于所述三维叠后地震资料，获取目标层位标志层深度数据和标志层的等时三维展布；
预处理模块，配置为基于所述测井资料，进行剔除异常值和标准化的预处理获得标准化测井数据；
扩频模拟模块，配置为基于所述三维叠后地震资料，进行地震子波扩频模拟，获得扩频模拟子波；
等时地层格架模型获取模块，配置为基于所述扩频模拟子波，构建目标层位的等时地层格架模型，并计算断溶体几何结构和空间分布；
Bimf分量获取模块，配置为基于所述三维叠后地震资料和扩频模拟子波，进行多尺度分解获得地震资料多尺度Bimf分量；
时深转换关系获取模块，配置为基于所述扩频模拟子波和标准化测井数据，进行井震标定并忽略地震资料多尺度Bimf分量中的小尺度Bimf分量，获得时深转换关系；
波阻抗反演数据体获取模块，配置为基于所述测井资料、提频后的地震资料和时深转换关系，在所述等时地层格架模型提供的地层趋势信息基础上，进行波形指示地震波阻抗反演获得波阻抗反演数据体；
未定地层波阻抗数据体获取模块，配置为基于所述三维叠后地震资料，计算三维方差属性体，并划定沉积稳定区，基于所述沉积稳定区选择虚拟井井位，并通过虚拟井井间波阻抗插值，获取稳定地层波阻抗数据体；
波阻抗异常数据体获取模块，配置为计算所述稳定地层波阻抗数据体与波阻抗反演数据体之差，获得波阻抗异常数据体；
断溶体波阻抗数据体获取模块，配置为通过去除三维方差属性体中低于平均值的区域，保留断溶体的空间几何轮廓内的波阻抗异常数据，获得包含断溶体几何结构与内部波阻抗特征的断溶体波阻抗数据体；
模型解释模块，配置为比对所述井旁地质解释结果和断溶体波阻抗数据体，划定洞穴储层特征值区间、裂缝破碎带特征值区间和围岩特征值区间，获得断溶体储层解释模型；
根据所述断溶体储层解释模型，计算二氧化碳注入后的运移储集空间，寻找具有预设尺寸阈值的圈闭构造，使得注入后的二氧化碳在断裂带旁侧致密碳酸盐岩封堵以及顶部泥灰岩盖层密闭条件下以高密度游离相存储的运移储集空间，选定为碳封存箱位置。
9.一种电子设备，其特征在于，包括：至少一个处理器；以及与至少一个所述处理器通信连接的存储器；其中，所述存储器存储有可被所述处理器执行的指令，所述指令用于被所述处理器执行以实现权利要求1-7任一项所述的基于多频带地震资料的碳封存场址优选方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于被所述计算机执行以实现权利要求1-7任一项所述的基于多频带地震资料的碳封存场址优选方法。
CN202210919676.5A 2022-08-02 2022-08-02 基于多频带地震资料的碳封存场址优选方法、系统和设备 Active CN114966856B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210919676.5A CN114966856B (zh) 2022-08-02 2022-08-02 基于多频带地震资料的碳封存场址优选方法、系统和设备
US17/902,842 US11852771B1 (en) 2022-08-02 2022-09-04 Method and system for optimally selecting carbon storage site based on multi-frequency band seismic data and equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210919676.5A CN114966856B (zh) 2022-08-02 2022-08-02 基于多频带地震资料的碳封存场址优选方法、系统和设备

Publications (2)

Publication Number Publication Date
CN114966856A CN114966856A (zh) 2022-08-30
CN114966856B true CN114966856B (zh) 2022-12-02

Family

ID=82969485

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210919676.5A Active CN114966856B (zh) 2022-08-02 2022-08-02 基于多频带地震资料的碳封存场址优选方法、系统和设备

Country Status (2)

Country Link
US (1) US11852771B1 (zh)
CN (1) CN114966856B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117434604A (zh) * 2023-02-20 2024-01-23 中国石油化工股份有限公司 地震数据处理方法、装置、存储介质及电子设备
CN117094233B (zh) * 2023-10-19 2024-01-30 中国科学院地质与地球物理研究所 深地油气精准导航三维速度场高精度建模方法与系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107329171A (zh) * 2017-06-07 2017-11-07 中国石油天然气股份有限公司 深度域储层地震反演方法及装置

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102466816B (zh) * 2010-11-04 2014-04-02 中国石油天然气集团公司 一种叠前地震数据地层弹性常数参数反演的方法
EP2917859A4 (en) * 2012-11-07 2016-08-24 Exxonmobil Upstream Res Co METHOD OF KNOWING KNOWLEDGE AND PATTERN RECOGNITION FOR DETECTION OF HYDROCARBON COLLECTIONS
CN104391319B (zh) 2014-12-09 2018-02-02 中国石油天然气集团公司 一种地震资料采集系统的确定方法及装置
CN104536042B (zh) 2014-12-24 2017-01-25 中国石油天然气股份有限公司 一种二维叠后地震资料振幅补偿方法及装置
CN106154323B (zh) * 2015-04-01 2018-08-17 中国石油化工股份有限公司 基于地震拓频处理的相控随机反演薄储层预测方法
CN107390271B (zh) 2017-07-06 2019-06-11 中国石油天然气集团公司 地震资料的振幅补偿方法和装置
CN111208561B (zh) * 2020-01-07 2020-09-01 自然资源部第一海洋研究所 基于时变子波与曲波变换约束的地震声波阻抗反演方法
CN111337980B (zh) * 2020-04-16 2021-03-16 中国矿业大学(北京) 基于时移全波形反演的二氧化碳封存监测方法和系统
CN111399048B (zh) * 2020-04-29 2020-12-08 四川杰瑞泰克科技有限公司 一种对断溶体计算相关属性及数据加权重构的方法
CN111551992B (zh) * 2020-05-12 2021-02-26 中国科学院地质与地球物理研究所 岩石储层构造表征方法、装置、计算机可读存储介质及电子设备
CN112698392A (zh) * 2020-11-20 2021-04-23 中国石油天然气股份有限公司 一种多层系储集体立体刻画与油气空间分布及定容方法
CN112698382B (zh) * 2020-12-04 2023-09-26 中国石油天然气股份有限公司 小尺度断控岩溶型储层地震预测方法及装置
CN113759424B (zh) * 2021-09-13 2022-03-08 中国科学院地质与地球物理研究所 基于频谱分解和机器学习的岩溶储层充填分析方法和系统

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107329171A (zh) * 2017-06-07 2017-11-07 中国石油天然气股份有限公司 深度域储层地震反演方法及装置

Also Published As

Publication number Publication date
CN114966856A (zh) 2022-08-30
US11852771B1 (en) 2023-12-26

Similar Documents

Publication Publication Date Title
CN113759425B (zh) 井震联合评价深层古岩溶储层充填特征的方法与系统
CN114966856B (zh) 基于多频带地震资料的碳封存场址优选方法、系统和设备
CN113759424B (zh) 基于频谱分解和机器学习的岩溶储层充填分析方法和系统
CN106842307B (zh) 一种基于正演约束下波形分类再检索的储层精细预测方法
US11054537B2 (en) Feature index-based feature detection
CA2818790C (en) Seismic trace attribute
US9952341B2 (en) Systems and methods for aligning a monitor seismic survey with a baseline seismic survey
CN114994759B (zh) 基于gan网络的碳封存箱智能识别方法和系统
CN114994758B (zh) 碳酸盐岩断控储层的波阻抗提取与结构表征方法和系统
CN103257363A (zh) 一种探测地下裂缝性储层中裂缝倾角的方法
Imam et al. Mapping the geological structures in the Ras El Ush field (Gulf of Suez, Egypt), based on seismic interpretation and 3D modeling techniques
EA030770B1 (ru) Система и способ адаптивной сейсмической оптики
CN114137616A (zh) 一种利用量子力学原理进行储层含气性检测的方法
CN106842291B (zh) 一种基于叠前地震射线阻抗反演的不整合圈闭储层岩性预测方法
Eikrem et al. Bayesian estimation of reservoir properties—effects of uncertainty quantification of 4D seismic data
EP3359983B1 (en) Seismic polynomial filter
CN105093330A (zh) 多地震道积累振幅差分解谱寻找井旁缝洞储集体的方法
CN112147700A (zh) 速度异常区的低频模型构建方法及系统
Kumar et al. A model-based approach for integration analysis of well log and seismic data for reservoir characterization
CN110568488A (zh) 一种基于非线性混沌算法的生物礁储层识别方法
CN113109875B (zh) 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法
Tian et al. Geology-geophysics-data mining solution enhances the identification of deep fault-karst carbonate reservoirs
Zeng et al. Intelligent Detection of Small Faults Using a Support Vector Machine. Energies 2021, 14, 6242
CN116953794A (zh) 深度域地震反演的低波数初始模型建立方法、装置及介质
CN115616667A (zh) 一种裂缝预测方法、系统、电子设备及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant