[发明专利]各向异性介质深度域成像的方法及装置

权利要求书

1.一种各向异性介质深度域成像的方法，其特征在于，包括：

根据工区的各向同性深度域解释层位数据，计算工区每一层位的偏移地层厚度；

根据工区每一层位的偏移地层厚度和井位处的井震误差，对工区每一层位进行插值计算，确定工区每一层位结合地层厚度后的井震误差曲面；

根据工区每一层位结合地层厚度后的井震误差曲面，建模得到工区的各向异性介质参数模型；

对工区的各向异性介质参数模型进行各向异性叠前深度偏移处理，确定工区的各向异性介质深度域成像数据。

2.如权利要求1所述的方法，其特征在于，根据工区的各向同性深度域解释层位数据，计算工区每一层位的偏移地层厚度，包括：

根据工区的各向同性深度域解释层位数据，确定工区每一层位的偏移地层深度；

根据工区每一层位的偏移地层深度，计算工区每一层位的偏移地层厚度。

3.如权利要求2所述的方法，其特征在于，按如下公式根据工区每一层位的偏移地层深度，计算工区每一层位的偏移地层厚度：

其中，为工区第i层位的偏移地层厚度；为工区第i层位的偏移地层深度，在i＝1时，

4.如权利要求1所述的方法，其特征在于，根据工区每一层位的偏移地层厚度和井位处的井震误差，对工区每一层位进行插值计算，确定工区每一层位结合地层厚度后的井震误差曲面，包括：

根据工区每一层位的偏移地层厚度和井位处的井震误差，确定工区每一层位的基于各向异性介质的Thomsen参数的δ参数的比例参数；

对工区每一层位的基于各向异性介质的Thomsen参数的δ参数的比例参数进行插值计算，确定工区每一层位结合地层厚度后的井震误差曲面。

5.如权利要求4所述的方法，其特征在于，还包括：根据工区每一层位结合地层厚度后的井震误差曲面，确定工区每一层位结合地层厚度后的每一井位处的井震误差；

按如下公式确定工区每一层位的基于各向异性介质的Thomsen参数的δ参数的比例参数：

其中，为工区第i层位的基于各向异性介质的Thomsen参数的δ参数的比例参数；Mⁱ(x_w,y_w)为工区第i层位的第w井位处的井震误差；为工区第i-1层位的结合地层厚度后的第w井位处的井震误差，i＝1时，为工区第i层位的偏移地层厚度。

6.如权利要求4所述的方法，其特征在于，按如下公式确定工区每一层位结合地层厚度后的井震误差曲面：

其中，为工区第i层位的结合地层厚度后的井震误差曲面；为工区第i层位的基于各向异性介质的Thomsen参数的δ参数的比例参数；为工区第i-1层位的结合地层厚度后的井震误差曲面，i＝1时，为工区第i层位的偏移地层厚度。

7.如权利要求1所述的方法，其特征在于，根据工区每一层位结合地层厚度后的井震误差曲面，建模得到工区的各向异性介质参数模型，包括：

获取工区的各向同性深度域速度；

从工区的各向同性深度域解释层位数据中，提取地层倾角和方位角，得到工区的倾角和方位角属性数据；

根据工区的倾角和方位角属性数据，每一层位结合地层厚度后的井震误差曲面，各向同性深度域速度和各向同性深度域解释层位数据，建立工区的网格层析反演矩阵；

对建立的网格层析反演矩阵，进行网格层析反演，得到工区的各向异性介质参数模型；其中，各向异性介质参数模型包括如下参数：v_P0参数、δ参数和ε参数，其中v_P0参数为纵波垂向速度，δ参数和ε参数为各向异性介质的Thomsen参数。