一种基于基追踪横向多道约束的各向异性参数反演方法

摘要

本发明公开了一种基于基追踪横向多道约束的各向异性参数反演方法，包括：处理地震数据，获取方位角度叠加道集；构建岩石物理模型，获取裂缝参数测井数据；对所述测井数据进行去相关处理；构建反演目标函数；基于所述反演目标函数，获得待反演参数反射系数序列。本发明的基于基追踪横向多道约束的各向异性参数反演方法，通过纳入方位角获得方位角度叠加道集，方位角度叠加道集的信道比高，且包括多重地震信息，进而较为准确地提供储层中的裂缝信息，通过对测井数据进行去相关处理，提高了反演结果的准确度，通过在目标函数中加入相关约束项，构建反演目标函数，提高了反演结果稳定性和横向连续性，以此降低未知参数估测的不确定性。

权利要求

1.一种基于基追踪横向多道约束的各向异性参数反演方法，其特征在于，包括：处理地震数据，获取方位角度叠加道集；

构建裂缝岩石物理模型，获取裂缝参数测井数据；

对所述测井数据进行去相关处理；

构建反演目标函数；

基于所述反演目标函数，获得待反演参数反射系数序列。

2.根据权利要求1所述的各向异性参数反演方法，其特征在于，所述对所述测井数据进行去相关处理包括：生成所述测井数据的协方差矩阵；

对所述协方差矩阵进行奇异值分解，获取去相关矩阵；

基于所述去相关矩阵，对公式d＝Gr进行变换，获得变换后的去相关公式，其中，d为所述方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数。

3.根据权利要求2所述的各向异性参数反演方法，其特征在于，所述变换后的去相关公式为：d＝G×r

其中，G＝G×V，r＝V-1×r，d为所述方位角度叠加道集的地震数据，G为去相关的子波矩阵，G为子波矩阵，V为去相关矩阵，r为去相关的待反演参数，r为待反演参数。

4.根据权利要求1所述的各向异性参数反演方法，其特征在于，所述构建反演目标函数包括：获得所述方位角度叠加道集的似然函数；

基于所述似然函数，获得目标函数；

在所述目标函数中加入相关约束项，构建所述反演目标函数。

5.根据权利要求4所述的各向异性参数反演方法，其特征在于，所述基于所述似然函数，获得目标函数包括：结合所述似然函数与先验分布，获得后验概率分布；

对所述后验概率分布进行最大化，获得所述目标函数。

6.根据权利要求4所述的各向异性参数反演方法，其特征在于，所述相关约束项包括：基于井数据平滑阻抗模型约束、L1范数约束和横向多道约束。

7.根据权利要求4所述的各向异性参数反演方法，其特征在于，所述反演目标函数为：其中，d为方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数，Jsmooth(r)为基于井数据平滑阻抗模型约束，r'为前一道待反演参数反演结果，λ、β为权系数。

8.根据权利要求1所述的各向异性参数反演方法，其特征在于，通过基追踪算法获得所述待反演参数反射系数序列。

9.根据权利要求8所述的各向异性参数反演方法，其特征在于，通过对所述反演目标函数求取最小值获得所述待反演参数反射系数序列。

10.根据权利要求1所述的各向异性参数反演方法，其特征在于，还包括：根据所述待反演参数反射系数序列，获得最终待反演参数，其中，所述最终待反演参数包括纵波阻抗、横波阻抗、裂缝岩石物理参数。

说明书

一种基于基追踪横向多道约束的各向异性参数反演方法

技术领域

[0001] 本发明属于地球物理反演领域，具体涉及一种基于基追踪横向多道约束的各向异性参数反演方法。

背景技术

[0002] 随着勘探目标转向复杂型油气藏，非常规油气勘探开发已成为研究重点。在页岩气的勘探中，有效识别裂缝发育程度是至关重要的，而裂缝型储层具有各向异性的特点，因此，准确获得能够反映裂缝发育程度的裂缝各向异性参数对地下裂缝预测具有重要意义。

