阅读说明：本技术 一种利用微测井数据构建品质因子体的方法 (method for constructing quality factor body by using micro-logging data ) 是由 刁瑞 王常波 尚新民 刘立彬 韩站一 冯玉苹 崔庆辉 于 2019-07-22 设计创作，主要内容包括：本发明公开了一种利用微测井数据构建品质因子体的方法,收集工区范围内的微测井测量数据,对微测井数据进行分组处理；在同一组微测井数据内,利用修正S域变换方法,对微测井数据进行二维S域频谱分析,计算得到每一道数据的主频、低频值、高频值、有效频宽、特定频带能量占比等参数属性,通过不同深度激发和不同位置接收的地震波数据计算主频、低频值、高频值、有效频宽、特定频带能量占比等参数属性,结合传播距离权值和参数属性比例,得到综合参数属性,最终通过克里金延拓插值方法将不同测量点位置的品质因子曲线内插为近地表体三维品质因子体,对野外施工井深设计、近地表吸收衰减补偿和高分辨率处理、真地表粘弹性介质偏移等具有重要作用。(The invention discloses a method for constructing a quality factor body by using micro-logging data, which collects the micro-logging measurement data in a work area range and carries out grouping processing on the micro-logging data; in the same set of micro-logging data, a correction S-domain transformation method is utilized to perform two-dimensional S-domain spectral analysis on the micro-logging data, the parameter attributes such as the dominant frequency, the low frequency value, the high frequency value, the effective bandwidth, the specific frequency band energy ratio and the like of each channel of data are obtained through calculation, the parameter attributes such as the dominant frequency, the low frequency value, the high frequency value, the effective bandwidth, the specific frequency band energy ratio and the like are calculated through seismic wave data excited at different depths and received at different positions, the comprehensive parameter attributes are obtained through combination of propagation distance weight and parameter attribute ratio, and finally, the quality factor curves at different measuring point positions are interpolated into near-surface body three-dimensional quality factor bodies through a kriging interpolation method, so that the method has important effects on field construction well depth design, near-surface absorption attenuation compensation and high resolution processing, real-surface viscoelastic medium migration.)

一种利用微测井数据构建品质因子体的方法

技术领域

本发明涉及油气田勘探技术领域，特别是一种利用微测井数据构建品质因子体的方法。

背景技术

我国大部分油田均属于典型的冲积平原，复杂近地表大致可分为低速层和降速层。由于复杂近地表低降速层的厚度变化和横向速度剧烈变化，地震波传播存在传播时间延迟和能量、频率剧烈衰减的特点。在以往的三维地震勘探中，对复杂近地表的问题研究不足，通常认为表层的潜水面基本是稳定的，横向变化相对较小，在进行野外数据采集阶段，全区仅选择典型的试验点位置，进行小折射或微测井测量，并采用谱比法计算品质因子，存在品质因子计算不稳地、数值差异大的问题，制约着复杂近地表品质因子体的准确构建，从而影响高品质原始数据的野外采集，并影响后续的地震资料处理，制约着三维地震精细勘探与开发。

由于复杂近地表纵横向速度和岩性变化都很大，存在剧烈的差异，需要从点到线，从线到面进行近地表精细探测，获得三维近地表品质因子体。但是由于现有技术的制约，复杂近地表品质因子体难于建立。

因此，亟待开发一种利用微测井数据构建品质因子体的方法。

发明内容

本发明的目的是要解决现有技术存在的不足，提供一种利用微测井数据构建品质因子体的方法。

为达到上述目的，本发明是按照以下技术方案实施的：

一种利用微测井数据构建品质因子体的方法，包括以下步骤：

S1、在工区范围内的浅井中不同深度激发，在地表采用多个检波器进行接收地震波数据即微测井数据，对微测井数据进行分组处理，每个激发深度的不同检波器接收数据为一组；

S2、在同一组微测井数据内，利用修正S域变换方法，对微测井数据进行二维S域频谱分析，计算得到每一道数据的如下参数属性：主频、低频值、高频值、有效频宽、特定频带能量占比；

S3、根据同一组内检波器的布设位置，进行传播距离赋权值，得到含有距离权值的参数属性，并进行不同检波器参数属性对比，剔除异常属性值；

S4、分别计算不同激发深度的主频、低频值、高频值、有效频宽、特定频带能量占比的参数属性，得到所有深度所对应的综合参数属性；

