[发明专利]火成岩约束速度建模方法、系统、计算机设备及存储介质

说明书

本申请涉及一种火成岩约束速度建模方法、系统、计算机设备和存储介质，该方法包括：利用火成岩地震相获取不同尺度火成岩的空间分布特征信息；将空间分布特征作为约束条件，通过常规层析反演的目标泛函中构建初始目标泛函，其中，空间分布特征用于对常规层析反演的目标泛函进行地震相约束；将从地震相中提取的火成岩区域速度区间信息应用于不等式正则化项，通过初始目标泛函得到新目标泛函，建立不等式约束的地震相控层析反演新目标泛函；将新目标泛函应用到层析反演过程中，以基于不等式正则化的地震相控火成岩速度建模；能够实现对不同岩性火成岩速度异常体的合理控制，提高火成岩建模的精度与稳定性。

技术领域

本申请涉及地球物理勘探领域，特别是涉及到一种火成岩约束速度建模方法、系统、计算机设备及存储介质。

背景技术

目前，在地震勘探领域，工业界中最为常用的方法依旧是基于射线的层析速度建模技术。即便是简单的基于射线理论近似的层析方法，依旧需要大量的人力、物力来构建初始速度并更新，而且构建的速度为低波数背景，精度较低，对于不同尺度、特别是小尺度的火成岩速度建模精度不够，反演结果较差，从而导致最终的火成岩区域成像精度较差。

基于此，不管是学术界还是工业界，诸多学者采取了一系列的手段提高火成岩区域速度建模精度。其中最为常用的方法是将已成像的火成岩进行人工雕刻，并根据经验填充速度，这种方法虽然看似有效，但过分依赖于处理解释人员的主观认识，一旦出现错误，往往会出现虚假构造，我们通常称这种方法为“硬约束”下的火成岩速度建模技术。相比较而言，另一种方法是在层析过程中引入构造约束，通过已有的数据进行正则化构造约束，完全基于数据驱动，反演火成岩区域的速度模型，可以对应成为“软约束”火成岩速度建模。这种方法虽然完全尊重数据本身，避免了处理人员的主观人事，但单纯的数据驱动并不能很好地反演小尺度的火成岩速度，火成岩速度建模效果不尽人意。

发明内容

基于此，有必要针对上述技术问题，提供一种能够实现对不同岩性火成岩速度异常体进行合理控制、提高火成岩建模的精度与稳定性的火成岩约束速度建模方法、系统、计算机设备及存储介质。

本发明的一方面公开了一种火成岩约束速度建模方法，所述方法包括：

利用火成岩地震相获取不同尺度火成岩的空间分布特征信息；

将所述空间分布特征作为约束条件，加入常规层析反演的目标泛函中构建初始目标泛函，其中，所述空间分布特征用于对所述常规层析反演的目标泛函进行地震相约束；

将从地震相中提取火成岩区域速度区间信息应用于不等式正则化项，通过所述初始目标泛函得到新目标泛函，建立不等式约束的地震相控层析反演所述新目标泛函；

将所述新目标泛函应用到层析反演过程中，形成基于不等式约束正则化的地震相控火成岩速度建模。

可选的，利用火成岩地震相获取不同尺度火成岩的空间分布特征信息，包括：

通过井震结合手段识别火成岩地震相模式，获得所述不同尺度火成岩的空间分布特征，其中，所述空间分布特征信息包括以下一项或两项：火成岩分布范围及构造顶底厚度。

可选的，常规层析反演的目标泛函包括：

其中，S(m)是道集拉平项，属于常规层析反演的准则。z^pick是拾取的反射点位置，z^true是真实的反射点位置。当层析速度准确时，S(m)最小。

可选的，所述初始目标泛函包括：

其中，W_mask是计算地震相得到的火成岩分布范围，Δv是上一轮的速度更新量。ε₁是约束权重因子，通常根据实际情况在0-1之间取值。