发明名称 高分辨率转换波裂缝预测方法
摘要 本发明提供了一种高分辨率转换波裂缝预测方法,所述方法包括:采集转换波地震数据,进行预处理和叠前偏移后得到转换波偏移剖面;在偏移剖面上建立稀疏约束和波阻抗约束的多薄层目标函数;对多薄层目标函数利用频域映射与L1范数联合优化算法求取高精度反射系数;构建高频子波同时保持该高频子波具有确定的频率宽度,用高频子波与高精度反射系数得到高分辨率转换波偏移剖面;利用改进的第三代特征值相干体计算法,对高分辨率转换波偏移剖面进行相干数据体切片提取,从而实现裂缝预测。本发明能够使得提取的裂缝预测剖面分辨率更高、裂缝走向更加清晰、细节刻画更加明显。
申请公布号 CN106443775A 申请公布日期 2017.02.22
申请号 CN201610352020.4 申请日期 2016.05.25
申请人 中国石油集团川庆钻探工程有限公司地球物理勘探公司 发明人 张华;何光明;周强;赵尧;金德刚;罗红明;谢小
分类号 G01V1/30(2006.01)I 主分类号 G01V1/30(2006.01)I
代理机构 成都中玺知识产权代理有限公司 51233 代理人 谭昌驰;薛义丹
主权项 一种高分辨率转换波裂缝预测方法,所述方法包括:采集转换波地震数据,进行预处理和叠前偏移后得到转换波偏移剖面,在转换波偏移剖面上取n道地震记录x(x,t),x=1,2,…,n;在偏移剖面上建立稀疏约束和波阻抗约束的多薄层目标函数;对多薄层目标函数利用频域映射与L1范数联合优化算法求取高精度反射系数;构建高频子波同时保持该高频子波具有确定的频率宽度,用高频子波与高精度反射系数得到高分辨率转换波偏移剖面;利用改进的第三代特征值相干体计算法,对高分辨率转换波偏移剖面进行相干数据体切片提取,从而实现裂缝预测;其中,所述目标函数为下式(1):<maths num="0001"><math><![CDATA[<mrow><mi>O</mi><mrow><mo>(</mo><msub><mi>r</mi><mi>e</mi></msub><mo>,</mo><msub><mi>r</mi><mi>o</mi></msub><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mo>|</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>a</mi><mi>e</mi></msub><mrow><mo>(</mo><msub><mi>b</mi><mi>e</mi></msub><mo>-</mo><msub><mi>A</mi><mi>e</mi></msub><mo>&times;</mo><msub><mi>r</mi><mi>e</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>a</mi><mi>o</mi></msub><mrow><mo>(</mo><msub><mi>b</mi><mi>o</mi></msub><mo>-</mo><msub><mi>A</mi><mi>o</mi></msub><mo>&times;</mo><msub><mi>r</mi><mi>o</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>|</mo><msubsup><mo>|</mo><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mi>&lambda;</mi><mo>|</mo><mo>|</mo><msub><mi>r</mi><mi>e</mi></msub><mo>+</mo><msub><mi>r</mi><mi>o</mi></msub><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>+</mo><mo>&part;</mo><mo>|</mo><mo>|</mo><mfrac><mrow><msup><mrow><mo>(</mo><mi>C</mi><mi>r</mi><mo>-</mo><msub><mi>&xi;</mi><mi>p</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mrow><mo>(</mo><mi>C</mi><mi>r</mi><mo>-</mo><msub><mi>&xi;</mi><mi>p</mi></msub><mo>)</mo></mrow></mrow><mn>2</mn></mfrac><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub><mo>,</mo></mrow>]]></math><img file="FDA0000999092840000011.GIF" wi="1541" he="183" /></maths>其中,r表示反射系数,a<sub>e</sub>代表反射系数的偶分量比例,a<sub>o</sub>代表反射系数的奇分量比例,r<sub>e</sub>为反射系数的偶分量,r<sub>o</sub>为反射系数的奇分量,A为变换矩阵,b为地震数据和子波相关矩阵,λ为稀疏因子,<img file="FDA0000999092840000013.GIF" wi="32" he="51" />为阻抗因子,C为积分算子矩阵,ξ<sub>p</sub>为纵波波阻抗矩阵。
地址 610213 四川省成都市双流县华阳镇华阳大道一段216号