[发明专利]HTI介质AVAZ定量反演方法及系统

说明书

公开了一种HTI介质AVAZ定量反演方法及系统。该方法可以包括：步骤1：基于各向异性项和方位角建立包括地震道N个时间采样的初始模型；步骤2：根据对应于不同入射角的不同方位角的地震记录差，获得观测数据；步骤3：根据方位角反射系数差矩阵与地震子波，获得地震合成记录，进而获得地震合成数据；步骤4：根据观测数据与地震合成数据，建立反演目标函数F(m+Δm)并进行反演；步骤5：当F(m+Δm)不满足设定精度要求时，将m修正为m+Δm，重复步骤3‑5；满足设定精度要求时，停止迭代，将迭代模型作为反演结果输出。本发明利用振幅的方位差异，对不同观测方向纵波反射系数差异定量反演各向异性项和方位角，通过定量反演提高裂缝预测精度。

技术领域

本发明涉及油气地球物理技术领域，更具体地，涉及一种HTI介质AVAZ定量反演方法及系统。

背景技术

裂缝作为油气运移的通道或者存储的空间，是储层甜点评价的重要参数之一。裂缝的存在导致地下介质呈现各向异性，同时由于上覆载荷的压实作用，在储层中主要分布垂直或者高角度裂缝，可等效为HTI(a transversely isotropic medium with ahorizontal axis of symmetry)介质。地震识别技术在裂缝预测方面发挥了越来越大的作用。早期人们依据裂缝引起的横波分裂特征进行裂缝描述；但是，多波多分量地震勘探的成本高昂，处理过程复杂，为了经济且有效地描述裂缝分布，展开运用纵波速度、振幅和旅行时等信息进行裂缝研究，其中利用纵波方位反射振幅信息即AVAZ(Amplitude variationwith incident angle and azimuth)裂缝预测的精度最高。根据调研，目前根据反射系数线性近似公式开展定量反演的研究比较少，薛娇于2016年发表的论文“基于等效介质模型的裂缝参数AVOA反演”能够直接对任意一点的地震时间序列进行定量预测的方法均直接对反射系数近似公式定量反演，但是，该方法涉及的参数多，速度、密度或者阻抗等参数的存在容易引起各向异性项预测结果的误差。因此，有必要开发一种HTI介质AVAZ定量反演方法及系统。

公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。

发明内容

本发明提出了一种HTI介质AVAZ定量反演方法及系统，其能够利用振幅的方位差异，对不同观测方向纵波反射系数差异定量反演各向异性项和方位角，通过定量反演提高裂缝预测精度。

根据本发明的一方面，提出了一种HTI介质AVAZ定量反演方法。所述方法可以包括：步骤1：基于各向异性项B_ani和方位角建立包括地震道N个时间采样的初始模型m；步骤2：根据对应于不同入射角的不同方位角的地震记录差，获得观测数据；步骤3：根据方位角反射系数差矩阵与地震子波，获得地震合成记录，进而获得地震合成数据；步骤4：根据所述观测数据与所述地震合成数据，建立反演目标函数F(m+Δm)并进行反演，其中，Δm是模型修正量；步骤5：当F(m+Δm)不满足设定精度要求时，将m修正为m+Δm，重复步骤3-5；满足设定精度要求时，停止迭代，将所述迭代模型作为反演结果输出。

优选地，所述初始模型m为：

其中，为第n个采样时间点的初始模型，n＝1，…，N，N为采样时间点的数量。

优选地，第n个采样时间点的初始模型m_n为：

其中，B_ani为各向异性项，为第n个采样时间点的方位角。

优选地，入射角为θ_i时的方位角记录差为：