一种确定致密储层含气性的方法及装置
阅读：596发布：2021-01-29
专利汇可以提供一种确定致密储层含气性的方法及装置专利检索，专利查询，专利分析的服务。并且本 申请 实施例 公开了一种确定致密储层含气性的方法及装置。所述方法包括：基于 测井 数据，确定目的储层的正演 共中心点 道集数据；基于正演 共中心点道集 数据和实测共中心点道集数据，确定目标划分数量；其中，目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量；确定实测共中心点道集数据中的最小纵波入射 角 度和最大纵波入射角度；基于最小纵波入射角度和最大纵波入射角度，从实测共中心点道集数据中确定目标划分数量的目标分道集数据；其中，目标分道集数据与纵波入射角度范围相对应；根据目标分道集数据中的振幅值，确定目的储层的含气性。本申请实施例提供的技术方案，可以提高所确定的致密储层含气性的准确度。，下面是一种确定致密储层含气性的方法及装置专利的具体信息内容。
1.一种确定致密储层含气性的方法，其特征在于，提供有目的储层的实测共中心点道集数据，以及所述目的储层中钻井位置处的测井数据；所述方法包括：
基于所述测井数据，确定所述目的储层的正演共中心点道集数据；
基于所述正演共中心点道集数据和所述实测共中心点道集数据，确定目标划分数量；
其中，所述目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量；
确定所述实测共中心点道集数据中的最小纵波入射角度和最大纵波入射角度；
基于所述最小纵波入射角度和所述最大纵波入射角度，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据；其中，所述目标分道集数据与纵波入射角度范围相对应；
根据所述目标分道集数据中的振幅值，确定所述目的储层的含气性。
2.根据权利要求1所述的方法，其特征在于，确定所述目的储层的正演共中心点道集数据，包括：
基于所述实测共中心点道集数据，确定所述目的储层的地震子波数据；
基于所述测井数据，确定所述目的储层的弹性参数；
根据预设纵波入射角和所述弹性参数，确定所述目的储层的地震波反射系数；
对所述地震子波数据和所述地震波反射系数进行正演处理，得到所述目的储层的正演共中心点道集数据。
3.根据权利要求1所述的方法，其特征在于，确定目标划分数量，包括：
针对预设多个划分数量中指定划分数量，按照纵波入射角度，将所述正演共中心点道集数据中指定纵波入射角度范围内的道集数据平均划分为所述指定划分数量的正演分道集数据，以及将所述实测共中心点道集数据中所述指定纵波入射角度范围内的道集数据平均划分为所述指定划分数量的实测分道集数据；
分别将各个所述正演分道集数据中的采样点的振幅值之和作为一个正演振幅值，得到一组正演数据；所述正演数据中包括所述指定划分数量的正演振幅值和所述正演振幅值对应的平均纵波入射角度；以及分别将所述实测分道集数据中的采样点的振幅值之和作为一个实测振幅值，得到一组实测数据；所述实测数据中包括所述指定划分数量的实测振幅值和所述实测振幅值对应的平均纵波入射角度；其中，所述平均纵波入射角度表示所述正演分道集数据或所述实测分道集数据的纵波入射角度的平均值；
计算所述正演数据中指定平均纵波入射角度对应的正演振幅与所述实测数据中所述指定平均纵波入射角度对应的实测振幅之间的差值，以及所述差值的平方；
将各个所述平均纵波入射角度分别对应的差值之和的绝对值，以及各个所述平均纵波入射角度分别对应的差值的平方之和分别作为所述指定划分数量对应的第一对比因子和第二对比因子；当所述指定划分数量对应的第一对比因子最小、且对应的第二对比因子最小时，将所述指定划分数量作为所述目标划分数量。
4.根据权利要求1所述的方法，其特征在于，确定所述实测共中心点道集数据中的最小纵波入射角度，包括：
按照指定纵波入射角度间隔，将所述实测共中心点道集数据划分为多个分道集数据，并对所述分道集数据进行叠加处理；
按照纵波入射角度由小到大，依次计算叠加处理后的分道集数据的信噪比，直至所述叠加处理后的分道集数据的信噪比大于或等于指定信噪比阈值，并将信噪比大于或等于指定信噪比阈值的叠加处理后的分道集数据对应的纵波入射角度范围中的最大入射角度作为所述最小纵波入射角度。
5.根据权利要求1所述的方法，其特征在于，所述方法还提供有所述目的储层的深度；
确定所述实测共中心点道集数据中的最大纵波入射角度，包括：
获取所述实测共中心点道集数据对应的最大偏移距；
根据所述最大偏移距和所述目的储层的深度，确定所述最大纵波入射角度。
6.根据权利要求1所述的方法，其特征在于，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据，包括：
从所述实测共中心点道集数据中获取纵波入射角度在所述最小纵波入射角和所述最大纵波入射角之间的目标共中心点道集数据；
根据所述目标共中心点道集数据，确定振幅随纵波入射角度变化的目标曲线；其中，所述目标曲线中包括多个目标数据点，所述目标数据点中包括目标振幅值和所述目标振幅值对应的目标纵波入射角度；所述目标振幅值表示所述目标共中心点道集数据中所述目标纵波入射角度对应的地震道数据的采样点的振幅值之和；
确定所述目标曲线的拐点，并将所述拐点对应的纵波入射角度作为用于划分共中心点道集数据的目标临界角度；
根据所述目标临界角度，将所述目标共中心点道集数据划分为所述目标划分数量的目标分道集数据。
7.根据权利要求1所述的方法，其特征在于，所述目标划分数量大于或等于3，确定所述目的储层的含气性，包括：
分别对所述目标划分数量的目标分道集数据中第一分道集数据、第二分道集数据和第三分道集数据进行叠加处理；其中，所述第一分道集数据和所述第三分道集数据用于表示所述目标划分数量的目标分道集数据中纵波入射角最小和纵波入射角最大的分道集数据；
所述第二分道集数据的纵波入射角度小于所述第三分道集数据的纵波入射角度、且大于所述第一分道集数据的纵波入射角度；
基于叠加处理后的第一分道集数据的最大振幅值、叠加处理后的第二分道集数据的最大振幅值和叠加处理后的第三分道集数据的最大振幅值，确定所述目的储层的含气等级。
8.根据权利要求7所述的方法，其特征在于，确定所述目的储层的含气等级，包括：
将所述叠加处理后的第一分道集数据的最大振幅值减去叠加处理后的第三分道集数据的最大振幅值，得到第一振幅差值，以及将所述叠加处理后的第二分道集数据的最大振幅值减去叠加处理后的第三分道集数据的最大振幅值，得到第二振幅差值；
将所述第一振幅差值与所述第二振幅差值相乘，并根据乘积结果，确定所述目的储层的含气等级。
9.根据权利要求8所述的方法，其特征在于，确定所述目的储层的含气等级，包括：
当所述乘积结果大于或等于指定阈值时，确定所述目的储层的含气等级为第一等级；
当所述乘积结果小于所述指定阈值时，确定所述目的储层的含气等级为第二等级。
10.一种确定致密储层含气性的装置，其特征在于，所述装置提供目的储层的实测共中心点道集数据，以及所述目的储层中钻井位置处的测井数据；所述装置包括：正演数据确定模块、目标划分数量确定模块、入射角度确定模块、分道集数据确定模块和含气性确定模块；其中，
所述正演数据确定模块，用于基于所述测井数据，确定所述目的储层的正演共中心点道集数据；
所述目标划分数量确定模块，用于基于所述正演共中心点道集数据和所述实测共中心点道集数据，确定目标划分数量；其中，所述目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量；
所述入射角度确定模块，用于确定所述实测共中心点道集数据中的最小纵波入射角度和最大纵波入射角度；
所述分道集数据确定模块，用于基于所述最小纵波入射角度和所述最大纵波入射角度，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据；其中，所述目标分道集数据与纵波入射角度范围相对应；
所述含气性确定模块，用于根据所述目标分道集数据中的振幅值，确定所述目的储层的含气性。
一种确定致密储层含气性的方法及装置
技术领域
[0001] 本申请涉及石油地球物理勘探技术领域，特别涉及一种确定致密储层含气性的方法及装置。
背景技术
[0002] 随着岩性勘探不断深入，天然气勘探开发对于储层含气性预测精度要求越来越高，随着AVO(Amplitude Versus Offset，振幅随偏移距的变化)技术的推广，AVO相关技术在储层特征描述中的作用日益凸显，充分发掘了叠前地震资料的有用信息，有效的降低勘探开发的风险与成本。
[0003] AVO技术是利用地震反射波振幅与炮检距或入射角的关系寻找油气的一种方法，相关主要技术有：岩石物理研究、AVO正演模拟、AVO属性、叠前弹性反演等。其中，叠前弹性反演是地震资料处理、地震解释、测井分析及岩石物理模拟等多学科一体化技术。该技术保留了地震反射振幅随骗移距不同或纵波入射角度不同而变化的特征，其理论基础是佐普利兹(Zoeppritz)方程，通过简化该方程，获取多种简化形式，通过线性和非线性反演可以得到很多叠后反演无法得到的参数，尤其适用于物性和流体变化较敏感的储层。
[0004] 然而，由于叠前弹性反演需要高保真、高分辨率和高信噪比的叠前地震数据，以及可靠的纵波测井数据和横波测井数据，同时不同反演软件存在几十至上百个反演参数的选择，这样会造成反演成果受地震数据、测井数据和人为因素的影响较大，尤其是在鄂尔多斯盆地西南部、天环和盆地东部地震低信噪比地区，采用叠前弹性反演技术预测致密储层含气性符合率较低。发明内容
[0005] 本申请实施例的目的是提供确定致密储层含气性的方法及装置，以提高所确定的致密储层含气性的准确度。
[0006] 为解决上述技术问题，本申请实施例提供一种确定致密储层含气性的方法及装置是这样实现的：
[0008] 基于所述测井数据，确定所述目的储层的正演共中心点道集数据；
[0009] 基于所述正演共中心点道集数据和所述实测共中心点道集数据，确定目标划分数量；其中，所述目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量；
[0010] 确定所述实测共中心点道集数据中的最小纵波入射角度和最大纵波入射角度；
[0011] 基于所述最小纵波入射角度和所述最大纵波入射角度，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据；其中，所述目标分道集数据与纵波入射角度范围相对应；
[0012] 根据所述目标分道集数据中的振幅值，确定所述目的储层的含气性。
[0013] 优选方案中，确定所述目的储层的正演共中心点道集数据，包括：
[0014] 基于所述实测共中心点道集数据，确定所述目的储层的地震子波数据；
[0015] 基于所述测井数据，确定所述目的储层的弹性参数；
[0017] 对所述地震子波数据和所述地震波反射系数进行正演处理，得到所述目的储层的正演共中心点道集数据。
[0018] 优选方案中，确定目标划分数量，包括：
[0019] 针对预设多个划分数量中指定划分数量，按照纵波入射角度，将所述正演共中心点道集数据中指定纵波入射角度范围内的道集数据平均划分为所述指定划分数量的正演分道集数据，以及将所述实测共中心点道集数据中所述指定纵波入射角度范围内的道集数据平均划分为所述指定划分数量的实测分道集数据；
[0020] 分别将各个所述正演分道集数据中的采样点的振幅值之和作为一个正演振幅值，得到一组正演数据；所述正演数据中包括所述指定划分数量的正演振幅值和所述正演振幅值对应的平均纵波入射角度；以及分别将所述实测分道集数据中的采样点的振幅值之和作为一个实测振幅值，得到一组实测数据；所述实测数据中包括所述指定划分数量的实测振幅值和所述实测振幅值对应的平均纵波入射角度；其中，所述平均纵波入射角度表示所述正演分道集数据或所述实测分道集数据的纵波入射角度的平均值；
[0021] 计算所述正演数据中指定平均纵波入射角度对应的正演振幅与所述实测数据中所述指定平均纵波入射角度对应的实测振幅之间的差值，以及所述差值的平方；
[0022] 将各个所述平均纵波入射角度分别对应的差值之和的绝对值，以及各个所述平均纵波入射角度分别对应的差值的平方之和分别作为所述指定划分数量对应的第一对比因子和第二对比因子；当所述指定划分数量对应的第一对比因子最小、且对应的第二对比因子最小时，将所述指定划分数量作为所述目标划分数量。
[0023] 优选方案中，确定所述实测共中心点道集数据中的最小纵波入射角度，包括：
[0024] 按照指定纵波入射角度间隔，将所述实测共中心点道集数据划分为多个分道集数据，并对所述分道集数据进行叠加处理；
[0025] 按照纵波入射角度由小到大，依次计算叠加处理后的分道集数据的信噪比，直至所述叠加处理后的分道集数据的信噪比大于或等于指定信噪比阈值，并将信噪比大于或等于指定信噪比阈值的叠加处理后的分道集数据对应的纵波入射角度范围中的最大入射角度作为所述最小纵波入射角度。
[0026] 优选方案中，所述方法还提供有所述目的储层的深度；确定所述实测共中心点道集数据中的最大纵波入射角度，包括：
[0027] 获取所述实测共中心点道集数据对应的最大偏移距；
[0028] 根据所述最大偏移距和所述目的储层的深度，确定所述最大纵波入射角度。
[0029] 优选方案中，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据，包括：
[0030] 从所述实测共中心点道集数据中获取纵波入射角度在所述最小纵波入射角和所述最大纵波入射角之间的目标共中心点道集数据；
[0031] 根据所述目标共中心点道集数据，确定振幅随纵波入射角度变化的目标曲线；其中，所述目标曲线中包括多个目标数据点，所述目标数据点中包括目标振幅值和所述目标振幅值对应的目标纵波入射角度；所述目标振幅值表示所述目标共中心点道集数据中所述目标纵波入射角度对应的地震道数据的采样点的振幅值之和；
[0032] 确定所述目标曲线的拐点，并将所述拐点对应的纵波入射角度作为用于划分共中心点道集数据的目标临界角度；
[0033] 根据所述目标临界角度，将所述目标共中心点道集数据划分为所述目标划分数量的目标分道集数据。
[0034] 优选方案中，所述目标划分数量大于或等于3，确定所述目的储层的含气性，包括：
[0035] 分别对所述目标划分数量的目标分道集数据中第一分道集数据、第二分道集数据和第三分道集数据进行叠加处理；其中，所述第一分道集数据和所述第三分道集数据用于表示所述目标划分数量的目标分道集数据中纵波入射角最小和纵波入射角最大的分道集数据；所述第二分道集数据的纵波入射角度小于所述第三分道集数据的纵波入射角度、且大于所述第一分道集数据的纵波入射角度；
[0036] 基于叠加处理后的第一分道集数据的最大振幅值、叠加处理后的第二分道集数据的最大振幅值和叠加处理后的第三分道集数据的最大振幅值，确定所述目的储层的含气等级。
[0037] 优选方案中，确定所述目的储层的含气等级，包括：
[0038] 将所述叠加处理后的第一分道集数据的最大振幅值减去叠加处理后的第三分道集数据的最大振幅值，得到第一振幅差值，以及将所述叠加处理后的第二分道集数据的最大振幅值减去叠加处理后的第三分道集数据的最大振幅值，得到第二振幅差值；
[0039] 将所述第一振幅差值与所述第二振幅差值相乘，并根据乘积结果，确定所述目的储层的含气等级。
[0040] 优选方案中，确定所述目的储层的含气等级，包括：
[0041] 当所述乘积结果大于或等于指定阈值时，确定所述目的储层的含气等级为第一等级；
[0042] 当所述乘积结果小于所述指定阈值时，确定所述目的储层的含气等级为第二等级。
[0043] 一种确定致密储层含气性的装置，所述装置提供目的储层的实测共中心点道集数据，以及所述目的储层中钻井位置处的测井数据；所述装置包括：正演数据确定模块、目标划分数量确定模块、入射角度确定模块、分道集数据确定模块和含气性确定模块；其中，[0044] 所述正演数据确定模块，用于基于所述测井数据，确定所述目的储层的正演共中心点道集数据；
[0045] 所述目标划分数量确定模块，用于基于所述正演共中心点道集数据和所述实测共中心点道集数据，确定目标划分数量；其中，所述目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量；
[0046] 所述入射角度确定模块，用于确定所述实测共中心点道集数据中的最小纵波入射角度和最大纵波入射角度；
[0047] 所述分道集数据确定模块，用于基于所述最小纵波入射角度和所述最大纵波入射角度，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据；其中，所述目标分道集数据与纵波入射角度范围相对应；
[0048] 所述含气性确定模块，用于根据所述目标分道集数据中的振幅值，确定所述目的储层的含气性。
[0049] 由以上本申请实施例提供的技术方案可见，本申请实施例提供的确定致密储层含气性的方法及装置，可以基于所述测井数据，确定所述目的储层的正演共中心点道集数据；可以基于所述正演共中心点道集数据和所述实测共中心点道集数据，确定目标划分数量；
其中，所述目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量；可以确定所述实测共中心点道集数据中的最小纵波入射角度和最大纵波入射角度；可以基于所述最小纵波入射角度和所述最大纵波入射角度，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据；其中，所述目标分道集数据与纵波入射角度范围相对应；可以根据所述目标分道集数据中的振幅值，确定所述目的储层的含气性。如此，针对实际应用中致密储层的实测地震数据的信噪比较低的情况，采用本申请的方法仍可以准确预测致密储层的含气性。
附图说明
其中，所述目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量；可以确定所述实测共中心点道集数据中的最小纵波入射角度和最大纵波入射角度；可以基于所述最小纵波入射角度和所述最大纵波入射角度，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据；其中，所述目标分道集数据与纵波入射角度范围相对应；可以根据所述目标分道集数据中的振幅值，确定所述目的储层的含气性。如此，针对实际应用中致密储层的实测地震数据的信噪比较低的情况，采用本申请的方法仍可以准确预测致密储层的含气性。
附图说明
[0050] 为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0051] 图1是本申请一种确定致密储层含气性的方法实施例的流程图；
[0052] 图2是本申请实施例中不同划分数量分别对应的正演数据与实测数据的对比示意图；
[0053] 图3是本申请确定致密储层含气性的装置的一种实施例的组成结构示意图；
[0054] 图4是本申请确定致密储层含气性的装置的另一种实施例的组成结构示意图。
具体实施方式
[0055] 本申请实施例提供一种确定致密储层含气性的方法及装置。
[0056] 为了使本技术领域的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。
[0057] 本申请实施例提供了一种确定致密储层含气性的方法。所述确定致密储层含气性的方法提供有目的储层的实测共中心点道集数据、所述目的储层中钻井位置处的测井数据，以及所述目的储层的深度。
[0058] 在本实施方式中，所述目的储层可以是含气性尚未确定的致密储层。所述致密储层例如可以是鄂尔多斯盆地西南部、天环和盆地东部地区中的致密储层。
[0060] 在本实施方式中，所述测井数据具体可以包括：所述目的储层中钻井位置处的测井纵波速度曲线数据、测井横波速度曲线数据和测井密度曲线数据等。所述测井数据可以用于确定所述目的储层中钻井位置处的地震波反射系数。其中，所述地震波反射系数可以包括纵波反射系数和横波反射系数。
[0061] 图1是本申请一种确定致密储层含气性的方法实施例的流程图。如图1所示，所述确定致密储层含气性的方法，包括以下步骤。
[0062] 步骤S101：基于所述测井数据，确定所述目的储层的正演共中心点道集数据。
[0063] 在本实施方式中，基于所述测井数据，确定所述目的储层的正演共中心点道集数据，具体可以包括，可以基于所述实测共中心点道集数据，确定所述目的储层的地震子波数据。可以基于所述测井数据，确定所述目的储层的弹性参数。可以根据预设纵波入射角和所述弹性参数，确定所述目的储层的地震波反射系数。其中，所述预设纵波入射角可以为0度～39度。可以对所述地震子波数据和所述地震波反射系数进行正演处理，得到所述目的储层的正演共中心点道集数据。在本实施方式中，井旁道是指所述实测共中心点道集数据中与所述目的储层中钻井位置距离最近的地震道。所述钻井位置可以是钻井在地表上的位置。
[0064] 在本实施方式中，在地震数据处理领域中，通常地震数据都是地震子波与地层的反射系数进行褶积后的结果。因此，根据已知的地层的反射系数和已知的地震数据可以计算得到地震子波。基于所述实测共中心点道集数据，确定所述目的储层的地震子波数据，具体可以包括，可以获取所述实测共中心点道集数据中的井旁道地震数据和所述目的储层中钻井位置处的地震波反射系数。利用所述实测共中心点道集数据中的井旁道地震数据和所述目的储层中钻井位置处的地震波反射系数，可以确定所述目的储层的地震子波数据。
[0065] 在本实施方式中，基于所述测井数据，确定所述目的储层的弹性参数，具体可以包括，基于所述测井数据，可以获取所述目的储层中钻井位置处的弹性参数。其中，所述弹性参数可以包括：纵波速度、横波速度和密度。基于所述目的储层中钻井位置处的弹性参数，对所述目的储层中非钻井的区域进行反距离加权处理，可以确定所述目的储层中非钻井的区域的弹性参数。根据所述目的储层中钻井位置处的弹性参数和所述目的储层中非钻井的区域的弹性参数，可以确定所述目的储层中的弹性参数。
[0066] 在本实施方式中，根据预设纵波入射角和所述弹性参数，确定所述目的储层的地震波反射系数，具体可以包括，基于佐普利兹(Zoeppritz)和斯奈尔定律，可以建立所述地震波反射系数与所述预设纵波入射角、弹性参数之间的映射关系。例如，对于一个预设纵波入射角和一个弹性参数，可以对应有一个地震波反射系数。根据预设纵波入射角和各个层位的弹性参数的初始值，以及所述建立的地震波反射系数与预设纵波入射角和弹性参数之间的映射关系，可以确定所述地震波反射系数。
[0067] 在本实施方式中，可以采用下述公式对所述地震子波数据和所述地震波反射系数进行正演处理，即褶积运算：
[0068] s'(θ)＝W*R(m,θ)
[0069] 其中，s'(θ)表示所述初始正演地震数据。R(m,θ)表示与弹性参数m和预设纵波入射角θ相关联的地震波反射系数。W表示所述地震子波数据。m表示所述目的储层中的弹性参数。
[0070] 步骤S102：基于所述正演共中心点道集数据和所述实测共中心点道集数据，确定目标划分数量；其中，所述目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量。
[0071] 在本实施方式中，基于所述正演共中心点道集数据和所述实测共中心点道集数据，可以确定目标划分数量，其中，所述目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量，具体可以包括以下步骤：
[0072] (1)针对预设多个划分数量中指定划分数量，按照纵波入射角度，将所述正演共中心点道集数据中指定纵波入射角度范围内的道集数据平均划分为所述指定划分数量的正演分道集数据，以及将所述实测共中心点道集数据中所述指定纵波入射角度范围内的道集数据平均划分为所述指定划分数量的实测分道集数据；
[0073] (2)分别将各个所述正演分道集数据中的采样点的振幅值之和作为一个正演振幅值，得到一组正演数据；所述正演数据中包括所述指定划分数量的正演振幅值和所述正演振幅值对应的平均纵波入射角度；以及分别将所述实测分道集数据中的采样点的振幅值之和作为一个实测振幅值，得到一组实测数据；所述实测数据中包括所述指定划分数量的实测振幅值和所述实测振幅值对应的平均纵波入射角度；其中，所述平均纵波入射角度表示所述正演分道集数据或所述实测分道集数据的纵波入射角度的平均值；
[0074] (3)计算所述正演数据中指定平均纵波入射角度对应的正演振幅与所述实测数据中所述指定平均纵波入射角度对应的实测振幅之间的差值，以及所述差值的平方；
[0075] (4)将各个所述平均纵波入射角度分别对应的差值之和的绝对值，以及各个所述平均纵波入射角度分别对应的差值的平方之和分别作为所述指定划分数量对应的第一对比因子和第二对比因子。这样，针对预设多个划分数量中的每一个划分数量，均可以得到对应的第一对比因子和第二对比因子。当所述指定划分数量对应的第一对比因子最小、且对应的第二对比因子最小时，即正演数据与实测数据的吻合度最好时，将所述指定划分数量作为所述目标划分数量。
[0076] 例如，图2是本申请实施例中不同划分数量分别对应的正演数据与实测数据的对比示意图。图2中(a)、(b)和(c)分别为划分数量分别为3、4、5时对应的正演数据与实测数据的对比示意图。图2中的横坐标和纵坐标分别为纵波入射角度和振幅，纵波入射角度的单位为度。图2中的黑色实心方块和黑色实心圆分别表示正演数据和实测数据。其中，所述最小纵波入射角度为3度，所述最大纵波入射角度为39度。从图2中，可以看出，当划分数量为3时，正演数据与实测数据的吻合度最好。
[0077] 步骤S103：确定所述实测共中心点道集数据中的最小纵波入射角度和最大纵波入射角度。
[0078] 在本实施方式中，确定所述实测共中心点道集数据中的最小纵波入射角度，具体可以包括，可以按照指定纵波入射角度间隔，将所述实测共中心点道集数据划分为多个分道集数据，并对所述分道集数据进行叠加处理。其中，所述指定纵波入射角度间隔的取值范围可以为0.4～0.8度。按照纵波入射角度由小到大，依次计算叠加处理后的分道集数据的信噪比，直至所述叠加处理后的分道集数据的信噪比大于或等于指定信噪比阈值，并将信噪比大于或等于指定信噪比阈值的叠加处理后的分道集数据对应的纵波入射角度范围中的最大入射角度作为所述最小纵波入射角度。其中，所述指定信噪比阈值的取值可以为2.5。例如，所述指定纵波入射角度间隔的取值为0.5度，按照指定纵波入射角度间隔，可以将所述实测共中心点道集数据划分为多个分道集数据，这些分道集数据对应的纵波入射角度范围分别为0～0.5度、0.5～1度、1～1.5度、…、2.5～3度、…。当纵波入射角度范围为2.5～3度的叠加处理后的分道集数据的信噪比大于或等于2.5时，可以将该分道集数据对应的纵波入射角度范围中的最大入射角度，即3度，作为所述最小纵波入射角度。
[0079] 在本实施方式中，确定所述实测共中心点道集数据中的最大纵波入射角度，具体可以包括，可以获取所述实测共中心点道集数据对应的最大偏移距。可以根据所述最大偏移距和所述目的储层的深度，确定所述最大纵波入射角度。具体地，例如，可以采用下述公式计算所述最大纵波入射角度θ：
[0080] θ＝arctan(Xmax/H)
[0081] 其中，Xmax表示所述最大偏移距，H表示所述目的储层的深度。
[0082] 步骤S104：基于所述最小纵波入射角度和所述最大纵波入射角度，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据；其中，所述目标分道集数据与纵波入射角度范围相对应。
[0083] 在本实施方式中，基于所述最小纵波入射角度和所述最大纵波入射角度，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据，具体可以包括以下步骤：
[0084] (1)从所述实测共中心点道集数据中获取纵波入射角度在所述最小纵波入射角和所述最大纵波入射角之间的目标共中心点道集数据；
[0085] (2)根据所述目标共中心点道集数据，确定振幅随纵波入射角度变化的目标曲线；其中，所述目标曲线中包括多个目标数据点，所述目标数据点中包括目标振幅值和所述目标振幅值对应的目标纵波入射角度；所述目标振幅值表示所述目标共中心点道集数据中所述目标纵波入射角度对应的地震道数据的采样点的振幅值之和；
[0086] (3)确定所述目标曲线的拐点，并将所述拐点对应的纵波入射角度作为用于划分共中心点道集数据的目标临界角度；
[0087] (4)根据所述目标临界角度，将所述目标共中心点道集数据划分为所述目标划分数量的目标分道集数据。其中，所述目标分道集数据与纵波入射角度范围相对应。
[0088] 在本实施方式中，所述目标曲线可以用于表征AVO(Amplitude Versus Offset，振幅随偏移距的变化)异常响应特征。其中，所述目标曲线通常情况下所表征的AVO异常响应特征的类型为第一类和第三类AVO异常响应。所述第一类和第三类AVO异常响应分别是正高阻抗差储层和负高阻抗差储层对应的AVO异常响应。
[0089] 在本实施方式中，所述目标曲线的拐点的数量通常可以包括至少两个。这样，在实际应用中，也至少有两个所述目标临界角度。例如，如果所述目标划分数量为3，那么，可以根据两个所述目标临界角度，将所述目标共中心点道集数据划分为3个的目标分道集数据。
[0090] 步骤S105：根据所述目标分道集数据中的振幅值，确定所述目的储层的含气性。
[0091] 在本实施方式中，所述目标划分数量通常可以是大于或等于3。根据所述目标分道集数据中的振幅值，确定所述目的储层的含气性，具体可以包括，可以分别对所述目标划分数量的目标分道集数据中第一分道集数据、第二分道集数据和第三分道集数据进行叠加处理。其中，所述第一分道集数据和所述第三分道集数据用于表示所述目标划分数量的目标分道集数据中纵波入射角最小和纵波入射角最大的分道集数据。所述第二分道集数据的纵波入射角度小于所述第三分道集数据的纵波入射角度、且大于所述第一分道集数据的纵波入射角度。例如，所述目标划分数量为3，所述目标曲线具有两个拐点，这两个拐点分别对应的纵波入射角度为11度和25度。那么，可以将所述目标共中心点道集数据划分为3个的目标分道集数据，按照纵波入射角度由小到大的顺序，分别为近入射角度的分道集数据、中入射角度的分道集数据和远入射角度的分道集数据，分别对应的入射角度范围可以为3～11度、11～25度和25～39度。可以基于叠加处理后的第一分道集数据的最大振幅值、叠加处理后的第二分道集数据的最大振幅值和叠加处理后的第三分道集数据的最大振幅值，确定所述目的储层的含气等级。其中，所述叠加处理具体可以是将分道集数据中的采样点的振幅值求和。
[0092] 在本实施方式中，基于叠加处理后的第一分道集数据的最大振幅值、叠加处理后的第二分道集数据的最大振幅值和叠加处理后的第三分道集数据的最大振幅值，确定所述目的储层的含气等级，具体可以包括，可以将所述叠加处理后的第一分道集数据的最大振幅值减去叠加处理后的第三分道集数据的最大振幅值，得到第一振幅差值，以及将所述叠加处理后的第二分道集数据的最大振幅值减去叠加处理后的第三分道集数据的最大振幅值，得到第二振幅差值。可以将所述第一振幅差值与所述第二振幅差值相乘，以放大不同纵波入射角对应的振幅之间的差异，并根据乘积结果，确定所述目的储层的含气等级。
[0093] 在本实施方式中，根据乘积结果，确定所述目的储层的含气等级，具体可以包括，当所述乘积结果大于或等于指定阈值时，确定所述目的储层的含气等级为第一等级。或者，当所述乘积结果小于所述指定阈值时，确定所述目的储层的含气等级为第二等级。其中，所述第一等级表征的含气性高于所述第二等级表征的含气性。所述第一等级表征所述目标储层的含气性较高，所述第二等级表征所述目的储层的含气性较低。
[0094] 在本实施方式中，所述指定阈值的取值可以为108。在实际应用中，所述指定阈值的取值可以根据实际勘探地区的储层地质情况来设定，此处不作限定。
[0095] 例如，在陇东地区选取33口井位，采用本申请的方法，预测成功率达到75.8％，应用效果较好，适用性广，可以为不同盆地的岩性勘探，尤其是低信噪比地区的岩性勘探提供重要的指导依据。
[0096] 所述确定致密储层含气性的方法实施例，可以基于所述测井数据，确定所述目的储层的正演共中心点道集数据；可以基于所述正演共中心点道集数据和所述实测共中心点道集数据，确定目标划分数量；其中，所述目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量；可以确定所述实测共中心点道集数据中的最小纵波入射角度和最大纵波入射角度；可以基于所述最小纵波入射角度和所述最大纵波入射角度，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据；其中，所述目标分道集数据与纵波入射角度范围相对应；可以根据所述目标分道集数据中的振幅值，确定所述目的储层的含气性。如此，针对实际应用中致密储层的实测地震数据的信噪比较低的情况，采用本申请的方法仍可以准确预测致密储层的含气性。不仅如此，采用本申请方法，基本无人为因素的影响，且计算步骤较少，进一步说明了本申请方法的可行性和实用性。
[0097] 图3是本申请确定致密储层含气性的装置的一种实施例的组成结构示意图。所述确定致密储层含气性的装置提供目的储层的实测共中心点道集数据，以及所述目的储层中钻井位置处的测井数据。如图3所示，所述确定致密储层含气性的装置可以包括：正演数据确定模块100、目标划分数量确定模块200、入射角度确定模块300、分道集数据确定模块400和含气性确定模块500。
[0098] 所述正演数据确定模块100，可以用于基于所述测井数据，确定所述目的储层的正演共中心点道集数据。
[0099] 所述目标划分数量确定模块200，可以用于基于所述正演共中心点道集数据和所述实测共中心点道集数据，确定目标划分数量；其中，所述目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量。
[0100] 所述入射角度确定模块300，可以用于确定所述实测共中心点道集数据中的最小纵波入射角度和最大纵波入射角度。
[0101] 所述分道集数据确定模块400，可以用于基于所述最小纵波入射角度和所述最大纵波入射角度，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据；其中，所述目标分道集数据与纵波入射角度范围相对应。
[0102] 所述含气性确定模块500，可以用于根据所述目标分道集数据中的振幅值，确定所述目的储层的含气性。
[0103] 图4是本申请确定致密储层含气性的装置的另一种实施例的组成结构示意图。如图4所示，所述确定致密储层含气性的装置可以包括：存储器、处理器、以及存储在所述存储器上的计算机程序，所述存储器上存储有目的储层的实测共中心点道集数据，以及所述目的储层中钻井位置处的测井数据，所述计算机程序被所述处理器运行时执行如下步骤：
[0104] 步骤S101：基于所述测井数据，确定所述目的储层的正演共中心点道集数据；
[0105] 步骤S102：基于所述正演共中心点道集数据和所述实测共中心点道集数据，确定目标划分数量；其中，所述目标划分数量用于表征从共中心点道集数据中划分出的分道集数据的数量；
[0106] 步骤S103：确定所述实测共中心点道集数据中的最小纵波入射角度和最大纵波入射角度；
[0107] 步骤S104：基于所述最小纵波入射角度和所述最大纵波入射角度，从所述实测共中心点道集数据中确定所述目标划分数量的目标分道集数据；其中，所述目标分道集数据与纵波入射角度范围相对应；
[0108] 步骤S105：根据所述目标分道集数据中的振幅值，确定所述目的储层的含气性。
[0109] 所述确定致密储层含气性的装置实施例与所述确定致密储层含气性的方法实施例相对应，可以实现确定致密储层含气性的方法实施例的技术方案，并取得方法实施例的技术效果。
[0110] 在20世纪90年代，对于一个技术的改进可以很明显地区分是硬件上的改进(例如，对二极管、晶体管、开关等电路结构的改进)还是软件上的改进(对于方法流程的改进)。然而，随着技术的发展，当今的很多方法流程的改进已经可以视为硬件电路结构的直接改进。设计人员几乎都通过将改进的方法流程编程到硬件电路中来得到相应的硬件电路结构。因此，不能说一个方法流程的改进就不能用硬件实体模块来实现。例如，可编程逻辑器件(Programmable Logic Device,PLD)(例如现场可编程门阵列(Field Programmable Gate Array，FPGA))就是这样一种集成电路，其逻辑功能由用户对器件编程来确定。由设计人员自行编程来把一个数字系统“集成”在一片PLD上，而不需要请芯片制造厂商来设计和制作专用的集成电路芯片。而且，如今，取代手工地制作集成电路芯片，这种编程也多半改用“逻辑编译器(logic compiler)”软件来实现，它与程序开发撰写时所用的软件编译器相类似，而要编译之前的原始代码也得用特定的编程语言来撰写，此称之为硬件描述语言(Hardware Description Language，HDL)，而HDL也并非仅有一种，而是有许多种，如ABEL(Advanced Boolean Expression Language)、AHDL(Altera Hardware Description Language)、Confluence、CUPL(Cornell University Programming Language)、HDCal、JHDL(Java Hardware Description Language)、Lava、Lola、MyHDL、PALASM、RHDL(Ruby Hardware Description Language)等，目前最普遍使用的是VHDL(Very-High-Speed Integrated Circuit Hardware Description Language)与Verilog2。本领域技术人员也应该清楚，只需要将方法流程用上述几种硬件描述语言稍作逻辑编程并编程到集成电路中，就可以很容易得到实现该逻辑方法流程的硬件电路。
[0111] 本领域技术人员也知道，除了以纯计算机可读程序代码方式实现控制器以外，完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至，可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
[0112] 上述实施例阐明的装置、模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。
[0113] 为了描述的方便，描述以上装置时以功能分为各种模块分别描述。当然，在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。
[0114] 通过以上的实施方式的描述可知，本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。该计算机软件产品可以包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。该计算机软件产品可以存储在内存中，内存可能包括计算机可读介质中的非永久性存储器，随机存取存储器(RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括短暂电脑可读媒体(transitory media)，如调制的数据信号和载波。
[0115] 本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。
[0116] 本申请可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
[0117] 本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种电磁波电阻率随钻测井仪 | 2020-05-18 | 283 |
| 随钻声波测井装置 | 2020-05-12 | 828 |
| 一种随钻测井仪器 | 2020-05-13 | 736 |
| 随钻测井仪器外壳及随钻测井仪器 | 2020-05-11 | 1010 |
| 一种多参数随钻测井仪 | 2020-05-13 | 199 |
| 随钻测井仪芯轴 | 2020-05-11 | 603 |
| 随钻核磁共振测井仪探头及随钻核磁共振测井仪 | 2020-05-15 | 117 |
| 随钻核磁共振测井仪探头及随钻核磁共振测井仪 | 2020-05-15 | 87 |
| 随钻偶极声波测井仪器以及声波测井方法 | 2020-05-16 | 506 |
| π随钻测井系统 | 2020-05-15 | 246 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
QQ群二维码
意见反馈
意见反馈