[0003] 随着宽方位地震数据采集和处理技术的快速发展，利用宽方位叠前地震数据预测地下裂缝已成为主要研究方法。目前基于方位叠前地震数据的各向异性参数反演方法主要利用地震波在裂缝型储层传播过程中反射系数随入射角和方位角的变化情况来进行研究，其中理论推导主要来源于Ruger在1998年提出的基于等效横向各向同性介质的反射系数近似公式，随后有相关学者对该方程进行了一系列简化，但在各向异性参数反演中仍存在可靠性和稳定性的问题。另外裂缝型储层的各向异性特征对地震波振幅产生的影响，使得在反演过程中出现横向不连续的情况。对此特别需要一种合理的各向异性参数反演方法，能提高反演的可靠性、稳定性和横向连续性，从而有助于开展利用地震信息测量地下各向异性强弱的方法研究，进而较为准确地提供储层中的裂缝信息，指导裂缝型储层的勘探和开发。

发明内容

[0004] 本发明的目的是提出一种提高反演准确性、可靠性、稳定性和横向连续性的基于基追踪横向多道约束的各向异性参数反演方法。

[0005] 为了实现上述目的，本发明提供一种基于基追踪横向多道约束的各向异性参数反演方法，包括：处理地震数据，获取方位角度叠加道集；构建裂缝岩石物理模型，获取裂缝参数测井数据；对所述测井数据进行去相关处理；构建反演目标函数；基于所述反演目标函数，获得待反演参数反射系数序列。

[0006] 优选的，所述对所述测井数据进行去相关处理包括：生成所述测井数据的协方差矩阵；对所述协方差矩阵进行奇异值分解，获得去相关矩阵；基于所述去相关矩阵，对公式d＝Gr进行变换，获得变换后的去相关公式，其中，d为所述方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数。

[0007] 优选的，所述变换后的去相关公式为：d＝G×r，其中，G＝G×V，r＝V-1×r，d为所述方位角度叠加道集的地震数据，G为去相关的子波矩阵，G为子波矩阵，V为去相关矩阵，r为去相关的待反演参数，r为待反演参数。

[0008] 优选的，所述构建反演目标函数包括：获得所述方位角度叠加道集的似然函数；基于所述似然函数，获得目标函数；在所述目标函数中加入相关约束项，构建所述反演目标函数。

[0009] 优选的，所述基于所述似然函数，获得目标函数包括：结合所述似然函数与先验分布，获得后验概率分布；对所述后验概率分布进行最大化，获得所述目标函数。

[0010] 优选的，所述相关约束项包括：基于井数据平滑阻抗模型约束、L1范数约束和横向多道约束。

[0011] 优选的，所述反演目标函数为：

[0012]

[0013] 其中，d为所述方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数，Jsmooth(r)为基于井数据平滑阻抗模型约束，r'为前一道待反演参数反演结果，λ、β为权系数。

[0014] 优选的，通过基追踪算法获得所述待反演参数反射系数序列。

[0015] 优选的，通过对所述反演目标函数求取最小值获得所述待反演参数反射系数序列。

[0016] 优选的，根据所述待反演参数反射系数序列，获得最终待反演参数，其中，所述最终待反演参数包括纵波阻抗、横波阻抗、裂缝岩石物理参数。

[0017] 本发明的有益效果在于：本发明的基于基追踪横向多道约束的各向异性参数反演方法，通过对地震数据进行处理获得方位角度叠加道集，通过构建岩石物理模型，获取裂缝参数测井数据，对测井数据之间的协方差矩阵进行去相关处理，构建反演目标函数，从而获得待反演参数反射系数序列，通过纳入方位角获得方位角度叠加道集，方位角度叠加道集的信道比高，且包括多重地震信息，进而较为准确地提供储层中的裂缝信息，指导裂缝型储层的勘探和开发。

[0018] 本发明的基于基追踪横向多道约束的各向异性参数反演方法，通过对测井数据进行去相关处理，提高了反演结果的准确度，为地下裂缝预测提供较为准确的信息。

