[发明专利]基于频谱分解和机器学习的岩溶储层充填分析方法和系统

权利要求书

1.一种基于频谱分解和机器学习的岩溶储层充填分析方法，其特征在于，所述方法包括：

步骤S100，获取原始地球物理测井资料：通过测井设备获取各样本井的原始测井数据，包括：通过测量电极测量各样本井的自然电位SP，通过自然伽马井下装置和自然伽马地面仪器测量各样本井的自然伽马GR，通过井径臂获取各样本井的井径CAL；通过传统测井设备获取电阻率曲线数据：深侧向测井RLLD、浅侧向测井RLLS和微侧向测井RLLM，物性表征曲线数据：补偿中子CNL、补偿声波曲线AC和密度曲线DEN；

基于成像测井信息、钻井信息、录井信息和岩心信息获得个别深度段的确定岩性信息和物性信息，进而确定目的层位标志层的深度数据；

步骤S200，获取地震数据，通过地震波激发装置与接收装置获取原始地震波反射信号数据，并根据原始地震波反射信号数据的波形获取目的层位标志层的等时三维展布；

步骤S300，原始地球物理测井数据预处理：基于所述样本井的所有原始测井数据，绘制测井曲线数据，并进行异常值处理和标准化处理，获得标准化测井曲线数据；

步骤S400，地震数据预处理：基于所述地震波反射信号数据，通过混合相位子波反褶积和扩散滤波，获得高精度的三维地震振幅数据体；

步骤S500，井震标定与特征参数选取：基于所述标准化测井曲线数据中的补偿声波曲线AC和密度曲线DEN获取样本井的波阻抗曲线，进而计算反射系数曲线，获取雷克子波的优选频率使其与高精度三维地震振幅数据体的主频保持一致，将所述雷克子波与所述反射系数曲线进行褶积运算，得到合成地震记录，并将所述目的层位标志层的深度数据与所述目的层位标志层的等时三维展布相比对进行井震标定，计算所述合成地震记录与井旁地震道波形的相关度，当相关度大于或等于预设的第一阈值时，判定井震标定结果合格，获得测井曲线数据与地震记录的时深转化关系和对储层敏感的特征参数；

所述对储层敏感的特征参数，其获得方法包括：

通过对井旁不同地质体产生的测井参数的响应绘制柱状统计图，当某一标准化测井曲线数据的数值可以将不同测井解释结论预设的第二阈值以上的数据点区分开时，则选定该标准化测井曲线数据为对储层敏感的特征参数；所述对储层敏感的特征参数至少包括波阻抗IMP，还可包括井径CAL、自然伽马GR、自然电位SP，电阻率曲线数据：深侧向测井RLLD、浅侧向测井RLLS、微侧向测井RLLM，物性表征曲线数据：补偿中子CNL、补偿声波曲线AC、密度曲线DEN中的一个或多个；

步骤S600，构建等时格架模型：基于高精度三维地震振幅数据体所反映的沉积地层规律和所述时深转化关系构建等时格架模型；

步骤S700，井间储层参数模拟：基于各样本井的标准化测井曲线数据，确定能够反映整体地质情况的最优样本数参量，选取最优样本数参量条地震波形相关性最高的样本井的对储层敏感的特征参数数据构建初始模型，不断更正初始模型的参数，输出高精度特征值模拟结果数据体，所述高精度特征值模拟结果数据体为与所述对储层敏感的特征参数一一对应的数据体；

步骤S800，刻画井间溶洞系统边界：基于对储层敏感的特征参数中的样本井的波阻抗曲线IMP，绘制波阻抗曲线直方图，确定区分古岩溶洞穴与围岩波阻抗门槛值；

基于所述古岩溶洞穴与围岩波阻抗门槛值，将高精度特征值模拟结果数据体中超过所述区分古岩溶洞穴与围岩波阻抗门槛值的区域定义为碳酸盐岩基岩，将低于所述区分古岩溶洞穴与围岩波阻抗门槛值的区域定义为古岩溶洞穴储层发育位置，获得古岩溶洞穴储层结构特征和空间分布规律；

步骤S900，刻画洞穴内部充填岩性物性边界：刻画洞穴内部充填岩性物性边界：基于所述对储层敏感的特征参数构建交会版图，根据所述个别深度段的岩性信息和物性信息在交汇版图中框选数据点作为测井解释样本点，获得充填程度与充填物类型的解释结论门槛值；

通过所述解释结论门槛值对所述高精度特征值模拟结果数据体进行交会分析获得洞穴内部充填程度与充填物组合三维空间形态特征；

步骤S1000，古岩溶洞穴结构与充填描述：基于所述古岩溶洞穴储层结构特征、空间分布规律、基于洞穴内部充填程度和充填物组合三维空间形态特征，采用岩性遮挡技术和三维雕刻技术雕刻出古岩溶洞穴空间展布和内部充填的发育特征。