CN119439271A 一种相位可控高效的逆时偏移成像方法及装置 （大庆油田有限责任公司）

一种相位可控高效的逆时偏移成像方法及本发明提供了一种相位可控高效的逆时偏后相位整形处理时无法解决子波空变和波场复2S3.对步骤S1中炮集数据中的任一炮地震数据，完成优化时间的震源波场正向延拓和S4.在步骤S3中优化时间的任一时间步，在该炮数据的成像空间中的任一成像点位置S6.根据步骤S5完成工区所有炮数据的逆时偏移处理，并通过对所有单炮逆时成像数所述地表一致性反褶积处理和地表一致性能量补偿处理均采用一个优选的，所述多信息约束网格层析速度反演方法，通过对12)完成第一轮深度域成像道集的剩余延迟谱计算和种子点拾取和，深度方向拾取密14)完成第二轮深度域成像道集的剩余延迟谱计算和种子点拾取和，深度方向拾取密模型与步骤12)中的表征大尺度构造特征的速度模型相加，由此获得表征中大尺度的复杂316)完成第三轮深度域成像道集的剩余延迟谱计算和种子点拾取，深度方向拾取密度17)当所述深度域速度模型对应的深度域成像数据体满足地质需求、速度模型与成像优选的，所述步骤17)满足所述地质需求的因素为构造成像位置准确，地质体成像清所述速度模型与成像结果构造匹配的条件为：深度域速度模型所述步骤S2所述产生震源波场的地震子波包括采用零相所述步骤S2所采用的地震波动方程类型，调21)若要使逆时成像结果与克希霍夫深度成像结果的波组特征保持一致，则需要对产22)若要使逆时成像结果具有近似零相位的波组特征，则产生震源波场的输入地震子优选的，所述采用的地震波动方程类型，包4所述步骤S3炮集数据中的任一炮地震数据，完成优化时间在步骤S2中地震波动方程x0,y0,z0位置加入产生震源波场的地震子波j(x,J,z,s,0),6.根据权利要求1或4所述的一种相位可控高效的述S3炮集数据中的任一炮地震数据，完成优化时间的检波点波场逆时延拓的具体过程如。更精确地，可以根据射线追踪法或基于程函方程的5优选的，所述振幅补偿处理的方法为通过对步骤S所述振幅补偿处理中作为分母的振幅能量数据体,需要对振幅能量数据体的绝对值数采用该振幅能量数据体的最大绝对值的[8098％]之间的值作为修改后振幅能量数采用该振幅能量数据体的最大绝对值的[110％]之间的值作为修改后振幅能量数上述处理可以消除因全区逆时偏移数据体与全区振幅能量数据体相除运算时分母为执行单元，炮集数据中的任一炮地震数据，完成优化时间的6第二获取单元，基于完成的炮集数据中任一单第一处理单元，炮集数据中任一单炮数据的逆时第二处理单元，用于对步骤S5完成工区所有炮7[0004]本发明所要解决的技术问题是克服背景技术中存在成像叠加结果上无法解决深度域地震子波深变和空变的相位校正的问题,而提供一种相位可控高效的逆时偏移成像方[0008]S3.对步骤S1中炮集数据中的任一炮地[0009]S4.在步骤S3中优化时间的任一时[0010]S5.根据步骤S3、步骤S4完成步骤S18[0011]S6.根据步骤S5完成工区所有炮[0013]所述地表一致性反褶积处理和地表一致性能量补偿处理均采用[0015]所述多信息约束网格层析速度反演方法，通过对所述深度域速度模型进行迭代,[0017]12)完成第一轮深度域成像道集的剩余延迟谱计算和种子点拾取和，深度方向拾[0019]14)完成第二轮深度域成像道集的剩余延迟谱计算和种子点拾取和，深度方向拾速度模型与步骤12)中的表征大尺度构造特征的速度模型相加，由此获得表征中大尺度的[0021]16)完成第三轮深度域成像道集的剩余延迟谱计算和种子点拾取，深度方向拾取度模型与步骤14)中的表征中大尺度的复杂构造速度模型相加，获得表征精细构造的速度[0022]17)当所述深度域速度模型对应的深度域成像数据体满足地质需求、速度模型与9[0024]所述成像道集的拉平标准要求为：深度域成像道集的剩余延迟谱上的能量团归[0026]优选的，所述步骤S2所述产生震源波场的地震子波包括采用零相位的雷克子波、[0028]21)若要使逆时成像结果与克希霍夫深度成像结果的波组特征保持一致，则需要[0029]22)若要使逆时成像结果具有近似零相位的波组特征，则产生震源波场的输入地[0030]23)若要使逆时成像结果具有相位的波组特征，可通过对产生震源波场的地震[0040]在步骤S2中地震波动方程x0,y0,z0位置加入产生震源波场的地震子波空间中的任一成像点位置应用成像条件，并对应用成像条件后的地震波场在0_T时间内进[0053]31)、检波点波场从步骤S1中炮集数据记录总道长的最大时刻开始进行逆时延拓[0065]所述振幅补偿处理中作为分母的振幅能量数据体,需要对振幅能量数据体的绝对[0067]采用该振幅能量数据体的最大绝对值的[8098％]之间的值作为修改后振幅能[0068]采用该振幅能量数据体的最大绝对值的[110％]之间的值作为修改后振幅能[0069]上述处理可以消除因全区逆时偏移数据体与全区振幅能量数据体相除运算时分[0092]对于炮集数据，需要经过地表一致性反褶积处理和地表一致性振幅能量补偿处[0096]12)完成第一轮深度域成像道集的剩余延迟谱计算和种子点拾取和，深度方向拾[0098]14)完成第二轮深度域成像道集的剩余延迟谱计算和种子点拾取和，深度方向拾速度模型与步骤12)中的表征大尺度构造特征的速度模型相加，由此获得表征中大尺度的[0100]16)完成第三轮深度域成像道集的剩余延迟谱计算和种子点拾取，深度方向拾取度模型与步骤14)中的表征中大尺度的复杂构造速度模型相加，获得表征精细构造的速度[0101]这里深度域速度模型迭代次数依据对应的深度域成像数据体是否满足地质需求[0107]3)若要使逆时成像结果具有相位[0118]在步骤S2中地震波动方程x0,y0,z0位置加入产生震源波场的地震子波在本地盘进行存储。[0130]1)、检波点波场从步骤S1中炮集数据记录的总道长最大时刻开始进行逆时延拓[0142]所述振幅补偿处理的方法为通过对步骤S6中全工区逆时偏移数据体与振幅能量数据体相除获得；所述振幅补偿处理中作为分母的振幅能量数据体,需要对振幅能量数据[0144]采用该振幅能量数据体的最大绝对值的[8098％]之间的值作为修改后振幅能[0145]采用该振幅能量数据体的最大绝对值的[110％]之间的值作为修改后振幅能[0146]上述处理可以消除因全区逆时偏移数据体与全区振幅能量数据体相除运算时分[0159]保真预处理后的炮集数据已经过地表一致性反褶积处理和地表一致性振幅能量[0160]高精度的速度模型采用多信息约束网格层析速度反演方法，采用了3轮速度从粗[0161]1)利用叠前时间偏移速度进行时深转换构建初始深度域速度模型(如图1所示)，密度为200米(即40倍的深度步长)；并完成立体网格层析速度反演，空间方向反演网格为密度为100米(即20倍的深度步长)；并完成立体网格层析速度反演，空间方向反演网格为倍的深度步长)，本步骤中立体网格层析反演速度模型与步骤4)中的表征中大尺度的复杂[0167]2.为研究震源波场地震子波对逆时偏移结果的波组相位特征的影响及其与常规零相位地震子波和90度相位旋转后的地震子波作为震源地震子波开展逆时成像处理应用，位旋转后的地震子波可以使逆时成像结果与克希霍夫深度成像结[0171]3.在步骤2中地震波动方程x0,y0,z0位置加入产生震源波场的地震子波在该炮数据的成像空间范围内，采用步骤2中的地震波动方程计算公式沿着时间变大的方向进行优化时间的震源波场的正向延拓，即从0时刻开始，一直传播到T=[0172]4.采用如下二阶形式的地震波动方程对步骤1炮集中的任一炮数据进行逆时延[0176]5.当步骤3中检波点波场从6s时刻一直