[0019] 本发明的基于基追踪横向多道约束的各向异性参数反演方法，通过在目标函数中加入相关约束项，构建反演目标函数，再对反演目标函数求取获得待反演参数反射系数序列，提高了反演结果稳定性和横向连续性，以此降低未知参数估测的不确定性，为地下裂缝预测提供合理可靠且有力的证据，解决了裂缝储层各向异性参数反演中存在的不稳定、横向不连续等问题，对于指导裂缝有效识别具有显著意义。

[0020] 本发明具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施例中将是显而易见的，或者将在并入本文中的附图和随后的具体实施例中进行详细陈述，这些附图和具体实施例共同用于解释本发明的特定原理。

附图说明

[0021] 通过结合附图对本发明示例性实施方式进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施方式中，相同的参考标号通常代表相同部件。

[0022] 图1示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的流程图。

[0023] 图2示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的方位角小角度叠加数据。

[0024] 图3示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的方位角中角度叠加数据。

[0025] 图4示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的方位角大角度叠加数据。

[0026] 图5示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的纵波阻抗反演结果。

[0027] 图6示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的横波阻抗反演结果。

[0028] 图7示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的裂缝岩石物理参数Γx反演结果。

[0029] 图8示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的裂缝岩石物理参数Γy反演结果。

[0030] 图9a示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的纵波阻抗IP反演结果对比。

[0031] 图9b示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的横波阻抗IS反演结果对比。

[0032] 图9c示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的裂缝岩石物理参数Γx反演结果对比。

[0033] 图9d示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的裂缝岩石物理参数Γy反演结果对比。

[0034] 图10a示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的纵波阻抗Ip横向多道对比。

[0035] 图10b示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的横波阻抗IS横向多道对比。

[0036] 图10c示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的裂缝岩石物理参数Γx横向多道对比。

[0037] 图10d示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的裂缝岩石物理参数Γy横向多道对比。

具体实施方式

[0038] 下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。

[0039] 根据本发明的基于基追踪横向多道约束的各向异性参数反演方法，包括：处理地震数据，获取方位角度叠加道集；构建裂缝岩石物理模型，获取裂缝参数测井数据；对测井数据进行去相关处理；构建反演目标函数；基于反演目标函数，获得待反演参数反射系数序列。

[0040] 具体地，对地震数据进行处理获得方位角度叠加道集，通过构建岩石物理模型，获取裂缝参数测井数据，测井数据包括：纵波阻抗、横波阻抗、裂缝岩石物理参数，再对测井数据之间进行去相关处理，构建反演目标函数，从而获得待反演参数反射系数序列。

[0041] 根据示例性的实施方式基于基追踪横向多道约束的各向异性参数反演方法，通过纳入方位角获得方位角度叠加道集，方位角度叠加道集的信道比高，且包括多重地震信息，进而较为准确地提供储层中的裂缝信息，指导裂缝型储层的勘探和开发。

[0042] 作为优选方案，对测井数据进行去相关处理包括：生成测井数据的协方差矩阵；对协方差矩阵进行奇异值分解，获得去相关矩阵；基于去相关矩阵，对公式d＝Gr进行变换，获得变换后的去相关公式，其中，d为方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数。

[0043] 作为优选方案，变换后的去相关公式为：d＝G×r，其中，G＝G×V，r＝V-1×r，d为方位角度叠加道集的地震数据，G为去相关的子波矩阵，G为子波矩阵，V为去相关矩阵，r为去相关的待反演参数，r为待反演参数。

[0044] 具体的，由于地层的纵波阻抗、横波阻抗和裂缝岩石物理参数之间是统计相关的，需要应用测井数据参数之间的协方差矩阵对测井数据进行去相关处理。首先将测井数据的协方差矩阵表示成：

[0045]

[0046] 其中： 为纵波阻抗反射系数的方差， 为纵波阻抗反射系数与横波阻抗反射系数之间的协方差， 为纵波阻抗反射系数与岩石物理参数之间的协方差，σIs2为横波阻抗反射系数的方差， 为横波阻抗反射系数与岩石物理参数之间的协方差， 为岩石物理参数之间的协方差， 为岩石物理参数的方差。

[0047] 一般来说，由于这四个变量间相互不独立，协方差矩阵Cm的非对角线上的元素不为零，对该协方差矩阵进行奇异值分解：

