CN105866860B - 一种确定目标区电阻率阈值的方法和装置 - Google Patents
一种确定目标区电阻率阈值的方法和装置 Download PDFInfo
- Publication number
- CN105866860B CN105866860B CN201610397481.3A CN201610397481A CN105866860B CN 105866860 B CN105866860 B CN 105866860B CN 201610397481 A CN201610397481 A CN 201610397481A CN 105866860 B CN105866860 B CN 105866860B
- Authority
- CN
- China
- Prior art keywords
- resistivity
- temperature
- preset standard
- preset
- actual measurement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V13/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Manufacturing & Machinery (AREA)
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本申请实施例提供一种确定目标区电阻率阈值的方法和装置。该方法包括:将获取的实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法计算。比较得到的结果数据的大小;最小的值所对应的预设温度作为第一温度;计算所述第一温度对应的预设标准电阻率与第一预设数量的实测电阻率的差值。获取按序有连续的第三预设数量的差值大于预设标准电阻率所对应的实测电阻率,将所对应的实测电阻率依次与预设标准中间电阻率进行最小二乘法计算。比较得到的结果数据的大小;最小的值所对应的预设温度作为第二温度;基于第一温度对应的预设标准电阻率和第二温度对应的预设标准电阻率确定目标区的电阻率阈值。
Description
技术领域
本发明涉及油气勘探技术领域,尤其涉及一种确定目标区电阻率阈值的方法和装置。
背景技术
地震数据采集是油气勘探的首要工作环节,如何保证地震数据资料野外采集质量是质量监控的关键。随着地震采集技术的发展,对采集接收设备性能要求也越来越高,这不但给野外采集接收设备的现场管理工作,带来越来越大的挑战,而且对如何更加科学的监控检波器各项性能参数,要求也越来越高。其中,电阻率参数与现场检波器排列的埋藏状态等有密切联系。因此,准确确定出检波器的电阻率参数可以确保野外检波器排列的安全,实现对检波器的有效质控。
现有技术中采用的检波器包括常规检波器和高灵敏度检波器,常规检波器不管埋置状况如何,即使在检波器外围温差变化很大的条件下,内部结构决定了电阻率阈值受温度的影响程度较小,阈值是固定的,所以可以对相同类型的检波器选择同一电阻率阈值作为质控标准。高灵敏度检波器的内部结构不同于常规检波器的构造,高灵敏度检波器由于其内部构造易受环境温度影响,电阻率阈值是随外界温度的变化而变化。在野外检波器的埋置条件复杂多变条件下,就导致相同的地表条件和外界温度,检波器埋置水平的不同导致内部温度存在差别,从而不同检波器之间电阻率阈值存在差异。这就要求现场施工期间需要重点关注温度变化导致的电阻率阈值的变化。
但现有的施工条件下检波器内部构造中不含温度记录功能,且相关仪器做检波器测试时,相关仪器也没有温度记录的功能,就不能记录到检波器测试时的温度值。同时野外地表情况复杂多变,施工人员的施工水平,熟练程度和各种外界因素会影响到检波器埋置效果,相同检波器由于埋置状态的不同,相同时间,同样地表温度的条件下,检波器内部的温度值也会有所差异。因此,现有技术不能准确的确定对应温度条件下的电阻率阈值,无法解决高灵敏度检波器的质控问题。
发明内容
本申请的目的是提供一种确定目标区电阻率阈值的方法和装置,可以准确确定出目标区电阻率阈值,解决高灵敏度检波器的质控问题。
本申请提供的确定目标区电阻率阈值的方法和装置是这样实现的:
一种确定目标区电阻率阈值的方法,所述方法包括:
获取目标区按序排列的第一预设数量的实测电阻率,所述实测电阻率与序号相对应;
将所述实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第一结果数据;
比较所述第一结果数据的数值大小,将所述第一结果数据中最小的值所对应的预设温度作为所述目标区的第一温度;
获取所述第一温度对应的预设标准电阻率,计算所述预设标准电阻率与所述第一预设数量的实测电阻率的差值,所述差值与所述实测电阻率的序号相对应;
比较所述差值与所述预设标准电阻率的大小,获取按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的实测电阻率,将所述所对应的实测电阻率依次与所述预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第二结果数据;
比较所述第二结果数据的数值大小,将所述第二结果数据中最小的值所对应的预设温度作为所述目标区的第二温度;
获取所述第二温度对应的预设标准电阻率,基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率确定所述目标区的电阻率阈值。
一种确定目标区电阻率阈值的装置,所述装置包括:
第一数据获取模块,用于获取目标区按序排列的第一预设数量的实测电阻率,所述实测电阻率与序号相对应;
第一计算模块,用于将所述实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第一结果数据;
第一数据处理模块,用于比较所述第一结果数据的数值大小,将所述第一结果数据中最小的值所对应的预设温度作为所述目标区的第一温度;
第二数据获取模块,用于获取所述第一温度对应的预设标准电阻率;
第二计算模块,用于计算所述预设标准电阻率与所述第一预设数量的实测电阻率的差值,所述差值与所述实测电阻率的序号相对应;
第二数据处理模块,用于比较所述差值与所述预设标准电阻率的大小,获取按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的实测电阻率,将所述所对应的实测电阻率依次与所述预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第二结果数据;
第三数据处理模块,用于比较所述第二结果数据的数值大小,将所述第二结果数据中最小的值所对应的预设温度作为所述目标区的第二温度;
第三数据获取模块,用于获取所述第二温度对应的预设标准电阻率;
电阻率阈值确定模块,用于基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率确定所述目标区的电阻率阈值。
本申请通过将获取的实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法的计算。然后,比较得到的结果数据的数值大小;最小的值所对应的预设温度作为所述目标区的第一温度;接着,计算所述第一温度对应的预设标准电阻率与所述第一预设数量的实测电阻率的差值。获取按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的实测电阻率,将所述实测电阻率依次与所述预设标准中间电阻率进行最小二乘法的计算。然后,比较得到的结果数据的数值大小;最小的值所对应的预设温度作为所述目标区的第二温度;最后,基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率可以准确确定所述目标区的电阻率阈值。与现有技术相比,利用本申请提供的技术方案可以在没有现场实测温度的情况下,可以准确确定出目标区不同温度下的电阻率阈值,为后续判断检波器埋藏状态提供数据支持,解决高灵敏度检波器的质控问题。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本申请提供的确定目标区电阻率阈值的方法的一种实施例的流程图;
图2是本申请提供的获取目标区按序排列的第一预设数量的实测电阻率的一种实施例的流程示意图;
图3是利用本申请实施例的技术方案确定的目标区的电阻率阈值与目标区的实测电阻率的一种比较示意图;
图4是本申请提供的确定目标区电阻率阈值的装置的一种实施例中的结构示意图;
图5是本申请提供的所述电阻率阈值确定模块的一种实施例的结构示意图;
图6是本申请提供的所述第一数据获取模块的一种实施例的结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
下面以几个具体的例子详细说明本申请实施例的具体实现。
以下首先介绍本申请一种确定目标区电阻率阈值的方法的一种实施例。图1是本申请提供的确定目标区电阻率阈值的方法的一种实施例的流程图,本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的系统或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。具体的如图1所示,所述方法可以包括:
S110:获取目标区按序排列的第一预设数量的实测电阻率,所述实测电阻率与序号相对应。
本申请实施例中,可以获取目标区按序排列的第一预设数量的实测电阻率,所述实测电阻率与序号相对应。具体的,图2是本申请提供的获取目标区按序排列的第一预设数量的实测电阻率的一种实施例的流程示意图,如图2所示,可以包括:
S111:根据目标区预设检波器排列的检波器排序的序号依次获取相应的实测电阻率。
具体的,在实际应用中,检波器会按序以排列的形式分布,且检波器会标注相应的序号。以排列为单位,按序获取每个排列中检波器的实测电阻率,所述实测电阻率与检波器的序号相对应。
S112:基于所述目标区勘探资料数据从所述实测电阻率中获取有效电阻率。
具体的,所述目标区勘探资料数据可以包括目标区的踏勘资料,排列资料,地表信息等数据。在实际应用中,由于检波器自身和数据传输存储等原因导致的部分实测电阻率为超限值或异常值,为了确保电阻率的有效性和合理性,可以去除相应的超限值或异常值,将去除相应的超限值或异常值之后的实测电阻率作为所述有效电阻率。
具体的,所述异常值可以包括结合目标区勘探资料数据确定出的与大多数电阻率有差异的电阻率。所述超限值可以包括超过预设值的电阻率,所述预设值可以结合目标区勘探资料数据确定。
S113:将所述有效电阻率按照检波器排序的序号的大小依次进行排序,将所述排序后的有效电阻率作为所述第一预设数量的实测电阻率。
具体的,假设某一排列中包括10个检波器,序号为1001的检波器对应电阻率为383.59欧姆/米;序号为1002的检波器对应电阻率为382.59欧姆/米;序号为1003的检波器对应电阻率为383.52欧姆/米;序号为1004的检波器对应电阻率为383.23欧姆/米;序号为1005的检波器对应电阻率为381.98欧姆/米;序号为1006的检波器对应电阻率为320.59欧姆/米;序号为1007的检波器对应电阻率为384.18欧姆/米;序号为1008的检波器对应电阻率为384.53欧姆/米;序号为1009的检波器对应电阻率为382.69欧姆/米;序号为1010的检波器对应电阻率为384.29欧姆/米。根据步骤S112的处理,假设序号为1006的检波器对应电阻率320.59欧姆/米为异常值。相应的,有效电阻率包括序号1001至1005以及序号1007至1010的检波器所对应的电阻率。相应的,按照检波器的序号大小,对所述有效电阻率重新进行排序,可以得到序号1至5依次对应检波器序号1001至1005的电阻率,序号6至9依次对应检波器序号1007至1010的电阻率。
此外,需要说明的是,本申请实施例所述序号的形式并不仅限于上述的形式,实际应用中,还可以包括其他形式,本申请实施例并不以此为限。
S120:将所述实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第一结果数据。
本申请实施例中,在步骤S110获取目标区按序排列的第一预设数量的实测电阻率之后,可以将所述实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第一结果数据,具体的,可以包括:
分别计算所述实测电阻率与第二预设数量的预设温度对应的预设标准中间电阻率的差值;
分别计算所述预设温度所对应的差值的平方和,将所述平方和作为相应的预设温度所对应的第一结果数据。
在实际应用中,高灵敏度检波器在不同温度下对应的电阻率与而不一样,且同一温度下,高灵敏度检波器的电阻率可以包括中间电阻率、最大电阻率、最小电阻率。具体的,本申请实施例中所述第二预设数量的预设温度可以包括结合目标区实际勘探环境确定的多个温度值。所述第二预设数量可以结合实际应用情况预先设置。相应的,所述预设标准中间电阻率可以包括某一温度对应的中间电阻率。
S130:比较所述第一结果数据的数值大小,将所述第一结果数据中最小的值所对应的预设温度作为所述目标区的第一温度。
本申请实施例中,在步骤S120得到与所述预设温度相对应的第一结果数据之后,可以比较所述第一结果数据的数值大小,将所述第一结果数据中最小的值所对应的预设温度作为所述目标区的第一温度。
S140:获取所述第一温度对应的预设标准电阻率,计算所述预设标准电阻率与所述第一预设数量的实测电阻率的差值,所述差值与所述实测电阻率的序号相对应。
本申请实施例中,在步骤S130得到所述目标区的第一温度之后,可以获取所述第一温度对应的预设标准电阻率,计算所述预设标准电阻率与所述第一预设数量的实测电阻率的差值,所述差值与所述实测电阻率的序号相对应。具体的,所述预设标准电阻率可以至少包括下述中的一种:
预设标准中间电阻率、预设标准最大电阻率、预设标准最小电阻率。
相应的,当所述预设标准电阻率包括预设标准中间电阻率、预设标准最大电阻率和预设标准最小电阻率时,可以分别计算预设标准中间电阻率、预设标准最大电阻率、预设标准最小电阻率与所述第一预设数量的实测电阻率的差值。相应的,由于所述第一预设数量的实测电阻率与相应的序号相对应,所述差值与所述实测电阻率的序号也相对应。
S150:比较所述差值与所述预设标准电阻率的大小,获取按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的实测电阻率,将所述所对应的实测电阻率依次与所述预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第二结果数据。
本申请实施例中,所述将所述所对应的实测电阻率依次与所述预设标准中间电阻率进行最小二乘法的计算,得到与预设温度相对应的第二结果数据可以包括:
分别计算所述差值大于所述预设标准电阻率所对应的实测电阻率与第二预设数量的预设温度对应的预设标准中间电阻率的差值;
分别计算所述预设温度对应的差值的平方和,将所述平方和作为相应的预设温度所对应的第二结果数据。
具体的,所述第三预设数量可以结合实际应用情况预先设置。
S160:比较所述第二结果数据的数值大小,将所述第二结果数据中最小的值所对应的预设温度作为所述目标区的第二温度。
具体的,实际应用中,会出现多段按序有连续的第三预设数量的差值大于所述预设标准电阻率的连续实测电阻率,相应的,本申请实施例中所述第二温度可以包括多个温度值。
S170:获取所述第二温度对应的预设标准电阻率,基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率确定所述目标区的电阻率阈值。
本申请实施例中,在步骤S160得到所述目标区的第二温度之后,可以获取所述第二温度对应的预设标准电阻率,然后,基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率确定所述目标区的电阻率阈值。具体的,所述基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率确定所述目标区的电阻率阈值可以包括:
确定所述按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的电阻率的序号;
以目标区预设检波器排列为单位,依次将所述目标区预设检波器排列中所述第一温度对应的预设标准电阻率中与所述所对应的电阻率的序号的电阻率设置为相应的第二温度对应的预设标准电阻率,将所述设置后的电阻率作为所述目标区的电阻率阈值。
具体的,所述电阻率阈值可以包括中间电阻率、最大电阻率以及最小电阻率。
进一步的,在一些实施例中,在步骤S170确定所述目标区的电阻率阈值之后,所述方法还可以包括:
根据所述目标区的电阻率阈值和所述第一预设数量的实测电阻率确定所述目标区的检波器的埋藏情况数据。
具体的,所述埋藏情况数据可以所述第一预设数量的实测电阻率所对应检波器埋藏是否异常的数据。在一个具体的实施例中,如图3所示,图3是利用本申请实施例的技术方案确定的目标区的电阻率阈值与目标区的实测电阻率的一种比较示意图。其中,横坐标为电阻率序号,纵坐标为电阻率,单位为欧姆/米。图中包括目标区最大电阻率曲线310、中间电阻率曲线320、实测电阻率曲线330、最小电阻率曲线340。从图中实测电阻率曲线330与目标区的电阻率阈值(最大电阻率曲线310、中间电阻率曲线320、最小电阻率曲线340)的比较,可以判断出与实测电阻率与目标区实际电阻率之间的差别,相应的,对于差别较大的电阻率,可以找到对应的检波器,进行埋藏情况的检查。
由以上本申请一种确定目标区电阻率阈值的方法的实施例可见,本申请通过将获取的实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法的计算。然后,比较得到的结果数据的数值大小;最小的值所对应的预设温度作为所述目标区的第一温度;接着,计算所述第一温度对应的预设标准电阻率与所述第一预设数量的实测电阻率的差值。获取按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的实测电阻率,将所述实测电阻率依次与所述预设标准中间电阻率进行最小二乘法的计算。然后,比较得到的结果数据的数值大小;最小的值所对应的预设温度作为所述目标区的第二温度;最后,基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率可以准确确定所述目标区的电阻率阈值。与现有技术相比,利用本申请提供的技术方案可以在没有现场实测温度的情况下,可以准确确定出目标区不同温度下的电阻率阈值,为后续判断检波器埋藏状态提供数据支持,解决高灵敏度检波器的质控问题。
本申请另一方面还提供一种确定目标区电阻率阈值的装置,图4是本申请提供的确定目标区电阻率阈值的装置的一种实施例中的结构示意图;如图4所示,所述装置400可以包括:
第一数据获取模块410,可以用于获取目标区按序排列的第一预设数量的实测电阻率,所述实测电阻率与序号相对应;
第一计算模块420,可以用于将所述实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第一结果数据;
第一数据处理模块430,可以用于比较所述第一结果数据的数值大小,将所述第一结果数据中最小的值所对应的预设温度作为所述目标区的第一温度;
第二数据获取模块440,可以用于获取所述第一温度对应的预设标准电阻率;
第二计算模块450,可以用于计算所述预设标准电阻率与所述第一预设数量的实测电阻率的差值,所述差值与所述实测电阻率的序号相对应;
第二数据处理模块460,可以用于比较所述差值与所述预设标准电阻率的大小,获取按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的实测电阻率,将所述所对应的实测电阻率依次与所述预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第二结果数据;
第三数据处理模块470,可以用于比较所述第二结果数据的数值大小,将所述第二结果数据中最小的值所对应的预设温度作为所述目标区的第二温度;
第三数据获取模块480,可以用于获取所述第二温度对应的预设标准电阻率;
电阻率阈值确定模块490,可以用于基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率确定所述目标区的电阻率阈值。
本申请还提供所述电阻率阈值确定模块490的具体实施例,图5是本申请提供的所述电阻率阈值确定模块的一种实施例的结构示意图,如图5所示,所述电阻率阈值确定模块490可以包括:
序号确定单元491,可以用于确定所述按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的电阻率的序号;
数据处理单元492,可以用于以目标区预设检波器排列为单位,依次将所述目标区预设检波器排列中所述第一温度对应的预设标准电阻率中与所述所对应的电阻率的序号的电阻率设置为相应的第二温度对应的预设标准电阻率,将所述设置后的电阻率作为所述目标区的电阻率阈值。
本申请还提供所述第一数据获取模块410的具体实施例,图6是本申请提供的所述第一数据获取模块的一种实施例的结构示意图,如图6所示,所述第一数据获取模块410可以包括:
第一实测电阻率获取单元411,可以用于根据目标区预设检波器排列的检波器排序的序号依次获取相应的实测电阻率;
有效电阻率获取单元412,可以用于基于所述目标区勘探资料数据从所述实测电阻率中获取有效电阻率;
第二实测电阻率获取单元413,可以用于将所述有效电阻率按照检波器排序的序号的大小依次进行排序,将所述排序后的有效电阻率作为所述第一预设数量的实测电阻率。
本申请还提供所述第一计算模块420的具体实施例,具体的,所述第一计算模块420可以包括:
第一计算单元,可以用于分别计算所述实测电阻率与第二预设数量的预设温度对应的预设标准中间电阻率的差值;
第二计算单元,可以用于分别计算所述预设温度所对应的差值的平方和,将所述平方和作为相应的预设温度所对应的第一结果数据。
本申请还提供所述第二数据处理模块460的具体实施例,具体的,所述第二数据处理模块460包括:
第三计算单元,可以用于分别计算所述差值大于所述预设标准电阻率所对应的实测电阻率与第二预设数量的预设温度对应的预设标准中间电阻率的差值;
第四计算单元,可以用于分别计算所述预设温度对应的差值的平方和,将所述平方和作为相应的预设温度所对应的第二结果数据。
在一些实施例中,所述预设标准电阻率可以至少包括下述中的一种:
预设标准中间电阻率、预设标准最大电阻率、预设标准最小电阻率。
在一些实施例中,所述装置400还可以包括:
埋藏情况数据确定模块,可以用于根据所述目标区的电阻率阈值和所述第一预设数量的实测电阻率确定所述目标区的检波器的埋藏情况数据。
由以上本申请一种确定目标区电阻率阈值的方法和装置的实施例可见,本申请通过将获取的实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法的计算。然后,比较得到的结果数据的数值大小;最小的值所对应的预设温度作为所述目标区的第一温度;接着,计算所述第一温度对应的预设标准电阻率与所述第一预设数量的实测电阻率的差值。获取按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的实测电阻率,将所述实测电阻率依次与所述预设标准中间电阻率进行最小二乘法的计算。然后,比较得到的结果数据的数值大小;最小的值所对应的预设温度作为所述目标区的第二温度;最后,基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率可以准确确定所述目标区的电阻率阈值。与现有技术相比,利用本申请提供的技术方案可以在没有现场实测温度的情况下,可以准确确定出目标区不同温度下的电阻率阈值,为后续判断检波器埋藏状态提供数据支持,解决高灵敏度检波器的质控问题。
本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。
Claims (14)
1.一种确定目标区电阻率阈值的方法,其特征在于,所述方法包括:
获取目标区按序排列的第一预设数量的实测电阻率,所述实测电阻率与序号相对应;
将所述实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第一结果数据;
比较所述第一结果数据的数值大小,将所述第一结果数据中最小的值所对应的预设温度作为所述目标区的第一温度;
获取所述第一温度对应的预设标准电阻率,计算所述预设标准电阻率与所述第一预设数量的实测电阻率的差值,所述差值与所述实测电阻率的序号相对应;
比较所述差值与所述预设标准电阻率的大小,获取按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的实测电阻率,将所述所对应的实测电阻率依次与所述预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第二结果数据;
比较所述第二结果数据的数值大小,将所述第二结果数据中最小的值所对应的预设温度作为所述目标区的第二温度;
获取所述第二温度对应的预设标准电阻率,基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率确定所述目标区的电阻率阈值。
2.根据权利要求1所述的方法,其特征在于,所述基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率确定所述目标区的电阻率阈值包括:
确定所述按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的预设标准电阻率的序号;
以目标区预设检波器排列为单位,依次将所述目标区预设检波器排列中所述第一温度对应的预设标准电阻率中与所述所对应的预设标准电阻率的序号的电阻率设置为相应的第二温度对应的预设标准电阻率,将所述设置后的电阻率作为所述目标区的电阻率阈值。
3.根据权利要求1或2任意一项所述的方法,其特征在于,所述获取目标区按序排列的第一预设数量的实测电阻率包括:
根据目标区预设检波器排列的检波器排序的序号依次获取相应的实测电阻率;
基于所述目标区勘探资料数据从所述实测电阻率中获取有效电阻率;
将所述有效电阻率按照检波器排序的序号的大小依次进行排序,将所述排序后的有效电阻率作为所述第一预设数量的实测电阻率。
4.根据权利要求1或2任意一项所述的方法,其特征在于,所述将所述实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第一结果数据包括:
分别计算所述实测电阻率与第二预设数量的预设温度对应的预设标准中间电阻率的差值;
分别计算所述预设温度所对应的差值的平方和,将所述平方和作为相应的预设温度所对应的第一结果数据。
5.根据权利要求1或2任意一项所述的方法,其特征在于,所述将所述所对应的实测电阻率依次与所述预设标准中间电阻率进行最小二乘法的计算,得到与预设温度相对应的第二结果数据包括:
分别计算所述差值大于所述预设标准电阻率所对应的实测电阻率与第二预设数量的预设温度对应的预设标准中间电阻率的差值;
分别计算所述预设温度对应的差值的平方和,将所述平方和作为相应的预设温度所对应的第二结果数据。
6.根据权利要求1或2任意一项所述的方法,其特征在于,所述预设标准电阻率至少包括下述中的一种:
预设标准中间电阻率、预设标准最大电阻率、预设标准最小电阻率。
7.根据权利要求1或2任意一项所述的方法,其特征在于,所述方法还包括:
根据所述目标区的电阻率阈值和所述第一预设数量的实测电阻率确定所述目标区的检波器的埋藏情况数据。
8.一种确定目标区电阻率阈值的装置,其特征在于,所述装置包括:
第一数据获取模块,用于获取目标区按序排列的第一预设数量的实测电阻率,所述实测电阻率与序号相对应;
第一计算模块,用于将所述实测电阻率依次与第二预设数量的预设温度对应的预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第一结果数据;
第一数据处理模块,用于比较所述第一结果数据的数值大小,将所述第一结果数据中最小的值所对应的预设温度作为所述目标区的第一温度;
第二数据获取模块,用于获取所述第一温度对应的预设标准电阻率;
第二计算模块,用于计算所述预设标准电阻率与所述第一预设数量的实测电阻率的差值,所述差值与所述实测电阻率的序号相对应;
第二数据处理模块,用于比较所述差值与所述预设标准电阻率的大小,获取按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的实测电阻率,将所述所对应的实测电阻率依次与所述预设标准中间电阻率进行最小二乘法的计算,得到与所述预设温度相对应的第二结果数据;
第三数据处理模块,用于比较所述第二结果数据的数值大小,将所述第二结果数据中最小的值所对应的预设温度作为所述目标区的第二温度;
第三数据获取模块,用于获取所述第二温度对应的预设标准电阻率;
电阻率阈值确定模块,用于基于所述第一温度对应的预设标准电阻率和所述第二温度对应的预设标准电阻率确定所述目标区的电阻率阈值。
9.根据权利要求8所述的装置,其特征在于,所述电阻率阈值确定模块包括:
序号确定单元,用于确定所述按序有连续的第三预设数量的差值大于所述预设标准电阻率所对应的预设标准电阻率的序号;
数据处理单元,用于以目标区预设检波器排列为单位,依次将所述目标区预设检波器排列中所述第一温度对应的预设标准电阻率中与所述所对应的预设标准电阻率的序号的电阻率设置为相应的第二温度对应的预设标准电阻率,将所述设置后的电阻率作为所述目标区的电阻率阈值。
10.根据权利要求8或9任意一项所述的装置,其特征在于,所述第一数据获取模块包括:
第一实测电阻率获取单元,用于根据目标区预设检波器排列的检波器排序的序号依次获取相应的实测电阻率;
有效电阻率获取单元,用于基于所述目标区勘探资料数据从所述实测电阻率中获取有效电阻率;
第二实测电阻率获取单元,用于将所述有效电阻率按照检波器排序的序号的大小依次进行排序,将所述排序后的有效电阻率作为所述第一预设数量的实测电阻率。
11.根据权利要求8或9任意一项所述的装置,其特征在于,所述第一计算模块包括:
第一计算单元,用于分别计算所述实测电阻率与第二预设数量的预设温度对应的预设标准中间电阻率的差值;
第二计算单元,用于分别计算所述预设温度所对应的差值的平方和,将所述平方和作为相应的预设温度所对应的第一结果数据。
12.根据权利要求8或9任意一项所述的装置,其特征在于,所述第二数据处理模块包括:
第三计算单元,用于分别计算所述差值大于所述预设标准电阻率所对应的实测电阻率与第二预设数量的预设温度对应的预设标准中间电阻率的差值;
第四计算单元,用于分别计算所述预设温度对应的差值的平方和,将所述平方和作为相应的预设温度所对应的第二结果数据。
13.根据权利要求8或9任意一项所述的装置,其特征在于,所述预设标准电阻率至少包括下述中的一种:
预设标准中间电阻率、预设标准最大电阻率、预设标准最小电阻率。
14.根据权利要求8或9任意一项所述的装置,其特征在于,所述装置还包括:
埋藏情况数据确定模块,用于根据所述目标区的电阻率阈值和所述第一预设数量的实测电阻率确定所述目标区的检波器的埋藏情况数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610397481.3A CN105866860B (zh) | 2016-06-07 | 2016-06-07 | 一种确定目标区电阻率阈值的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610397481.3A CN105866860B (zh) | 2016-06-07 | 2016-06-07 | 一种确定目标区电阻率阈值的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105866860A CN105866860A (zh) | 2016-08-17 |
CN105866860B true CN105866860B (zh) | 2018-03-13 |
Family
ID=56677160
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610397481.3A Active CN105866860B (zh) | 2016-06-07 | 2016-06-07 | 一种确定目标区电阻率阈值的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105866860B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6256587B1 (en) * | 1998-11-17 | 2001-07-03 | Baker Hughes, Inc. | Method for correcting well log data for effects of changes in instrument velocity (cable yo-yo) |
CN202614956U (zh) * | 2012-05-11 | 2012-12-19 | 威海双丰物探设备股份有限公司 | 一种检波器参数测试装置 |
WO2014130342A1 (en) * | 2013-02-20 | 2014-08-28 | Apache Corporation | Methods for determining well log attributes for formation characterization |
CN104459804A (zh) * | 2014-12-18 | 2015-03-25 | 上海艾都能源科技有限公司 | 野外异常区边缘圈定的快速电磁勘探方法 |
CN104656145A (zh) * | 2014-11-20 | 2015-05-27 | 河南理工大学 | 基于数字电极技术的电法测量系统 |
CN204405873U (zh) * | 2015-02-14 | 2015-06-17 | 合肥国为电子有限公司 | 兼具电磁和地震数据采集功能的地球物理勘探仪器 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015012870A1 (en) * | 2013-07-26 | 2015-01-29 | Halliburton Energy Services Inc. | Computer program for wellbore resistivity logging calibration |
AU2014388388B2 (en) * | 2014-03-25 | 2017-11-23 | Halliburton Energy Services, Inc. | Permanent EM monitoring systems using capacitively coupled source electrodes |
US10459107B2 (en) * | 2014-11-13 | 2019-10-29 | Halliburton Energy Services, Inc. | Well monitoring with autonomous robotic diver |
-
2016
- 2016-06-07 CN CN201610397481.3A patent/CN105866860B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6256587B1 (en) * | 1998-11-17 | 2001-07-03 | Baker Hughes, Inc. | Method for correcting well log data for effects of changes in instrument velocity (cable yo-yo) |
CN202614956U (zh) * | 2012-05-11 | 2012-12-19 | 威海双丰物探设备股份有限公司 | 一种检波器参数测试装置 |
WO2014130342A1 (en) * | 2013-02-20 | 2014-08-28 | Apache Corporation | Methods for determining well log attributes for formation characterization |
CN104656145A (zh) * | 2014-11-20 | 2015-05-27 | 河南理工大学 | 基于数字电极技术的电法测量系统 |
CN104459804A (zh) * | 2014-12-18 | 2015-03-25 | 上海艾都能源科技有限公司 | 野外异常区边缘圈定的快速电磁勘探方法 |
CN204405873U (zh) * | 2015-02-14 | 2015-06-17 | 合肥国为电子有限公司 | 兼具电磁和地震数据采集功能的地球物理勘探仪器 |
Non-Patent Citations (1)
Title |
---|
基于检波器的短周期地震计的改进;张晓鹏 等;《地球物理学进展》;20141231;第29卷(第5期);第2466-2471页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105866860A (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Woessner et al. | Assessing the quality of earthquake catalogues: Estimating the magnitude of completeness and its uncertainty | |
Chousianitis et al. | Predictive model of Arias intensity and Newmark displacement for regional scale evaluation of earthquake-induced landslide hazard in Greece | |
Dunn et al. | Pairwise homogeneity assessment of HadISD | |
WO2016141598A1 (zh) | 共姿态道集地表一致性振幅补偿方法及装置 | |
CN111914457B (zh) | 一种输电线路塔基边坡稳定性判别方法、装置和存储介质 | |
CN104298923B (zh) | 漏洞类型识别方法以及装置 | |
CN103258128A (zh) | 一种地磁场空间延拓算法的评估方法 | |
CN106706715A (zh) | 基于三维高密度电阻率法的污染土检测方法 | |
Twiddle et al. | Pollen productivity estimates for a pine woodland in eastern Scotland: The influence of sampling design and vegetation patterning | |
CN115660407A (zh) | 污染场地风险等级预测模型的训练方法和训练装置 | |
CN103678778A (zh) | 一种放射性物化探信息综合方法 | |
CN106706119B (zh) | 一种基于信号频域特征的振动源识别方法及系统 | |
Jiang et al. | 3D seismic geometry quality control and corrections by applying machine learning | |
CN105866860B (zh) | 一种确定目标区电阻率阈值的方法和装置 | |
CN104296845B (zh) | 基于最佳底面压强测量点的粮仓储粮重量检测方法与装置 | |
Ommi et al. | Analyses of seismicity parameters of the August 11th, 2012, Ahar-Varzaghan earthquakes in north-western Iran | |
CN106285661A (zh) | 一种基于判定指数的低阻油层识别方法及装置 | |
CN104819382A (zh) | 一种用于光纤预警系统的自适应恒虚警率振源检测方法 | |
CN116070907A (zh) | 一种基于层次分析的岩溶塌陷易发性评估方法及系统 | |
CN111221043B (zh) | 一种被动源电场法的工作参数优化方法 | |
CN109001799A (zh) | 一种自动识别异常道的方法及系统 | |
CN205139068U (zh) | 基于高密度电阻率法的污染土检测测线布置结构 | |
CN105510986B (zh) | 一种确定海绿石砂岩中海绿石含量的方法和装置 | |
CN107194187B (zh) | 顾及空间组成与构型的地表覆盖样本量计算方法 | |
RU2509890C1 (ru) | Способ контроля и определения координат опасного состояния массива горных пород при подземных горных работах |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |