CN113759425B 井震联合评价深层古岩溶储层充填特征的方法与系统 （中国科学院地质与地球物理研究所）

所号井震联合评价深层古岩溶储层充填特征的本发明属于数据识别、记录载体的处理领大范围内复杂盆地碳酸盐岩洞穴型储层的发育得高精度的三维地震振幅数据体；通过PCA方法过解释结论门槛值进行交会分析获得古岩溶洞穴空间展布和不同填充物类型的洞穴内部充填2基于成像测井信息、钻井信息、录井信息和岩心信息步骤S200，获取地震数据，通过地震波激发装步骤S700，构建等时格架模型：基于高精度三维3基于所述岩性信息和物性信息获取第一组各类地质体离散数据点分布图的解释结论，基于所述第一组PCA数据门槛值，对所述第一组高精度特征值模拟结果数据体进行交步骤S1000，刻画洞穴内部充填岩性物性边界：基洞穴数据点，将所述岩溶洞穴数据点所对应的对储层敏感的特征参数通过PCA方法进行分基于所述第二组PCA数据，通过如S800所述方法生成第二组高精度特征值模拟结果数基于所述岩性信息和物性信息获取第二组各类地质体离散数据点分布图的解释结论，基于所述第二组PCA数据门槛值，对所述第二组高精度特征值模拟结果数据体进行交W为地层水电阻率值；通过探测器获取自然伽马射线，并将所述自然伽马射线所述地面仪器把每分钟形成的电脉冲计数转4基于所述协方差矩阵的特征值和特征向量计算包括p种特j表示第j个PCA数据向量，nj表示第计算方差累计贡献量：5步骤S420，将所述频率域地震记录褶积模型的$(t)=⃞(t)+(t)其中，(t)表示地震波形记录的复赛谱序列，表示地震子波的复赛谱序列，(t)表最大相位分量v(ttheo)所对应的最小相位函数的复赛谱，i(-teheo)为地震子波最小相位分步骤S470，基于所述振幅谱中的复赛谱确定一组具有相同振幅谱的混合相位子波集60表示梯度向量张量积；v和v3对应的三个特征值；步骤S490，基于所述结构张量的特征向量分别计算线状步骤S4100，重复步骤S480-S490的步骤，直至将所述地震信号振幅谱的极大值和所述模拟地震子波振幅谱通过最小二乘法拟合的将所述地震信号振幅谱的极大值除以所拟合出的相应频7步骤S820，根据波形相似性原则选取数量为样本数参量的样本井第一组PCA数据与参考目标井标第一组PCA数据进行相关性分析获得样本数参量-参考目标井的第一组PCA数据数据相关性值连接,获得参考井第一组PCA数据相关性随样本个参考井第一组PCA数据相关性随样本数参量变化的相关性曲线，将所有参考井的第一组PCA数据波形相关性随样本数参量变化的相关性曲线拟合为整体相关性曲线，选取所述整步骤S850，基于所述高精度的三维地震振幅数据体和等步骤S860，基于所述初始模型，选取最优样本数步骤S880，基于所述最大似然函数和先验信息步骤S890，以所述目标函数作为所述初始模参数对应于所述第一组PCA数据一一对应的特步骤S881，利用白噪声满足高斯分布的规律，8步骤S893，由M中的某一点P0出发，通过马尔科夫链进行抽样模拟，产生点序列：n；＝P(h-1)；迭代次数hs，获得后验样本p(3),p(2),p(3)…p(s),进而计算后验分布的各阶矩阵获得期步骤S330，基于所述去除离群点的测井曲线数据，叠合工区所述原始地球物理测井资料获取模块，配置为通过测9基于成像测井信息、钻井信息、录井信息和岩心信息所述地震数据获取模块，配置为通过地震波激发装置所述原始地球物理测井数据预处理模块，配置为基于所述地震数据预处理模块，配置为基于所述地震波反射信所述井震标定与特征参数选取模块，配置为井震标定化测井曲线数据中的补偿声波曲线和密度曲线DEN获取样本井的波阻抗曲线，进而计算反射系数曲线，获取雷克子波的优选频率使其与高精度三维地震振幅数据体的主频保持一所述合成地震记录与井旁地震道波形的相关度，当相关度大于或等于预设的第一阈值时，所述第一组PCA数据获取模块，配置为从所述对储层敏感的特征参数中优选预设条数所述等时格架模型构建模块，配置为基于高精度三维地数据体，所述第一组高精度特征值模拟结果数据体为与所述第一组PCA数据一一对应的数所述刻画井间溶洞系统边界模块，配置为基于样本井的第一组PCA数据，进行交会分基于所述岩性信息和物性信息获取第一组各类地质体离散数据点分布图的解释结论，基于所述第一组PCA数据门槛值，对所述第一组高精度特征值模拟结果数据体进行交所述洞穴内部充填岩性物性边界刻画模块，配置为刻画洞穴内部充填岩性物性边界：所述岩溶洞穴数据点，将所述岩溶洞穴数据点所对应的对储层敏感的特征参数通过PCA方基于所述第二组PCA数据，通过如井间储层参数模拟模块的功能生成第二组高精度特基于所述岩性信息和物性信息获取第二组类地基于所述第二组PCA数据门槛值，对所述第二组高精度特征值模拟结果数据体进行交所述古岩溶洞穴结构与充填描述模块，配置为：基于所雕刻出古岩溶洞穴空间展布和不同填充物类型的[0002]含油气盆地深层(埋深超过4500m)蕴含丰富的油气资源，已经成为石油工业界勘探的重要领域。碳酸盐岩油藏占全球已探明石油储量的52产量的60是深层油气勘探测井曲线数据与地震记录的时深转化关系和数据体，所述第一组高精度特征值模拟结果数据体为与所述第一组PCA数据一一对应的数[0019]基于所述岩性信息和物性信息获取第一组各类地质体离散数据点分布图的解释[0020]基于所述第一组PCA数据门槛值，对所述第一组高精度特征值模拟结果数据体进[0022]基于所述第二组PCA数据，通过如S800所述方法生成第二组高精度特征值模拟结[0024]基于所述岩性信息和物性信息获取第二组各类地质体离散数据点分布图的解释[0025]基于所述第二组PCA数据门槛值，对所述第二组高精度特征值模拟结果数据体进[0045]步骤B300，基于所述归一化样本特征参数数据xS，计算样本特征参数协方差Cov[0061]当k从1开始不断增加时，选取Ck首次大于预设的第三阈值时对应的前k个协方差[0066]步骤S420，将所述频率域地震记录褶积模型的等式两边取对[0069]$(t)=⃞(t)+(t)[0070]其中，5(r)表示地震波形记录的复赛谱序列，表示地震子波的复赛谱序列，(t)表示地层反射系数的复赛谱序列，t表示地震波形记录时间；[0079]27(t)=i(t)+⃞(t)+i(-t)+(-t)大相位分量v(t)所对应的最小相位函数的复赛谱，i(-t)为地震子波最小相位分量u(t)所[0105]将所述地震信号振幅谱的极大值和所述模拟地震子波振幅谱通过最小二乘法拟[0110]在一些优选的实施方式中，所述步骤S800，具体包括，具体包括步骤S810-步骤[0112]步骤S820，根据波形相似性原则选取数量为样本数参量的样本井第一组PCA数据与参考目标井标第一组PCA数据进行相关性分析获得样本数参量-参考目标井的第一组PCA本数参量-参考目标井的第一组PCA数据相关性值，将所有样本数参量-参考目标井的第一得多个参考井第一组PCA数据相关性随样本数参量变化的相关性曲线，将所有参考井的第一组PCA数据波形相关性随样本数参量变化的相关性曲线拟合为整体相关性曲线，选取所方法(MarkovchainMonteCarlomethod，MCMC)和Metropolis-Hastings抽样准则对后[0129]其中，Zsim表示待模拟的特征参数，J2表示与地质和测井资料先验信息有关的函n；＝P(h-1)；预设的迭代次数hs，获得后验样本p(1),p(2),p3)…p(s),进而计算后验分布的各阶矩阵[0146]本发明的另一方面，提出了一种井震联合评价深层古岩溶储层洞穴内部充填岩性物性边界模块和古岩溶洞[0150]所述原始地球物理测井数据预处理模块，配置为基于所述样本井的原始测井数算所述合成地震记录与井旁地震道波形的相关度，当相关度大于或等于预设的第一阈值[0155]所述第一组PCA数据获取模块，配置为从所述对储层敏感的特征参数中优选预设拟结果数据体，所述第一组高精度特征值模拟结果数据体为与所述第一组PCA数据一一对[0159]基于所述岩性信息和物性信息获取第一组各类地质体离散数据点分布图的解释[0160]基于所述第一组PCA数据门槛值，对所述第一组高精度特征值模拟结果数据体进[0162]基于所述第二组PCA数据，通过如井间储层参数模拟模块所述功能生成第二组高[0164]基于所述岩性信息和物性信息获取第二组各类地质体离散数据点分布图的解释[0165]基于所述第二组PCA数据门槛值，对所述第二组高精度特征值模拟结果数据体进技术雕刻出古岩溶洞穴空间展布和不同填充物类型[0174]本发明通过利用高精度的三维地震波形数据之间的关联性在初始等时地质格架模型基础上井震联合进行三维参数模拟，并运用两次PCA方法(PCA)对多条经预处理后的[0175]本发明还通过基于PCA数据所构建的二维交会图很好地区分岩溶洞穴及其内部充用至三维参数模拟资料中，得到PCA交会分析约束下的岩溶洞穴纵横向结构及其内部充填[0176]本发明在贝叶斯理论指导下的马尔科夫链抽样准则与蒙特卡罗估计法建立地球[0178]图1是井震联合评价深层古岩溶储层充填特征的方法与系统实施例的流程示意[0183]图6是本发明实施例中将所有参考井的波形相关性随样本数参量变化的相关性曲[0184]图7是本发明实施例中第一组PCA数据交会分析抓取岩溶洞穴数据点并框选门槛[0186]图9是本发明实施例中第二组PCA数据交会分析划分岩溶洞穴内部充填物并框选洞穴充填物组合有关(Yuetal.,2018)。Xu等(2016)通过物理模型模拟同样可以发现不[0219]补偿中子测井由一个中子发射源与两个接收探测器组成。通过对两个探测器计电流时，可动电阻的某一固定端与滑动端之间的电位差将随着其间电阻值的变化而变化。子波反褶积(mixed-phasewaveletestimationandmaximumposteriorideconvolution)与扩散滤波(enhancingdiffusionfiltering)进行三维地震数据的拓频[0234]混合相位子波反褶积是一种在确保处理后的地震数据具有较高保真度高的前提下，通过拓宽有效频带从而提高地震信号分辨率的数据处理方法，相当于本实施例中的[0239]步骤S420，将所述频率域地震记录褶积模型的等式两边取对[0242]$(t)=⃞(t)+(t)[0243]其中，5(t)表示地震波形记录的复赛谱序列，表示地震子波的复赛谱序列，(t)表示地层反射系数的复赛谱序列，t表示地震波形记录时间；[0250]将所述地震信号振幅谱的极大值和所述模拟地震子波振幅谱通过最小二乘法拟(-tr子波最大相位分量v(ttheo)所对应的最小相位函数的复赛谱，i(-teheo)为地震子波最小相[0260]Fhemers和Hocker于2003年首次将扩散滤波技术应用在地震资料的处理解释中。选取85%;[0301]不同测井曲线反映了地质体不同物理学参数，由于各种测井曲线的探测方法有[0305]步骤B300，基于所述归一化样本特征参数数据xS，计算样本特征参数协方差Cov[0321]当k从1开始不断增加时，选取Ck首次大于预设的第三阈值时对应的前k个协方差的取85％。前k个特征参数对应的PCA数据反映了表征储层各项参数的Ck部分信息，优选[0324]PC1＝0.3859NGR-0.64[0325]PC2＝0.1425NGR+0.67中的第一数据；PC2表示第一次对所有数据点特征参数曲线进行PCA方法的第一组PCA数据数据体，所述第一组高精度特征值模拟结果数据体为与所述第一组PCA数据一一对应的数[0331]步骤S820，根据波形相似性原则选取数量为样本数参量的样本井第一组PCA数据与参考目标井标第一组PCA数据进行相关性分析获得样本数参量-参考目标井的第一组PCA本数参量-参考目标井的第一组PCA数据相关性值，将所有样本数参量-参考目标井的第一组PCA数据相关性值连接,获得参考井第一组PCA数据相关性随样本数参量变化的相关性曲得多个参考井第一组PCA数据相关性随样本数参量变化的相关性曲线，将所有参考井的第一组PCA数据波形相关性随样本数参量变化的相关性曲线拟合为整体相关性曲线，选取所本数参量；将所有参考井的第一组PCA数据波形相关性随样本数参量变化的相关性曲线拟[0349]其中，Zsim表示待模拟的特征参数，J2表示与地质和测井资料先验信息有关的函方法MCMC和Metropolis-Hastings抽样准则对后验概率分布抽样，不断优化初始模型的参[0352]基于所述岩性信息和物性信息获取第一组类地质体离散数据点分布图的解释结[0354]基于所述第一组PCA数据门槛值，对所述第一组高精度特征值模拟结果数据体进[0356]基于所述第二组PCA数据，通过如S700所述方法生成第二组高精度特征值模拟结[0358]基于所述岩性信息和物性信息获取第二组类地质体离散数据点分布图的解释结[0359]基于所述第二组PCA数据门槛值，对所述第二组高精度特征值模拟结果数据体进所对应的五条特征参数曲线再次进行PCA方法，获得对充填物响应敏感的第二组PCA数据[0362]F1＝0.5574NGR-0.48[0363]F2＝-0.4732NGR-0.46[0364]其中，F1表示第二此对洞穴表征数据点特征参数曲线进行PCA方法的第二组PCAIMP为洞穴表征数据点的波阻抗值[0365]通过所述岩溶洞穴内部沉积充填与垮塌充填划分的解释结论门槛值对所述第二组PCA数据模拟结果数据体进行交会分析获得岩溶洞穴内部沉积充填与垮塌充填的三维空[0370]所述原始地球物理测井数据预处理模块，配置为基于所述样本井的原始测井数算所述合成地震记录与井旁地震道波形的相关度，当相关度大于或等于预设的第一阈值[0376]所述第一组PCA数据获