[0048]

[0049] 其中， 为纵波阻抗的方差， 为横波阻抗的方差， 为岩石物理参数的方差。

[0050] 假设 在假定反射系数不变且独立的情况下将其扩展为N个时间采样，则去相关矩阵V：

[0051]

[0052] 对公式d＝Gr进行变换：

[0053]

[0054] 得：

[0055] d＝G×r   (5)

[0056] 其中，d为方位角度叠加道集的地震数据，G为子波矩阵，G为去相关的子波矩阵，V为去相关矩阵，r为待反演参数，r为去相关的待反演参数。

[0057] 由公式(5)可以得出，去相关处理后的d方位角度叠加道集的地震数据与原来的d方位角度叠加道集的地震数据相同，因此，去相关处理对方位角度叠加道集的地震数据没有影响。

[0058] 所以，经过变换后的四参数协方差矩阵Cm'等于：

[0059]

[0060] 可以看出，Cm'非对角线上的元素都为零，表明经过变换后的四个参数之间是相互独立的。

[0061] 根据示例性的实施方式的基于基追踪横向多道约束的各向异性参数反演方法，通过对测井数据进行去相关处理，提高了反演结果的准确度，为地下裂缝预测提供较为准确的信息。

[0062] 作为优选方案，构建反演目标函数包括：获得所述方位角度叠加道集的似然函数；基于似然函数，获得目标函数；在目标函数中加入相关约束项，构建反演目标函数。

[0063] 具体的，假设地震资料背景噪声服从高斯分布，则可得到方位角度叠加道集的似然函数：

[0064]

[0065] 其中，d为所述方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数，σn2是噪音的方差。

[0066] 基于似然函数获得目标函数，同时将去相关处理后获得的参数代入目标函数中。

[0067] 通常通过对目标函数进行反演从而获得可靠的待反演参数，但是仅靠地震数据与反演结果间的误差最小二乘拟合不能得到精确反演解，在此基础上在目标函数中加入了相关约束项来增加反演稳定性以及横向连续性，从而减少反演结果的多解性，进而得到有意义的地层反射系数序列，最终获得可靠的待反演参数。

[0068] 根据示例性的实施方式的基于基追踪横向多道约束的各向异性参数反演方法，通过在目标函数中加入相关约束项，构建反演目标函数，再对反演目标函数求取获得待反演参数反射系数序列，提高了反演结果稳定性和横向连续性，以此降低未知参数估测的不确定性，为地下裂缝预测提供合理可靠且有力的证据，解决了裂缝储层各向异性参数反演中存在的不稳定、横向不连续等问题，对于指导裂缝有效识别具有显著意义。

[0069] 作为优选方案，基于似然函数，获得目标函数包括：结合似然函数与先验分布，获得后验概率分布；对后验概率分布进行最大化，获得目标函数。

[0070] 具体的，将似然函数与以前的数据经验结合起来得到后验概率分布，将后验概率分布最大化，得到目标函数：

[0071] J(r)＝(d-Gr)T(d-Gr)   (8)

[0072] 其中，d为所述方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数。

[0073] 作为优选方案，相关约束项包括：基于井数据平滑阻抗模型约束、L1范数约束和横向多道约束。

[0074] 具体的，基于井数据平滑阻抗模型约束是在有井资料时，通过数学变换把测井信息加入到反演方程中进一步约束反射系数，使反演问题变得更加适定，更有意义，同时，也可以地提高了反演的精度和分辨率。

[0075] 地下地层反射系数为稀疏序列，通过加入L1范数约束使反射系数稀疏化，反演结果更稳定。

[0076] 在地层平稳沉积中，地层横向基本连续，表现在地震反演中即道与道之间具有相关性，通过加入道与道之间的约束，即横向多道约束，来提高横向连续性。

[0077] 作为优选方案，反演目标函数为：

[0078]

[0079] 其中， d为方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数，Jsmooth(r)为基于井数据平滑阻抗模型约束，r'为前一道待反演参数反演结果，λ、β为权系数，αp为纵波阻抗反射系数约束对应的系数项，αs为横波阻抗反射系数约束对应的系数项， 为裂缝岩石物理参数反射系数约束对应的系数项，Cp为纵波阻抗对应的积分矩阵，Cs为横波阻抗对应的积分矩阵， 为裂缝岩石物理参数对应的积分矩阵，ξP为纵波阻抗初始模型，ξs为横波阻抗初始模型， 为裂缝岩石物理参数初始模型，E为单位矩阵，ξ'为前一道待反演参数反演结果。

[0080] 具体的，Jsmooth(r)为基于井数据平滑阻抗模型约束，λ||r||1为L1范数约束，为横向多道约束。

[0081] 在地震数据与反演结果间的误差最小二乘拟合下，使用公式r＝V-1×r，其中，V为去相关矩阵，r为去相关的待反演参数，r为待反演参数，把测井数据的纵波阻抗、横波阻抗、裂缝岩石物理参数这四个待反演参数乘以V-1，获得去相关处理后的待反演参数加入到反演目标函数中进一步约束反射系数。

[0082] 以纵波阻抗为例，一般有：

[0083]

[0084] 其中，Ip(t)为纵波阻抗，ΔIp(t)为上下地层纵波阻抗差，t为某个采样点。

[0085] 对反射系数积分有：

[0086]

[0087] 其中，t＝[t0,...,tM-1]T，Ip(t0)为初始纵波波阻抗值，rp(η)为反射系数。

[0088] 写成矩阵形式，即：

[0089] ξp＝Cprp   (12)

[0090] 其中，rp为反射系数， 且 为积分算子矩阵，其离散形式可表示为

[0091]

[0092] 所以基于井数据平滑阻抗模型约束项可表示为：

[0093]

[0094] 其中，αp为纵波阻抗反射系数约束对应的系数项，αs为横波阻抗反射系数约束对应的系数项， 为裂缝岩石物理参数反射系数约束对应的系数项，Cp为纵波阻抗对应的积分矩阵，Cs为横波阻抗对应的积分矩阵， 为裂缝岩石物理参数对应的积分矩阵，ξP为纵波阻抗初始模型，ξs为横波阻抗初始模型， 为裂缝岩石物理参数初始模型。

[0095] 反演目标函数可表示为：

[0096]

[0097] 其中， d为方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数，Jsmooth(r)为基于井数据平滑阻抗模型约束，r'为前一道待反演参数反演结果，λ、β为权系数，αp为纵波阻抗反射系数约束对应的系数项，αs为横波阻抗反射系数约束对应的系数项， 为裂缝岩石物理参数反射系数约束对应的系数项，Cp为纵波阻抗对应的积分矩阵，Cs为横波阻抗对应的积分矩阵， 为裂缝岩石物理参数对应的积分矩阵，ξP为纵波阻抗初始模型，ξs为横波阻抗初始模型， 为裂缝岩石物理参数初始模型，E为单位矩阵，ξ'为前一道待反演参数反演结果。。

[0098] 作为优选方案，通过基追踪算法获得待反演参数反射系数序列。

[0099] 作为优选方案，通过对反演目标函数求取最小值获得待反演参数反射系数序列。

[0100] 作为优选方案，根据待反演参数反射系数序列，获得最终待反演参数，其中，最终待反演参数包括纵波阻抗、横波阻抗、裂缝岩石物理参数。

[0101] 具体的，对得到的目标函数求取最小值得到待反演参数反射系数序列，进而推导得到最终待反演参数，其中包括纵波阻抗、横波阻抗、裂缝岩石物理参数。

[0102] 实施例

[0103] 图1示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的流程图。

[0104] 如图1所示，基于基追踪横向多道约束的各向异性参数反演方法，包括：

[0105] S102：处理地震数据，获取方位角度叠加道集；

[0106] S104：构建裂缝岩石物理模型，获取裂缝参数测井数据；

[0107] 根据实际工期及经验公式构建裂缝岩石物理模型，基于裂缝岩石物理模型获取裂缝参数与弹性参数之间的关系，得到裂缝参数的响应，即获得了裂缝参数测井数据，其中，测井数据包括：纵波阻抗、横波阻抗、裂缝岩石物理参数。

[0108] S106：对测井数据进行去相关处理；

[0109] 其中，对测井数据进行去相关处理包括：生成测井数据的协方差矩阵；对协方差矩阵进行奇异值分解，获得去相关矩阵；基于去相关矩阵，对公式d＝Gr进行变换，获得变换后的去相关公式，其中，d为方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数。

[0110] 其中，变换后的去相关公式为：d＝G×r，其中，G＝G×V，r＝V-1×r，d为方位角度叠加道集的地震数据，G为去相关的子波矩阵，G为子波矩阵，V为去相关矩阵，r为去相关的待反演参数，r为待反演参数。

[0111] 具体的，由于地层的纵波阻抗、横波阻抗和裂缝岩石物理参数之间是统计相关的，需要应用测井数据参数之间的协方差矩阵对测井数据进行去相关处理。首先将测井数据的协方差矩阵表示成：

[0112]

[0113] 其中： 为纵波阻抗反射系数的方差， 为纵波阻抗反射系数与横波阻抗反射系数之间的协方差， 为纵波阻抗反射系数与岩石物理参数之间的协方差，σIs2为横波阻抗反射系数的方差， 为横波阻抗反射系数与岩石物理参数之间的协方差， 为岩石物理参数之间的协方差， 为岩石物理参数的方差。

[0114] 一般来说，由于这四个变量间相互不独立，协方差矩阵Cm的非对角线上的元素不为零，对该协方差矩阵进行奇异值分解：

[0115]

[0116] 其中， 为纵波阻抗的方差， 为横波阻抗的方差， 为岩石物理参数的方差。

[0117] 假设 在假定反射系数不变且独立的情况下将其扩展为N个时间采样，则去相关矩阵V：

[0118]

[0119] 对公式d＝Gr进行变换：

[0120]

[0121] 得：

[0122] d＝G×r   (5)

[0123] 其中，d为方位角度叠加道集的地震数据，G为子波矩阵，G为去相关的子波矩阵，V为去相关矩阵，r为待反演参数，r为去相关的待反演参数。

[0124] 由公式(5)可以得出，去相关处理后的d方位角度叠加道集的地震数据与原来的d方位角度叠加道集的地震数据相同，因此，去相关处理对方位角度叠加道集的地震数据没有影响。

[0125] 所以，经过变换后的四参数协方差矩阵Cm'等于：

[0126]

[0127] 可以看出，Cm'非对角线上的元素都为零，表明经过变换后的四个参数之间是相互独立的。

[0128] S108：构建反演目标函数；

[0129] S108构建反演目标函数包括子步骤S1081-S1084：

[0130] S1081：假设地震资料背景噪声服从高斯分布，则可得到方位角度叠加道集的似然函数：

[0131]

[0132] 其中，d为方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数，σn2是噪音的方差。

[0133] S1082：将似然函数与先验分布结合起来得到后验概率分布，将后验概率分布最大化，得到目标函数：

[0134] J(r)＝(d-Gr)T(d-Gr)   (8)

[0135] S1083：将去相关处理后获得的参数代入目标函数。

[0136] S1084：在目标函数中加入相关约束项，构建反演目标函数。

[0137] 一般来说，仅靠地震数据与反演结果间的误差最小二乘拟合不能得到精确反演解，在此基础上本发明在目标函数中加入了相关约束项来增加反演稳定性以及横向连续性，从而减少反演结果的多解性，进而得到有意义的地层反射系数序列，最终获得可靠的待反演参数。这里的相关约束项包括基于井数据平滑阻抗模型约束、L1范数约束和横向多道约束，最终获得的反演目标函数为：

[0138]

[0139] 其中Jsmooth(r)为基于井数据平滑阻抗模型约束项，λ||r||1为L1范数约束，为横向多道约束，r'为前一道待反演参数反演结果，λ、β为权系数。

[0140] 具体的，使用公式r＝V-1×r，其中，V为去相关矩阵，r为去相关的待反演参数，r为待反演参数，把测井数据的纵波阻抗、横波阻抗、裂缝岩石物理参数这四个待反演参数乘以V-1，获得去相关处理后的待反演参数加入到反演目标函数中进一步约束反射系数。

