主权项 |
一种高分辨率转换波裂缝预测方法,所述方法包括:采集转换波地震数据,进行预处理和叠前偏移后得到转换波偏移剖面,在转换波偏移剖面上取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>×</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>×</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>λ</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>∂</mo><mo>|</mo><mo>|</mo><mfrac><mrow><msup><mrow><mo>(</mo><mi>C</mi><mi>r</mi><mo>-</mo><msub><mi>ξ</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>ξ</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>为纵波波阻抗矩阵。 |