S5、根据激发深度推算传播距离，结合不同激发深度所对应的综合参数属性的变化，得到微测井测量位置处的品质因子曲线；

S6、计算所有微测井测量位置处的品质因子曲线后，采用克里金延拓插值方法，将工区范围内的测量位置散点内插为一个三维近地表品质因子体。

进一步，所述S1中，炮点激发深度用i表示，即i＝1，2，......，n，每个激发深度的不同检波器接收数据为一组，每个激发深度采用5道或9道或13道检波器进行接收，同一组数据内不同的检波器用j表示，即j＝1，2，......，m。

进一步，所述S2中，在S域信噪比曲线和S域纯度曲线的控制下，通过S域频谱求取主频Mf_i，j、低频值Lf_i，j、高频值Hf_i，j、有效频宽Bf_i，j、特定频带能量占比Ef_i，j参数属性，其中：i为炮点激发深度，i＝1，2，......，n；j为同一组数据内不同的检波器，j＝1，2，......，m；其中信噪比曲线是近地表数据在S域频谱内求取的有效信号与噪音干扰的比值；其中纯度曲线是近地表数据在S域频谱内求取的各个频带信号的纯度。

进一步，所述S3中，所述传播距离赋权值记为δ_j，其中：j＝1，2，......，m，表示不同的检波器；得到含有距离权值的参数属性分别为：主频δ_j×Mf_i，j、低频值δ_j×Lf_i，j、高频值δ_j×Hf_i，j、有效频宽δ_j×Bf_i，j、特定频带能量占比δ_j×Ef_i，j，其中：i为炮点激发深度，i＝1，2，......，n；j为同一组数据内不同的检波器，j＝1，2，......，m；对不同的检波器j＝1，2，......，m进行参数属性比较，剔除所有m个检波器中的异常属性值。

进一步，所述S4中，计算得到对应的综合参数属性的具体步骤为：进行计算综合参数属性时，每一种属性赋予不同的比例,综合参数属性其中：δ_j为传播距离权值，M_j为主频Mf_i，j的比例，L_j为低频值Lf_i，j的比例，H_j为高频值Hf_i，j的比例，B_j为有效频宽Bf_i，j的比例，E_j为特定频带能量占比Ef_i，j的比例，i＝1，2，......，n，j＝1，2，......，m。

进一步，在步骤5中，通过激发深度推算传播距离，根据不同激发深度的综合参数属性θ_i，i＝1，2，......，n，计算得到微测井测量位置处的品质因子曲线。

作为优选技术方案，在步骤6中，计算所有微测井测量位置处的品质因子曲线后，采用克里金延拓插值方法，将工区范围内的测量位置散点内插为一个三维体。在步骤5中已经计算得到了每个测量点位置处的品质因子曲线，根据克里金延拓插值方法，可以将变化曲线内插为三维数据体，即近地表体三维品质因子体。

与现有技术相比，本发明充分发挥微测井数据的特点，采用S域信噪比曲线、S域纯度曲线约束条件，通过不同深度激发和不同位置接收的地震波数据计算主频、低频值、高频值、有效频宽、特定频带能量占比等参数属性，结合传播距离权值和参数属性比例，得到综合参数属性，最终通过克里金延拓插值方法将不同测量点位置的品质因子曲线内插为近地表体三维品质因子体，对野外施工井深设计、近地表吸收衰减补偿和高分辨率处理、真地表粘弹性介质偏移等具有重要作用。

附图说明

图1为本发明的流程图。

图2为本发明实施例中复杂近地表三维品质因子体示意图。

具体实施方式

为使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步的详细说明。此处所描述的具体实施例仅用于解释本发明，并不用于限定发明。

参照图1-图2，本实施例的一种利用微测井数据构建品质因子体的方法，具体操作过程如下：

(1)收集工区范围内的微测井测量数据，微测井测量数据是在浅井中不同深度激发，在地表采用多个检波器进行接收的地震波数据，炮点激发深度用i表示，即i＝1，2，......，n，每个激发深度的不同检波器接收数据为一组，通常采用5道或9道或13道检波器进行接收，同一组数据内不同的检波器用j表示，即j＝1，2，......，m。

(2)对微测井数据进行分组处理，首先，不同测量点位置处的微测井数据为一组；其次，不同炮点激发深度的微测井数据为一组，接下来对每个激发深度的不同检波器接收数据进行综合参数属性求取，如图1所示。

(3)在相同激发深度的不同检波器微测井数据内，利用修正S域变换方法，对每一道微测井数据进行二维S域频谱分析，修正S域变换方法采用能量归一化变宽度窗函数，对微测井数据进行分时窗的频谱分析，在不同频率段采用不同宽度窗函数，计算得到随时间变化的频率特征，该变换方法具有分辨率自动调整、各频率段能量均衡的优点。通过S域频谱求取主频Mf_i，j、低频值Lf_i，j、高频值Hf_i，j、有效频宽Bf_i，j、特定频带能量占比Ef_i，j等参数属性，其中：i为炮点激发深度，i＝1，2，......，n；j为同一组数据内不同的检波器，j＝1，2，......，m。

(4)对微测井数据进行二维S域频谱分析，在S域信噪比曲线和S域纯度曲线的控制下，计算每一道数据的主频、低频值、高频值、有效频宽、特定频带能量占比等参数属性，其中信噪比曲线是近地表数据在S域频谱内求取的有效信号与噪音干扰的比值，其中纯度曲线是近地表数据在S域频谱内求取的各个频带信号的纯度。

(5)根据同一组内检波器的布设位置，进行传播距离赋权值δ_j，其中：j＝1，2，......，m,表示不同的检波器，传播距离的远近决定了权值δ_j的大小。

(6)进行传播距离赋权值后，这样就可以得到含有距离权值的参数属性，分别为主频δ_j×Mf_i，j、低频值δ_j×Lf_i，j、高频值δ_j×Hf_i，j、有效频宽δ_j×Bf_i，j、特定频带能量占比δ_j×Ef_i，j，其中：i为炮点激发深度，i＝1，2，......，n；j为同一组数据内不同的检波器，j＝1，2，......，m。

(7)进行不同检波器参数属性对比，即对不同的检波器j＝1，2，......，m进行参数属性比较，剔除所有m个检波器中的异常属性值，从而排除异常道的影响，提高参数属性的准确性和可靠度。

(8)异常属性剔除后进行综合参数属性计算，根据参数属性的重要性给每一种属性赋予不同的比例,综合参数属性其中：δ_j为传播距离权值，M_j为主频Mf_i，j的比例，L_j为低频值Lf_i，j的比例，H_j为高频值Hf_i，j的比例，B_j为有效频宽Bf_i，j的比例，E_j为特定频带能量占比Ef_i，j的比例，i＝1，2，......，n，j＝1，2，......，m。

(9)根据该激发深度的主频、低频值、高频值、有效频宽、特定频带能量占比等参数属性，计算得到了一个炮点激发深度的综合参数属性，如果还有炮点激发深度没有计算完成，则返回，对不同激发深度的微测井数据进行相同流程的计算。

(10)分别计算不同激发深度的主频Mf_i，j、低频值Lf_i，j、高频值Hf_i，j、有效频宽Bf_i，j、特定频带能量占比Ef_i，j等参数属性，并得到对应深度的综合参数属性θ_i，i＝1，2，......，n，如果所有炮点激发深度都已计算完成，则流程继续进行。

(11)通过激发深度推算传播距离，根据不同激发深度的综合参数属性θ_i，i＝1，2，......，n，计算得到微测井测量位置处的品质因子曲线。品质因子值的大小由激发深度和综合参数属性决定，根据综合参数属性在当前和下一个激发深度的变化情况，通过振幅和频率衰减关系公式，可以计算得到该激发深度的品质因子值，所有激发深度的品质因子值计算完成后，可以通过激发深度和品质因子值计算得到测量位置处的品质因子曲线。

(12)如果还有不同测量点位置的微测井数据没有计算完成，则返回，对不同测量点位置的微测井数据进行相同流程的计算，得到所有测量点位置处的品质因子曲线，如果所有测量点都已计算完成，则流程继续进行。

(13)计算所有微测井测量位置处的品质因子曲线后，采用克里金延拓插值方法，将工区范围内的测量位置散点内插为一个三维体，克里金延拓插值方法基于最小二乘算法，通过数学协方差函数进行加权赋值，实现三维空间内曲线的三维立体插值，从而完成由曲线到三维体的延拓建模。采用克里金延拓插值方法以将品质因子曲线内插为三维数据体，即近地表体三维品质因子体，如图2所示。

本发明的技术方案不限于上述具体实施例的限制，凡是根据本发明的技术方案做出的技术变形，均落入本发明的保护范围之内。