[0141] 以纵波阻抗为例，一般有：

[0142]

[0143] 对反射系数积分有：

[0144]

[0145] 其中，t＝[t0,...,tM-1]T，Ip(t0)为初始纵波波阻抗值。写成矩阵形式，即：

[0146] ξp＝Cprp   (12)

[0147] 其中， 且 为积分算子矩阵，其离散形式可表示为

[0148] 所以基于井数据平滑阻抗模型约束项可表示为：

[0149]

[0150] 反演目标函数可表示为：

[0151]

[0152] 其中， d为方位角度叠加道集的地震数据，G为子波矩阵，r为待反演参数，Jsmooth(r)为基于井数据平滑阻抗模型约束，r'为前一道待反演参数反演结果，λ、β为权系数，αp为纵波阻抗反射系数约束对应的系数项，αs为横波阻抗反射系数约束对应的系数项， 为裂缝岩石物理参数反射系数约束对应的系数项，Cp为纵波阻抗对应的积分矩阵，Cs为横波阻抗对应的积分矩阵， 为裂缝岩石物理参数对应的积分矩阵，ξP为纵波阻抗初始模型，ξs为横波阻抗初始模型， 为裂缝岩石物理参数初始模型，E为单位矩阵，ξ'为前一道待反演参数反演结果。

[0153] S110：基于反演目标函数，获得待反演参数反射系数序列。

[0154] 通过基追踪算法对反演目标函数求取最小值获得待反演参数反射系数序列。

[0155] 最后的步骤，根据待反演参数反射系数序列，获得最终待反演参数，其中，最终待反演参数包括纵波阻抗、横波阻抗、裂缝岩石物理参数。

[0156] 该基于基追踪横向多道约束的各向异性参数反演方法的工作过程如下：对地震数据进行处理得到方位角度叠加道集，根据工区构建岩石物理模型进而计算获取裂缝参数测井数据，由于地层的纵波阻抗、横波阻抗和裂缝岩石物理参数之间是统计相关的，应用它们之间的协方差矩阵对测井数据进行去相关处理，去相关处理时首先生成待反演的弹性参数及裂缝岩石物理参数之间的协方差矩阵，对该协方差矩阵进行奇异值分解，获得去相关矩阵V，再对公式d＝Gr进行变换，得到变换后的去相关公式d＝G×r，可知去相关处理对方位角度叠加道集的地震数据没有影响；基于方位角度叠加道集生成似然函数，将似然函数结合以前的数据经验，获得后验概率分布，对后验概率分布进行最大化，从而获得目标函数，然后将测井数据参数乘以V-1，获得去相关处理后的参数，将去相关后的参数代入目标函数，再将基于井数据平滑阻抗模型约束、L1范数约束和横向多道约束的参数也加入目标函数，最后获得反演目标函数，通过基追踪算法对反演目标函数求取最小值获得待反演参数反射系数序列，获得最终待反演参数，其中，最终待反演参数包括纵波阻抗、横波阻抗、裂缝岩石物理参数。

[0157] 图2、图3和图4分别示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的方位角小角度叠加数据、方位角中角度叠加数据、方位角大角度叠加数据。在图2-4中，横坐标表示地震道数，纵坐标表示采样点数。

[0158] 如图2、图3、图4所示，对地震数据信息通过纳入方位角获得的方位角度叠加道集，该方位角度叠加道集反应的是真实的地震数据。

[0159] 图5、图6分别示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的纵波阻抗IP反演结果、横波阻抗IS反演结果；图7、图8分别示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的裂缝岩石物理参数Γx反演结果、裂缝岩石物理参数Γy反演结果。在图5-8中，横坐标表示地震道数，纵坐标表示采样点数。

[0160] 如图5、图6、图7和图8所示，通过基于基追踪横向多道约束的各向异性参数反演方法得到的纵波阻抗IP反演结果、横波阻抗IS反演结果、裂缝岩石物理参数Γx反演结果，以及裂缝岩石物理参数Γy反演结果，由图可以看出纵波阻抗IP反演结果、横波阻抗IS反演结果、裂缝岩石物理参数Γx反演结果，以及裂缝岩石物理参数Γy反演结果都与地震数据吻合，基本能反映裂缝发育程度。

