CN104278990B - 测井数据质量恢复方法及装置 - Google Patents
测井数据质量恢复方法及装置 Download PDFInfo
- Publication number
- CN104278990B CN104278990B CN201310274809.9A CN201310274809A CN104278990B CN 104278990 B CN104278990 B CN 104278990B CN 201310274809 A CN201310274809 A CN 201310274809A CN 104278990 B CN104278990 B CN 104278990B
- Authority
- CN
- China
- Prior art keywords
- tension
- depth
- data
- abnormal
- point
- 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
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000002159 abnormal effect Effects 0.000 claims abstract description 127
- 238000011084 recovery Methods 0.000 claims abstract description 41
- 238000005070 sampling Methods 0.000 claims description 118
- 230000008859 change Effects 0.000 claims description 25
- 230000005856 abnormality Effects 0.000 claims description 7
- 238000003384 imaging method Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 abstract description 10
- 230000001133 acceleration Effects 0.000 description 30
- 238000010586 diagram Methods 0.000 description 10
- 238000012952 Resampling Methods 0.000 description 8
- 238000009499 grossing Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 230000003068 static effect Effects 0.000 description 2
- 208000012661 Dyskinesia Diseases 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/26—Storing data down-hole, e.g. in a memory or on a record carrier
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种测井数据质量恢复方法及装置，属于测井技术领域。所述方法包括：获取测井仪器测量得到的测井曲线及拉动所述测井仪器的绳索的张力曲线；根据绳索的张力曲线上相邻采样点的张力数据，获取所述张力曲线发生异常的异常深度段以及发生异常的张力数据变化幅度；根据所述异常深度段以及发生异常的张力数据变化幅度，获取张力曲线发生异常时异常深度段中各个采样点对应的真实深度；对所述测井曲线上与各个采样点对应的真实深度对应的测井数据进行恢复。本发明通过在测井作业中采用普遍常用的张力曲线作为数据恢复依据，不需要测量额外数据即可实现测井数据恢复，简单易行，同时也不需要增加额外设备，避免成本增加。
Description
技术领域
本发明涉及测井技术领域，特别涉及一种测井数据质量恢复方法及装置。
背景技术
测井数据是测井仪器在穿透地层的过程中，通过测量地层的电学、声学、放射性等物理响应得到的地层物理参数，主要用于判断地层的含油气特征。通常情况下，测井仪器在测量测井数据时需要处于匀速运动状态，以确保整个测量系统的误差稳定。
测井曲线用于记录测量过程中的井深与测量得到的测井数据之间对应关系。目前在测井领域，由于钻井越来越深，井孔控制的难度加大，造成井壁的不规则情况加剧，使得测井仪器在运动过程中经常出现非匀速运动状态，从而引起测井曲线的畸变，影响后期的处理和解释。
在现有技术中，一般利用仪器的加速度曲线对仪器在非匀速运动状态下测量得到的测井曲线进行恢复，这种方法精度相对较高。
在实现本发明的过程中，发明人发现现有技术至少存在以下问题：
一是仪器加速度曲线的测量需要增加专用的加速度测量仪器，使得仪器的总长度加长，导致仪器遇卡风险加大，同时也导致成本增加；二是加速度测量仪器操作比较繁琐，影响校正质量。
发明内容
为了解决现有技术的问题，本发明实施例提供了一种测井数据质量恢复方法。所述技术方案如下：
一方面，提供了一种测井数据质量恢复方法，所述方法包括：
获取测井仪器测量得到的测井曲线及拉动所述测井仪器的绳索的张力曲线；
根据绳索的张力曲线上相邻采样点的张力数据，获取所述张力曲线发生异常的异常深度段以及发生异常的张力数据变化幅度；
根据所述异常深度段以及发生异常的张力数据变化幅度，获取张力曲线发生异常时异常深度段中各个采样点对应的真实深度；
对所述测井曲线上与各个采样点对应的真实深度对应的测井数据进行恢复。
另一方面，提供了一种测井数据质量恢复装置，所述装置包括：
曲线数据获取模块，用于获取测井仪器测量得到的测井曲线及拉动所述测井仪器的绳索的张力曲线；
张力异常判断模块，用于根据绳索的张力曲线上相邻采样点的张力数据，获取所述张力曲线发生异常的异常深度段以及发生异常的张力数据变化幅度；
真实深度获取模块，用于根据所述异常深度段以及发生异常的张力数据变化幅度，获取张力曲线发生异常时异常深度段中各个采样点对应的真实深度；
测井数据恢复模块，用于对所述测井曲线上与各个采样点对应的真实深度对应的测井数据进行恢复。
本发明实施例提供的技术方案带来的有益效果是：
本发明通过在测井作业中采用普遍常用的张力曲线作为数据恢复依据，不需要测量额外数据即可实现测井数据恢复，简单易行，同时也不需要增加额外设备，避免成本增加。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的测井数据质量恢复方法流程图；
图2是本发明实施例提供的仪器异常时加速度曲线、张力曲线以及测量曲线的对应关系示意图；
图3是本发明实施例提供的异常深度段对比示意图；
图4是本发明实施例提供的异常测井曲线校正前后对比示意图；
图5是本发明实施例提供的正常仪器测井曲线与异常仪器校正后测井曲线对比示意图；
图6是本发明实施例提供的测井数据质量恢复方法流程图；
图7是本发明实施例提供的张力曲线中异常深度段起始深度点和结束深度点识别示意图；
图8是本发明实施例提供的测井曲线异常深度段的起始深度点A处校正前和未做平滑处理出现跳变台阶的对比示意图；
图9是本发明实施例提供的测井数据质量恢复装置结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
图1是本发明实施例提供的测井数据质量恢复方法流程图。参见图1，该实施例包括：
101、获取测井仪器测量得到的测井曲线及拉动测井仪器的绳索的张力曲线；
张力曲线用于记录测量过程中拉动仪器的绳索的张力数据与井深之间的对应关系，在现有技术中，张力曲线主要用于人为监控仪器的异常，但是不作为测井数据恢复的依据。
如图2所示，加速度曲线、张力曲线及测井曲线有一致性特点。当仪器被卡住时，仪器的加速度为零，拉动仪器的绳索的张力增加，仪器测量同一深度测井数据。仪器解卡后，随着深度的变化，仪器的加速度曲线首先有很大的波动，然后逐渐恢复正常，拉动仪器的绳索的张力逐渐减小，最后恢复至正常，测井曲线上的测井数据开始变动。由于加速度曲线、张力曲线以及测井曲线有一致性特点，因此利用张力曲线可以对测井数据进行恢复。
102、根据绳索的张力曲线上相邻采样点的张力数据，获取该张力曲线发生异常的异常深度段以及发生异常的张力数据变化幅度；
本发明实施例分别获取拉动仪器的绳索在各个采样点的张力数据，并计算相邻采样点张力的差值，根据相邻采样点张力数据的变化量是否大于预设阈值，来确定仪器是否正常。当相邻采样点张力数据的变化量大于预设阈值时，认为仪器异常。
如图3所示，图3(a)中实线为张力曲线，虚线为张力曲线的异常识别曲线，从识别结果可以看出，在张力曲线被识别曲线框住部分出现异常，同时图3(b)中的测井曲线在相同深度段也出现异常。因此本发明实施例通过张力曲线可以准确识别异常深度段的范围。
103、根据该异常深度段以及发生异常的张力数据变化幅度，获取张力曲线发生异常时异常深度段中各个采样点对应的真实深度；
由于张力曲线与测井曲线具有一致性，当张力曲线异常时，测量仪器被卡住，仪器测量的是同一深度的测井数据，测井数据与实际采样点深度不符。因此本发明实施例需要获取异常深度段中各个采样点的真实深度。
104、对测井曲线上与各个采样点对应的真实深度对应的测井数据进行恢复。
如图4所示，图4(a)为测井数据校正恢复前测井曲线，图4(b)为测井数据校正恢复后测井曲线，异常深度段的测井数据校正后消除了因仪器运动异常产生的层位，提高了数据的质量
如图5所示，图5(b)是仪器2在仪器1测量的同一口井、同一深度段的测量数据，仪器2测量时没有出现异常，数据正常。将图5(a)中校正后数据和图5(b)中正常数据相比，二者有很好的一致性，因此本发明实施例能够恢复真实的地质信息，具有较大的应用价值。
本发明通过在测井作业中采用普遍常用的张力曲线作为数据恢复依据，不需要测量额外数据即可实现测井数据恢复，简单易行，同时也不需要增加额外设备，避免成本增加。
可选地，根据绳索的张力曲线上相邻采样点的张力数据，获取所述张力曲线发生异常的异常深度段以及发生异常的张力数据变化幅度包括：
当第i点张力数据变化量的绝对值大于第一预设阈值时，将第i对应深度确定为仪器解卡开始运行深度，异常张力数据变化幅度确定为第i点与相邻采样点的张力数据差值，其中，张力数据变化量的计算公式为(Ten[i]-Ten[i-1])/Ten[i]，变化幅度的计算公式为Ten[i]-Ten[i-1]，Ten[i]表示第i个采样点对应的张力数据，i∈N；
从第i点对应深度沿张力曲线向上寻找，当第j点张力数据变化量的绝对值小于第二预设阈值时，将第j点对应深度确定为异常深度段结束深度，其中，第二预设阈值小于第一预设阈值，张力数据变化量的计算公式为(Ten[j]-Ten[j-1])/Ten[j]，Ten[j]表示第j个采样点对应的张力数据，j∈N；
从第i点对应深度沿张力曲线向下寻找，当第k点张力数据变化量的绝对值小于第二预设阈值时，将第k点对应深度确定为异常深度段起始深度，其中，第二预设阈值小于第一预设阈值，张力数据变化量的计算公式为(Ten[k]-Ten[k-1])/Ten[k]，Ten[k]表示第k个采样点对应的张力数据，k∈N。
可选地，根据所述异常深度段以及发生异常的张力数据变化幅度，获取张力曲线发生异常时异常深度段中各个采样点对应的真实深度包括：
根据异常深度段中各个采样点对应的预设窗长内多个采样点的张力数据，应用公式一，获取所述各个采样点的平均张力，其中，公式一为TenAvg[i]为第i个采样点的平均张力，2*WLen为窗长，i+n表示第i+n个采样点，i∈N，n∈Z；
根据所述各个采样点的平均张力，应用公式二获取所述各个采样点处仪器的加速度，其中，公式二为Az[i]=(Ten[i]-TenAvg[i])/m，Az[i]为异常段中仪器在第i个采样点的加速度，m为仪器质量；
根据所述各个采样点处仪器的加速度、应用公式三、公式四以及公式五，获取所述异常深度段中各个采样点对应的真实深度，其中公式三为Dep[N]-SDep=TotalDep，公式四为公式五为Dep含有m、Time常数，Dep[N]为异常深度段结束深度，SDep为异常深度段起始深度，TotalDep为异常深度段深度，Time为测井数据采样时间间隔，N为异常变化结束时的采样点，Speed[i]表示异常段中第i个采样点处仪器的速度，i,n∈N,且n<i。
可选地，对所述测井曲线上与各个采样点对应的真实深度对应的测井数据进行恢复包括：
根据所述异常深度段中各个采样点对应的真实深度以及公式六，对异常深度段中各个采样点真实深度对应的测井数据重采样插值进行初步恢复，其中公式六为OutData[(Dep[i]-SDep)/DepLevel]=InData[i]，SDep为异常深度段起始深度，DepLevel为采样深度间隔，InData[i]为第i个采样点的测井数据，OutData[i]为第i个采样点的重采样测井数据，i∈N。
可选地，对所述测井曲线上与各个采样点对应的真实深度对应的测井数据进行恢复还包括：在所述初步恢复后，对异常深度段的起始深度点进行平滑处理。
可选地，所述方法应用于常规测井数据质量恢复或成像测井数据质量恢复。
图6是本发明实施例提供的测井数据质量恢复方法流程图。参见图6，该实施例包括：
601、获取测井仪器测量得到的测井曲线及拉动测井仪器的绳索的张力曲线；步骤601与步骤201实现方法相同，在这里不再重复叙述。
602、当第i点张力数据变化量的绝对值大于第一预设阈值时，将第i对应深度确定为仪器解卡开始运行深度，异常张力数据变化幅度确定为第i点与相邻采样点的张力数据差值，其中，张力数据变化量的计算公式为(Ten[i]-Ten[i-1])/Ten[i]，变化幅度的计算公式为Ten[i]-Ten[i-1]，Ten[i]表示第i个采样点对应的张力数据，i∈N；
在本发明实施例中，分别获取拉动仪器的绳索在各个采样点处的张力数据Ten[i]，并计算相邻采样点处张力的差值Ten[i]-Ten[i-1]，根据张力数据的变化量是否大于第一预设阈值来确定张力曲线是否异常，当张力数据(Ten[i]-Ten[i-1])/Ten[i]的绝对值大于第一预设阈值时，认为仪器解卡开始运行，将第i点与相邻采样点张力数据差值Ten[i]-Ten[i-1]确定为异常张力数据变化幅度。
如图7所示，深度B对应的采样点处的张力数据变化量大于第一预设阈值，认为仪器解卡开始运行，将深度点B处采样点i与相邻采样点的张力数据差值确定为张力数据变化幅度。
603、从第i点对应深度沿张力曲线向上寻找，当第j点张力数据变化量的绝对值小于第二预设阈值时，将第j点对应深度确定为异常深度段结束深度，其中，第二预设阈值小于第一预设阈值，张力数据变化量的计算公式为(Ten[j]-Ten[j-1])/Ten[j]，Ten[j]表示第j个采样点对应的张力数据，j∈N；
在本发明实施例中，当异常张力数据变化幅度确定为第i点与相邻采样点的张力数据差值时，从第i点对应深度沿张力曲线向上依次获取采样点处张力数据Ten[j]，并根据相邻采样点的张力数据变化量是否小于第二预设阈值来确定仪器是否恢复正常运行，即当张力数据(Ten[j]-Ten[j-1])/Ten[j]的绝对值小于第二预设阈值时，认为仪器恢复正常运行，将第j点对应深度确定为异常深度段结束深度，其中，第二预设阈值小于第一预设阈值。
如图7所示，从深度B沿张力曲线向上寻找，深度点C对应的采样点处的张力数据变化量小于第二预设阈值，认为仪器恢复正常运行，将深度点B确定为异常深度段结束深度。
604、从第i点对应深度沿张力曲线向下寻找，当第k点张力数据变化量的绝对值小于第二预设阈值时，将第k点对应深度确定为异常深度段起始深度，其中，张力数据变化量的计算公式为(Ten[k]-Ten[k-1])/Ten[k]，Ten[k]表示第k个采样点对应的张力数据，k∈N；
在本发明实施例中，当异常张力数据变化幅度确定为i点与相邻采样点的张力数据差值时，从第i点对应深度沿张力曲线向下依次获取采样点处张力数据Ten[k]，并根据相邻采样点的张力数据变化量是否小于第二预设阈值来确定仪器是否发生异常，即当张力数据(Ten[k]-Ten[k-1])/Ten[k]的绝对值小于第二预设阈值时，认为仪器开始发生异常，将第k点对应深度确定为异常深度段开始深度，其中，第二预设阈值小于第一预设阈值。
如图7所示，从深度B沿张力曲线向下寻找，深度点A对应的采样点处的张力数据变化量小于第二预设阈值，认为仪器开始发生异常，将深度点A确定为异常深度段起始深度。
在本发明实施例中，也可以先执行步骤604，再执行步骤603，本发明不做限定，由开发人员自行安排。
605、根据异常深度段中各个采样点对应的预设窗长内多个采样点的张力数据，应用公式一，获取各个采样点的平均张力，其中，公式一为TenAvg[i]为第i个采样点的平均张力，2*WLen为窗长，i+n表示第i+n个采样点，i∈N，n∈Z；
窗长是为了获取采样点处绳索的平均张力所选取的采样点间距个数。窗长的取值可由开发人员自行决定，本发明实施例不做限定。
本发明实施例为了提高精度，以窗长内多个采样点的平均张力作为第i个采样点处绳索的张力，实际上也可以直接使用第i点处绳索张力。
606、根据各个采样点的平均张力，应用公式二获取各个采样点处仪器的加速度，其中，公式二为Az[i]=(Ten[i]-TenAvg[i])/m，Az[i]为异常段中仪器在第i个采样点的加速度，m为仪器质量；
在本发明实施例中，异常深度段包括仪器遇卡段和仪器解卡后的恢复段，仪器在遇卡段时，仪器被卡住没有运动，加速度为0。仪器在恢复段为加速度运动状态，根据各个采样点的平均张力TenAvg[i]，仪器质量m，应用公式二Az[i]=(Ten[i]-TenAvg[i])/m，可获取各个采样点处仪器的加速度，其中，Az[i]为异常段中仪器在第i个采样点的加速度。
607、根据各个采样点处仪器的加速度、应用公式三、公式四以及公式五，获取异常深度段中各个采样点对应的真实深度，其中公式三为Dep[N]-SDep=TotalDep，公式四为公式五为Dep含有m、Time常数，Dep[N]为异常深度段结束深度，SDep为异常深度段起始深度，TotalDep为异常深度段深度，Time为测井数据采样时间间隔，N为异常变化结束时的采样点，Speed[i]表示异常段中第i个采样点处仪器的速度，i,n∈N,且n<i。
在本发明实施例中，Dep表示深度，含有m、Time常数，Time为测井数据采样时间间隔，Dep[N]为异常深度段结束深度，N为异常变化结束时的采样点，SDep为异常深度段起始深度，TotalDep为异常深度段深度，由于总深度不变，应用公式三Dep[N]-SDep=TotalDep，可获取常数m和Time的总影响。
又Speed[i]表示异常段中第i个采样点处仪器的速度，应用公式四可获取第i个采样点处仪器的速度，其中Az[n]为异常段中仪器在第n个采样点的加速度。
根据获取的常数m和Time的总影响以及异常段中第i个采样点处仪器的速度Speed[i]，应用公式五可获取异常深度段中各个采样点对应的真实深度Dep[i]。
608、根据异常深度段中各个采样点对应的真实深度以及公式六，对异常深度段中各个采样点真实深度对应的测井数据重采样插值进行初步恢复，其中公式六为OutData[(Dep[i]-SDep)/DepLevel]=InData[i]，SDep为异常深度段起始深度，DepLevel为采样深度间隔，InData[i]为第i个采样点的测井数据，OutData[i]为第i个采样点的重采样测井数据，i∈N。
在本发明实施例中，通过对异常深度段测井数据重采样插值进行初步恢复，例如，假设仪器正常时在异常深度段采样点有10个，但仪器异常时，由于仪器在遇卡段不动，重复采集同一深度的测井数据，而仪器解卡后为加速度运动状态，在遇卡段仪器只能采集7个采样点的测井数据，本发明实施例将恢复段采集的7个采样点的测井数据恢复到整个异常段，并且对空白值进行插值，对测井数据进行初步恢复。
在本发明实施例中，对测井曲线上测井数据进行插值可以采用平均值法、拉格朗日插值法以及牛顿插值法等方法。
609、在初步恢复后，对异常深度段的起始深度点进行平滑处理。
由于插值后异常深度段的起始深度点处易出现跳变，为了减少数据突变，本发明实施例需要对异常深度段的起始深度点处作平滑处理。
如图8所示，图8(a)为测井曲线异常深度段的起始深度点A处校正前示意图，图8(b)为测井曲线异常深度段的起始深度点A处插值重采样未平滑处理时出现跳变台阶示意图。
本发明实施例不仅能够对常规测井曲线进行质量恢复处理，还能够推广应用到成像测井等多维测井曲线的质量恢复处理，静止测量为测井仪器在同一深度点重复测量很多次，因此本发明实施例不适用于静止测量的测井数据质量恢复处理。
图9是本发明实施例提供的测井数据质量恢复装置结构示意图。参见图9，该实施例包括：曲线数据获取模块91、张力异常判断模块92、真实深度获取模块93以及测井数据恢复模块94。
曲线数据获取模块91用于获取测井仪器测量得到的测井曲线及拉动该测井仪器的绳索的张力曲线。张力异常判断模块92与曲线数据获取模块91相连，张力异常判断模块92用于根据绳索的张力曲线上相邻采样点的张力数据，获取该张力曲线发生异常的异常深度段以及发生异常的张力数据变化幅度。真实深度获取模块93与张力异常判断模块92相连，真实深度获取模块93用于根据该异常深度段以及发生异常的张力数据变化幅度，获取张力曲线发生异常时异常深度段中各个采样点对应的真实深度。测井数据恢复模块94与真实深度获取模块93相连，测井数据恢复模块94用于对该测井曲线上与各个采样点对应的真实深度对应的测井数据进行恢复。
可选地，张力异常判断模块92包括：变化幅度获取单元921、结束深度获取单元922以及起始深度获取单元923。
变化幅度获取单元921用于当第i点张力数据变化量的绝对值大于第一预设阈值时，将第i对应深度确定为仪器解卡开始运行深度，异常张力数据变化幅度确定为第i点与相邻采样点的张力数据差值，其中，张力数据变化量的计算公式为(Ten[i]-Ten[i-1])/Ten[i]，变化幅度的计算公式为Ten[i]-Ten[i-1]，Ten[i]表示第i个采样点对应的张力数据，i∈N。结束深度获取单元922与变化幅度获取单元921相连，结束深度获取单元922用于从第i点对应深度沿张力曲线向上寻找，当第j点张力数据变化量的绝对值小于第二预设阈值时，将第j点对应深度确定为异常深度段结束深度，其中，第二预设阈值小于第一预设阈值，张力数据变化量的计算公式为(Ten[j]-Ten[j-1])/Ten[j]，Ten[j]表示第j个采样点对应的张力数据，j∈N。起始深度获取单元923与变化幅度获取单元921相连，起始深度获取单元923用于从第i点对应深度沿张力曲线向下寻找，当第k点张力数据变化量的绝对值小于第二预设阈值时，将第k点对应深度确定为异常深度段起始深度，其中，第二预设阈值小于第一预设阈值，张力数据变化量的计算公式为(Ten[k]-Ten[k-1])/Ten[k]，Ten[k]表示第k个采样点对应的张力数据，k∈N。
可选地，真实深度获取模块93包括：平均张力获取单元931、加速度获取单元932以及真实深度获取单元933。
平均张力获取单元931用于根据异常深度段中各个采样点对应的预设窗长内多个采样点的张力数据，应用公式一，获取该各个采样点的平均张力，其中，公式一为TenAvg[i]为第i个采样点的平均张力，2*WLen为窗长，i+n表示第i+n个采样点，i∈N，n∈Z。加速度获取单元932与平均张力获取单元931相连，加速度获取单元932用于根据该各个采样点的平均张力，应用公式二获取该各个采样点处仪器的加速度，其中，公式二为Az[i]=(Ten[i]-TenAvg[i])/m，Az[i]为异常段中仪器在第i个采样点的加速度，m为仪器质量。真实深度获取单元933与加速度获取单元932相连，真实深度获取单元933用于根据该各个采样点处仪器的加速度、应用公式三、公式四以及公式五，获取该异常深度段中各个采样点对应的真实深度，其中公式三为Dep[N]-SDep=TotalDep，公式四为公式五为Dep含有m、Time常数，Dep[N]为异常深度段结束深度，SDep为异常深度段起始深度，TotalDep为异常深度段深度，Time为测井数据采样时间间隔，N为异常变化结束时的采样点，Speed[i]表示异常段中第i个采样点处仪器的速度，i,n∈N,且n<i。
可选地，测井数据恢复模块94包括初步恢复单元941。
初步恢复单元941用于根据该异常深度段中各个采样点对应的真实深度以及公式六，对异常深度段中各个采样点真实深度对应的测井数据重采样插值进行初步恢复，其中公式六为OutData[(Dep[i]-SDep)/DepLevel]=InData[i]，SDep为异常深度段起始深度，DepLevel为采样深度间隔，InData[i]为第i个采样点的测井数据，OutData[i]为第i个采样点的重采样测井数据，i∈N。
可选地，测井数据恢复模块94还包括平滑处理单元942。
平滑处理单元942与上述初步恢复单元941相连，平滑处理单元942用于在该初步恢复后，对异常深度段的起始深度点进行平滑处理。
可选地，该测井数据质量恢复装置可以应用于常规测井数据质量恢复，也可以应用于成像测井数据质量恢复。
本发明通过在测井作业中采用普遍常用的张力曲线作为数据恢复依据，不需要测量额外数据即可实现测井数据恢复，简单易行，同时也不需要增加额外设备，避免成本增加。
需要说明的是：上述实施例提供的测井数据质量恢复装置在恢复测井数据质量时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的测井数据质量恢复的方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种测井数据质量恢复方法，其特征在于，所述方法包括：
获取测井仪器测量得到的测井曲线及拉动所述测井仪器的绳索的张力曲线；
根据绳索的张力曲线上相邻采样点的张力数据，获取所述张力曲线发生异常的异常深度段以及发生异常的张力数据变化幅度；
根据所述异常深度段以及发生异常的张力数据变化幅度，获取张力曲线发生异常时异常深度段中各个采样点对应的真实深度；
对所述测井曲线上与各个采样点对应的真实深度对应的测井数据进行恢复；
所述根据绳索的张力曲线上相邻采样点的张力数据，获取所述张力曲线发生异常的异常深度段以及发生异常的张力数据变化幅度包括：
当第i点张力数据变化量的绝对值大于第一预设阈值时，将第i对应深度确定为仪器解卡开始运行深度，异常张力数据变化幅度确定为第i点与相邻采样点的张力数据差值，其中，张力数据变化量的计算公式为(Ten[i]-Ten[i-1])/Ten[i]，变化幅度的计算公式为Ten[i]-Ten[i-1]，Ten[i]表示第i个采样点对应的张力数据，i∈N；
从第i点对应深度沿张力曲线向上寻找，当第j点张力数据变化量的绝对值小于第二预设阈值时，将第j点对应深度确定为异常深度段结束深度，其中，第二预设阈值小于第一预设阈值，张力数据变化量的计算公式为(Ten[j]-Ten[j-1])/Ten[j]，Ten[j]表示第j个采样点对应的张力数据，j∈N；
从第i点对应深度沿张力曲线向下寻找，当第k点张力数据变化量的绝对值小于第二预设阈值时，将第k点对应深度确定为异常深度段起始深度，其中，第二预设阈值小于第一预设阈值，张力数据变化量的计算公式为 (Ten[k]-Ten[k-1])/Ten[k]，Ten[k]表示第k个采样点对应的张力数据，k∈N；
其中，N为采样总数。
2.如权利要求1所述的方法，其特征在于，所述方法应用于常规测井数据质量恢复或成像测井数据质量恢复。
3.一种测井数据质量恢复装置，其特征在于，所述装置包括：
曲线数据获取模块，用于获取测井仪器测量得到的测井曲线及拉动所述测井仪器的绳索的张力曲线；
张力异常判断模块，用于根据绳索的张力曲线上相邻采样点的张力数据，获取所述张力曲线发生异常的异常深度段以及发生异常的张力数据变化幅度；
真实深度获取模块，用于根据所述异常深度段以及发生异常的张力数据变化幅度，获取张力曲线发生异常时异常深度段中各个采样点对应的真实深度；
测井数据恢复模块，用于对所述测井曲线上与各个采样点对应的真实深度对应的测井数据进行恢复；
所述张力异常判断模块包括：
变化幅度获取单元，用于当第i点张力数据变化量的绝对值大于第一预设阈值时，将第i对应深度确定为仪器解卡开始运行深度，异常张力数据变化幅度确定为第i点与相邻采样点的张力数据差值，其中，张力数据变化量的计算公式为(Ten[i]-Ten[i-1])/Ten[i]，变化幅度的计算公式为Ten[i]-Ten[i-1]，Ten[i]表示第i个采样点对应的张力数据，i∈N；
结束深度获取单元，用于从第i点对应深度沿张力曲线向上寻找，当第j点张力数据变化量的绝对值小于第二预设阈值时，将第j点对应深度确定为异常深度段结束深度，其中，第二预设阈值小于第一预设阈值，张力数据变化量的计算公式为(Ten[j]-Ten[j-1])/Ten[j]，Ten[j]表示第j个采样点对应的张力数据，j∈N；
起始深度获取单元，用于从第i点对应深度沿张力曲线向下寻找，当第k点 张力数据变化量的绝对值小于第二预设阈值时，将第k点对应深度确定为异常深度段起始深度，其中，第二预设阈值小于第一预设阈值，张力数据变化量的计算公式为(Ten[k]-Ten[k-1])/Ten[k]，Ten[k]表示第k个采样点对应的张力数据，k∈N；其中，N为采样总数。
4.如权利要求3所述的装置，其特征在于，所述装置应用于常规测井数据质量恢复或成像测井数据质量恢复。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310274809.9A CN104278990B (zh) | 2013-07-02 | 2013-07-02 | 测井数据质量恢复方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310274809.9A CN104278990B (zh) | 2013-07-02 | 2013-07-02 | 测井数据质量恢复方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104278990A CN104278990A (zh) | 2015-01-14 |
| CN104278990B true CN104278990B (zh) | 2017-06-13 |
Family
ID=52254281
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310274809.9A Active CN104278990B (zh) | 2013-07-02 | 2013-07-02 | 测井数据质量恢复方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104278990B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111026715A (zh) * | 2019-11-14 | 2020-04-17 | 西北工业大学 | 数据处理方法及装置、计算机存储介质及电子设备 |
| CN113052195A (zh) * | 2019-12-27 | 2021-06-29 | 北京国双科技有限公司 | 测井曲线异常值标注方法、装置、计算机设备和存储介质 |
| CN115880482B (zh) * | 2023-02-17 | 2023-05-09 | 中海油田服务股份有限公司 | 一种测井图像遇卡识别与倾角校正方法、装置及计算设备 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1959665A (zh) * | 2005-11-01 | 2007-05-09 | 中国石油大学(北京) | 一种确定等间隔顺序采样数据的突变界面的方法 |
| CN103806899A (zh) * | 2012-11-14 | 2014-05-21 | 中国石油集团长城钻探工程有限公司 | 一种用于确定仪器在井中的深度的方法和装置 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4297879A (en) * | 1979-07-02 | 1981-11-03 | Howells Anthony P | Well logging correlation method and apparatus |
| CN2175826Y (zh) * | 1993-04-14 | 1994-08-31 | 王连生 | 带有数字智能绞车面板的测井深度测量器 |
| CN101519960B (zh) * | 2009-02-23 | 2011-02-09 | 中国石化集团胜利石油管理局钻井工艺研究院 | 一种自适应钻井深度跟踪校正方法 |
| CN102003172B (zh) * | 2010-10-11 | 2013-04-03 | 青岛杰瑞自动化有限公司 | 测井系统中的测井深度校正计算方法 |
| CN201934090U (zh) * | 2010-12-29 | 2011-08-17 | 廊坊开发区新赛浦石油设备有限公司 | 液压测井绞车的智能控制系统 |
| CN202266255U (zh) * | 2011-08-29 | 2012-06-06 | 天津市泰华科技有限公司 | 在石油测井系统中实现深度测量的装置 |
-
2013
- 2013-07-02 CN CN201310274809.9A patent/CN104278990B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1959665A (zh) * | 2005-11-01 | 2007-05-09 | 中国石油大学(北京) | 一种确定等间隔顺序采样数据的突变界面的方法 |
| CN103806899A (zh) * | 2012-11-14 | 2014-05-21 | 中国石油集团长城钻探工程有限公司 | 一种用于确定仪器在井中的深度的方法和装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104278990A (zh) | 2015-01-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106197970B (zh) | 一种基于优化张紧弦模型的桥索监测方法及系统 | |
| CN110781592B (zh) | 一种设备剩余寿命的预测方法及系统 | |
| CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
| WO2017167191A1 (zh) | 地震数据处理方法和装置 | |
| CN105588883B (zh) | 三维岩石力学参数获取方法和系统 | |
| CN107678064B (zh) | 一种声波时差实时提取方法 | |
| CN110702786B (zh) | 基于多尺度奇异吸引子预测误差的梁结构损伤识别方法 | |
| CN104278990B (zh) | 测井数据质量恢复方法及装置 | |
| US20110166842A1 (en) | Layer stripping method | |
| CN113050189B (zh) | 一种测井曲线的重构方法、装置、设备及存储介质 | |
| CN114861737B (zh) | 分布式光纤传感系统的远端扰动特征提取方法及系统 | |
| CN112528227A (zh) | 一种基于数理统计的传感器异常数据识别方法 | |
| CN112578435B (zh) | 岩石超声测试初至拾取方法及系统 | |
| CN112577664A (zh) | 传感器故障检测方法、装置及相关产品 | |
| CN111737868A (zh) | 基于时间序列算法的随钻自然伽马曲线预测方法 | |
| CN117060379A (zh) | 一种基于爬坡特征的深远海风电功率短期预测方法及系统 | |
| CN104832161B (zh) | 基于双尺度相关对比的自动深度校正方法及系统 | |
| CN111950605A (zh) | 表计识别模型的学习方法、装置、设备和表计识别方法 | |
| CN116662800A (zh) | 一种基于自适应注意力机制的滚动轴承故障诊断方法 | |
| CN114969994B (zh) | 基于单测点多向数据融合的rv减速器故障诊断方法及系统 | |
| US10487643B2 (en) | Two-dimensional imaging with multi-stage processing | |
| CN114722860A (zh) | 一种基于多特征融合卷积神经网络的弱磁异常自适应探测方法 | |
| CN110824555B (zh) | 地震能量均衡方法及装置、计算机可读存储介质 | |
| CN118013636B (zh) | 一种砌体结构抗压性能检测设备及检测方法 | |
| CN113822329B (zh) | 一种水电机组主轴摆度信号处理方法及装置 |
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 |