[发明专利]针对激发井深设计的微测井层析成像方法

说明书

本发明提供一种针对激发井深设计的微测井层析成像方法，包括：输入微测井信息、微测井初至走时数据、微测井层位解释信息；确定反演深度和反演网格；基于确定的反演网格，对微测井层位解释信息进行网格化，获得与模型尺度相一致的解释约束项；构建方程组，并添加速度平滑项以及微测井解释约束项；求解方程组，获得反演结果；对结果重采样，获得符合要求的速度模型。该针对激发井深设计的微测井层析成像方法降低了人为影响，所得速度分布包含更多的近地表构造细节，便于得到准确的最佳激发深度，确保地震数据的采集质量，为之后的地震数据处理奠定了基础。

技术领域

本发明涉及气勘探地震资料处理技术领域，特别是一种针对激发井深设计的微测井层析成像方法。

背景技术

近地表问题是陆上地震勘探最常见的问题之一。为研究近地表结构，近地表调查(小折射、微测井等)是目前常规的生产环节。微测井解释方法能为地震数据采集选择最佳激发岩性和井深，确保地震数据的采集质量。此外，利用微测井资料还可得到工区内表层速度和厚度在横向上和纵向上的变化规律，为室内地震数据处理提供静校正量。对于山前带等复杂地区的地震勘探，如何计算近地表速度分布，并有效确定表层厚度、激发井深等，成为十分重要的研究课题，而近地表调查发挥了举足轻重的作用。

微测井常规解释首先钻井穿过低降速带，从井底到井口依次激发，根据接收到的透射地震波记录，拾取地震波初至走时并转换成垂向时间，与对应的深度绘在时深坐标系内，根据直达波在地层中的传播规律，若不同深度点位于同一速度层内，时深曲线为一直线，不同速度层对应的直线斜率不同，根据其分布规律，划分出各层的位置，直线的斜率为介质的层速度，两直线的交点为介质的分界面。拟合时深曲线可得层速度及其厚度。

常规微测井解释过程中，层位的划分需要凭借处理人员的视觉和经验来判定，参数选取的正确性和可靠性受人为因素的影响大，难以真实客观的评价处理方法和参数的选择。结果只能粗略描述近地表结构，难以描述复杂地表的细节，不易准确给定激发井深，降低近地表调查结果的可信度，这制约了微测井技术的应用。为此，我们发明了一种针对激发井深设计的微测井层析成像方法，解决了以上技术问题。

发明内容

本发明的目的是提供一种针对微测井解释人工分层的问题，创新提出了微测井层析成像，对近地表速度分布进行有效评价的针对激发井深设计的微测井层析成像方法。

本发明的目的可通过如下技术措施来实现：针对激发井深设计的微测井层析成像方法，该针对激发井深设计的微测井层析成像方法包括：步骤一，输入微测井信息、微测井采集的深度与初至走时数据、微测井层位解释信息；步骤二，确定最大反演深度和反演速度网格，设定反演的模型参数S；步骤三，基于已定的反演网格，对微测井层位解释信息进行网格化，获得与模型尺度相一致的解释约束项V_mdl；步骤四，构建方程组，并添加速度平滑项以及微测井解释约束项；步骤五，应用SIRT算法求解方程组，获得反演结果；步骤六，对结果重采样，获得符合要求的速度模型。

本发明的目的还可通过如下技术措施来实现：

在步骤一中，输入的微测井信息包含站号、横纵坐标、高程，微测井深度与初至走时数据包含微测井各点的激发深度和初至走时信息，微测井层位解释信息包括各层位的深度及速度。

在步骤一中，将初至的走时斜率转化为速度，并通过设定地表最小速度和最大速度，剔除速度过低或过高的微测井数据。

在步骤二中，反演所采用的深度大于所用微测井的最大采集深度。

在步骤二中，设定初始模型参数为初至的走时斜率。

在步骤三中，为实现微测井层位解释信息的网格化，从地表到最大反演深度逐步循环，当网格深度位于某解释层内时，采用该层的解释速度作为解释约束项。

在步骤四中，构建的方程组如下：