[0161] 图9a、图9b、图9c和图9d分别示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的纵波阻抗IP、横波阻抗IS、裂缝岩石物理参数Γx、裂缝岩石物理参数Γy的反演结果对比。

[0162] 图10a、图10b、图10c和图10d分别示出了根据本发明的一个实施例的基于基追踪横向多道约束的各向异性参数反演方法的纵波阻抗Ip、横波阻抗IS、裂缝岩石物理参数Γx、裂缝岩石物理参数Γy的横向多道对比。

[0163] 如9a所示，通过基于基追踪横向多道约束的各向异性参数反演方法抽取的纵波阻抗IP反演结果对比，黑色点线为反演结果，灰色直线为模型数据，横坐标的数值是纵波阻抗(kg/m3)·(m/s)，由图可以看出纵波阻抗Ip反演结果与模型数据比较吻合。如9b所示，通过基于基追踪横向多道约束的各向异性参数反演方法抽取的横波阻抗IS反演结果对比，黑色3
点线为反演结果，灰色直线为模型数据，横坐标的数值是横波阻抗(kg/m)·(m/s)，由图可以看出横波阻抗IS反演结果稍不理想，这也是下一步需要改进的地方。如9c所示，通过基于基追踪横向多道约束的各向异性参数反演方法抽取的裂缝岩石物理参数Γx反演结果对比，黑色点线为反演结果，灰色直线为模型数据，由图可以看出裂缝岩石物理参数Γx反演结果与模型数据比较吻合。如9d所示，通过基于基追踪横向多道约束的各向异性参数反演方法抽取的裂缝岩石物理参数Γy反演结果对比，黑色点线为反演结果，灰色直线为模型数据，由图可以看出裂缝岩石物理参数Γy反演结果与模型数据比较吻合。在图9a，纵坐标表示采样点数，横坐标表示纵波阻抗(kg/m2*s)，在图9b，纵坐标表示采样点数，横坐标表示横波阻抗(kg/m2*s)，在图9c，纵坐标表示采样点数，横坐标表示岩石物理参数τx，在图9d，纵坐标表示采样点数，横坐标表示岩石物理参数τy。

[0164] 图10a所示，通过基于基追踪横向多道约束的各向异性参数反演方法抽取的纵波阻抗Ip横向多道对比，黑色点线为反演结果，灰色直线为模型数据，横坐标的数值是纵波阻抗(kg/m3)·(m/s)，由图可以看出纵波阻抗IP横向上反演结果也吻合模型数据，横向连续性也得到了提高。图10b所示，通过基于基追踪横向多道约束的各向异性参数反演方法抽取的横波阻抗IS横向多道对比，黑色点线为反演结果，灰色直线为模型数据，横坐标的数值是横波阻抗(kg/m3)·(m/s)，由图可以看出横波阻抗IS横向上反演结果也吻合模型数据，横向连续性也得到了提高。如10c所示通过基于基追踪横向多道约束的各向异性参数反演方法抽取的裂缝岩石物理参数Γx横向多道对比，黑色点线为反演结果，灰色直线为模型数据，由图可以看出裂缝岩石物理参数Γx横向上反演结果也吻合模型数据，横向连续性也得到了提高。如图10d所示，通过基于基追踪横向多道约束的各向异性参数反演方法抽取的裂缝岩石物理参数Γy横向多道对比，黑色点线为反演结果，灰色直线为模型数据，由图可以看出裂缝岩石物理参数Γy横向上反演结果也吻合模型数据，横向连续性也得到了提高。在图10a，纵坐标表示纵波阻抗(kg/m2*s)，横坐标表示采样点数，在图10b，纵坐标表示横波阻抗(kg/m2*s)，横坐标表示采样点数，在图10c，纵坐标表示岩石物理参数τx，横坐标表示采样点数，在图10d，纵坐标表示岩石物理参数τy，横坐标表示采样点数。

[0165] 以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。