CN104462792B - 一种测井数据岩性层数值归约方法 - Google Patents
一种测井数据岩性层数值归约方法 Download PDFInfo
- Publication number
- CN104462792B CN104462792B CN201410682792.5A CN201410682792A CN104462792B CN 104462792 B CN104462792 B CN 104462792B CN 201410682792 A CN201410682792 A CN 201410682792A CN 104462792 B CN104462792 B CN 104462792B
- Authority
- CN
- China
- Prior art keywords
- log data
- stratum
- data
- extreme point
- value
- 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
Landscapes
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
本发明涉及一种测井数据岩性层数值归约方法，包括以下步骤：1)采用五点二次法对原始测井数据进行平滑处理，得到滤波后的平滑测井数据；2)采用活度分层法对平滑测井数据进行岩性分层，计算平滑测井数据的活度值，并根据活度截止值得到各个地层的分界点，相邻分界点形成一个岩性地层，进而得到岩性地层分层数据；3)根据岩性地层分层数据对原始测井数据进行分段，在每一段原始测井数据上，根据测井数据形态和测井仪器的纵向分辨率，结合岩性分层数据，提取各岩性地层的特征值数据；4)对所有岩性地层的特征值进行数值归约，将特征值相同的相邻岩性层归为一类，生成测井数据岩性层数值归约数据表。本发明可以广泛应用于油气勘探、开发测井储层评价、精细油藏描述等领域。
Description
技术领域
本发明涉及一种石油测井曲线数据预处理方法，特别是关于一种测井数据岩性层数值归约方法。
背景技术
石油测井数据岩性层数值归约(也称测井数据岩性层分层取值)是石油测井储层评价和储量计算不可缺少的重要环节。复杂的数据分析和数据挖掘技术也要求对海量石油测井数据进行归约表示。经过归约处理后的石油测井数据集要比原有的石油测井数据集小很多，而且能够基本保持原有石油测井数据的完整性。这样就使得原来针对海量石油测井数据无法实现的数据处理得以进行，并且产生相同或几乎相同的分析结果。
目前，测井数据岩性层数值归约方法包括两种，一种是人工分层取值，另一种是自动分层取值。其中人工分层取值通常采用的方法为首先确定取值深度，然后从深度曲线上人工读值，读值时根据人对取值深度处相邻深度曲线值大小的综合认识进行的，当人的认识不同时读取的具体数值也会不同，该方法劳动强度大，且读值依赖于人的经验，无法保证取值精度。自动分层取值分为自动分层和自动取值两步，其中，自动分层方法有微分法、活度法、句法分析法和小波分析法等；自动取值方法大多是基于曲线形态进行取值，但是其并没有考虑测井方法本身的技术参数，例如探头的源距、间距等。
发明内容
针对上述问题，本发明的目的是提供一种充分考虑测井方法本身技术参数，效率高、取值精度高的测井数据岩性层数值归约方法。
为实现上述目的，本发明采取以下技术方案：一种测井数据岩性层数值归约方法，包括以下步骤：1)采用五点二次法对原始测井数据进行平滑处理，得到滤波后的平滑测井数据；2)采用活度分层法对平滑测井数据进行岩性分层，计算所述平滑测井数据的活度值，并根据活度截止值得到各个地层的分界点，相邻分界点形成一个岩性地层，进而得到岩性地层分层数据；3)根据岩性地层分层数据对原始测井数据进行分段，在每一段原始测井数据上，根据测井数据形态和测井仪器的纵向分辨率，结合岩性分层数据，提取各岩性地层的特征值数据；4)对所有岩性地层的特征值进行数值归约，将特征值相同的相邻岩性层归为一类，生成测井数据岩性层数值归约数据表。
所述步骤2)中，所述平滑测井数据中第i个滑动平均值的活度值Ei为：
其中，N为窗长L内平滑测井数据的存取样点数，L为给定活度分层窗长，且L∈N+，为滑动平均值yi前、后各N/2范围内所有滑动平均值的平均值，即为：
所述步骤3)中，提取各岩性地层特征值数据的规则为：
①当h<＝d时，取该岩性地层测井数据的极值点作为特征值，其中h表示所述步骤2)中得到的各岩性地层的厚度，d表示测井仪器的纵向分辨率；
②当h>d时，根据该岩性地层测井数据的形态采用不同的取值规则，测井数据形态分为三类：
a、第一类：该岩性地层测井数据的形态特征为其只有一个极值点，此时取该极值点的值作为该岩性地层的特征值；
b、第二类：该岩性地层测井数据的形态特征为其第一个极值点及最后一个极值点同为极大或极小值点，且第一个极值点与最后一个极值点差值大于该层测井数据的最大值与最小值的差值的1/3，此时取第一个极值点与最后一个极值点中较大或较小的值作为该岩性地层的特征值；
c、第三类：该岩性地层测井数据的形态特征为第一个极值点与最后一个极值点同为极大或极小值，或者呈现相反的极性，且各极值点在一个值上下摆动，也就是各个极值点的两两差值都小于该段测井数据的最大值与最小值差值的1/3；则取第一个极值点与最后一个极值点之间所有极值点的平均值作为该岩性地层的特征值。
本发明由于采取以上技术方案，其具有以下优点：1、本发明由于采用五点二次法平滑滤波和活度分层相结合的方法对测井曲线进行自动分层，获得了效率高、满足精度要求的分层取值效果。2、本发明由于针对不同的测井仪器和不同的测井曲线特征采用不同的取值原则，有效的保证了对测井曲线分层取值的精度，能够得到真实的地层测井值。本发明由于将将滤波、分层、取值三个技术环节配套应用，可以获得理想的岩性层测井值，因而本发明可以广泛应用于油气勘探、开发测井储层评价、精细油藏描述等领域。
附图说明
图1是本发明方法流程示意图，
图2是本发明涉及到的三类测井曲线形态，其中横坐标是测井数据取值，纵坐标是测井数据的深度值，竖线为判断极值点差值的基准线。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示，本发明测井数据岩性层数值归约方法包括以下步骤：
1)采用五点二次法对原始测井数据进行平滑处理，得到滤波后的平滑测井数据。
对原始测井数据进行采样，得到测井值xi，其中i＝1,2,...,n，n为总采样点数；则第i个采样点测井值xi的滑动平均值yi为：
2)采用活度分层法对平滑测井数据进行岩性分层，根据活度截止值得到各个地层的分界点，相邻分界点形成一个岩性地层，进而得到岩性地层分层数据。
首先根据平滑测井数据计算活度值，第i个滑动平均值yi的活度值Ei为：
其中，N为窗长L内平滑测井数据的存取样点数，L为给定活度分层窗长，且
L∈N+，为滑动平均值yi前、后各N/2范围内所有滑动平均值的平均值，即为：
根据实际应用中对分层精度的要求，分段设置活度截止值ζ(ζ>＝0)，根据活度截止值对离散测井曲线进行分层，得到各个地层之间的分界点集合O为：
O＝{i|Ei>Ei-1且Ei>Ei+1且Ei>＝ζ,i＝1,.....n} (4)
分界点集合中的每一个分界点是原始测井数据采样点的一个序号i，根据不同采样点序号对平滑测井数据进行分层，得到岩性层分层数据。
3)如图2所示，根据步骤2)得到的岩性地层分层数据对原始测井数据进行分段，在每一段原始测井数据上，根据测井数据形态和测井仪器的纵向分辨率，结合岩性分层数据，提取各岩性地层的特征值数据。提取特征值数据的规则如下：
①当h<＝d时，取该岩性地层测井数据的极值点作为特征值，其中h表示步骤2)中得到的各岩性地层的厚度，d表示测井仪器的纵向分辨率。
②当h>d时，根据该岩性地层测井数据的形态采用不同的取值规则，测井数据形态可以分为三类，具体的：
a、第一类：该岩性地层测井数据的形态特征为其只有一个极值点，此时取该极值点的值作为该岩性地层的特征值(如图2中的①②所示)。
b、第二类：该岩性地层测井数据的形态特征为其第一个极值点及最后一个极值点同为极大或极小值点，且第一个极值点与最后一个极值点差值大于该层测井数据的最大值与最小值的差值的1/3，测井数据总体有一个明显的变化趋势(如图2中的③④所示)，此时取第一个极值点与最后一个极值点中较大或较小的值作为该岩性地层的特征值。
c、第三类：该岩性地层测井数据的形态特征为第一个极值点与最后一个极值点同为极大或极小值，或者两个极值点呈现相反的极性，且该岩性地层各极值点在一个值上下摆动，也就是各个极值点的两两差值都小于该段测井数据的最大值与最小值差值的1/3。则取第一个极值点与最后一个极值点之间所有极值点的平均值作为该岩性地层的特征值(如图2中的⑤⑥⑦⑧所示)。
4)对所有岩性地层的特征值进行数值归约，将特征值相同的相邻岩性层归为一类，生成测井数据岩性层数值归约数据表。
测井数据岩性层数值归约数据表包括岩性层顶界深度、岩性层底界深度和岩性层测井特征值，其中岩性层顶界深度、岩性层底界深度是经过测井数据分层、特征值计算和归约后的分界点对应的深度值。
实施例：
1)采用五点二次法对原始测井数据进行平滑处理，得到滤波后的平滑测井数据；
2)采用活度分层法对平滑测井数据进行岩性分层，根据活度截止值得到各个地层的分界点，相邻分界点形成一个岩性地层，进而得到岩性地层分层数据。
本实施例中，给定活度分层窗长L＝10，利用式(2)和式(3)计算平滑测井数据的活度值，得到活度值为(0.000124，0.000075，0.000062，0.000061，0.000063，0.00008，…)。根据实际应用中对分层精度的要求，分段设置活度截止值ζ，并根据式(4)得到各个地层之间的分界点集合O。
3)根据步骤2)得到的岩性分层数据对原始测井数据进行分段，在每一段原始测井数据上，根据测井数据形态和测井仪器的纵向分辨率，结合岩性分层数据，提取各岩性地层的特征值数据。
本实施例中测井仪器的纵向分辨率d＝0.61m，得到的各岩性地层的特征值数据如下表(如表1所示)，其中depth表示测井深度记录，smo_ri表示滤波后的平滑测井数据，vartual_ri表示完成分层取值后各岩性层的特征值。
表1测井数据分层取值数据表
| depth | smo_ri | virtual-ri |
| 3680 | 8.9688 | 11.336 |
| 3680.1 | 9.8025 | 11.336 |
| 3680.3 | 10.562 | 11.336 |
| 3680.4 | 11.027 | 11.336 |
| 3680.5 | 11.336 | 11.336 |
| 3680.6 | 11.083 | 11.336 |
| 3680.8 | 10.597 | 11.336 |
| 3680.9 | 10.299 | 11.336 |
| 3681 | 10.415 | 11.336 |
| 3681.1 | 10.361 | 11.336 |
| 3681.3 | 10.226 | 11.336 |
| 3681.4 | 10.158 | 11.336 |
| 3681.5 | 10.109 | 11.336 |
| 3681.6 | 9.5757 | 7.4447 |
| 3681.8 | 8.4705 | 7.4447 |
| 3681.9 | 7.4621 | 7.4447 |
| 3682 | 7.4447 | 7.4447 |
| 3682.1 | 8.047 | 7.4447 |
| 3682.3 | 8.5824 | 7.4447 |
| 3682.4 | 8.7534 | 7.4447 |
| 3682.5 | 8.9132 | 7.4447 |
| 3682.6 | 9.2871 | 7.4447 |
| 3682.8 | 10.022 | 11.334 |
| 3682.9 | 10.927 | 11.334 |
| 3683 | 11.631 | 11.334 |
| 3683.1 | 11.907 | 11.334 |
| 3683.3 | 11.846 | 11.334 |
| 3683.4 | 11.751 | 11.334 |
| 3683.5 | 11.477 | 11.334 |
| 3683.6 | 11.199 | 11.334 |
| 3683.8 | 11.088 | 11.334 |
| 3683.9 | 11.295 | 11.334 |
| 3684 | 11.377 | 11.334 |
| 3684.1 | 11.281 | 11.334 |
| 3684.3 | 11.15 | 11.334 |
| 3684.4 | 11.015 | 11.334 |
| 3684.5 | 10.861 | 11.334 |
| 3684.6 | 10.819 | 11.334 |
| 3684.8 | 11.332 | 11.334 |
| 3684.9 | 11.607 | 11.334 |
| 3685 | 11.436 | 11.334 |
| 3685.1 | 11.021 | 11.334 |
| 3685.3 | 10.999 | 11.334 |
| 3685.4 | 10.934 | 11.334 |
| 3685.5 | 10.774 | 11.334 |
| 3685.6 | 10.025 | 7.923 |
| 3685.8 | 8.9471 | 7.923 |
| 3685.9 | 8.1016 | 7.923 |
| 3686 | 7.923 | 7.923 |
| 3686.1 | 7.9853 | 7.923 |
| 3686.3 | 8.1748 | 7.923 |
| 3686.4 | 8.7123 | 12.735 |
| 3686.5 | 10.113 | 12.735 |
| 3686.6 | 11.942 | 12.735 |
| 3686.8 | 12.735 | 12.735 |
| 3686.9 | 11.61 | 12.735 |
| 3687 | 9.8997 | 12.735 |
| 3687.1 | 8.611 | 12.735 |
| 3687.3 | 7.7093 | 12.735 |
| 3687.4 | 6.7585 | 6.1855 |
| 3687.5 | 6.2154 | 6.1855 |
| 3687.6 | 6.1855 | 6.1855 |
| 3687.8 | 6.3612 | 6.1855 |
| 3687.9 | 6.63 | 6.1855 |
| 3688 | 7.3045 | 6.1855 |
| 3688.1 | 8.6267 | 10.317 |
| 3688.3 | 9.8223 | 10.317 |
| 3688.4 | 10.317 | 10.317 |
| 3688.5 | 10.076 | 10.317 |
| 3688.6 | 9.4552 | 10.317 |
| 3688.8 | 8.5322 | 10.317 |
| 3688.9 | 7.6121 | 10.317 |
| 3689 | 6.9543 | 10.317 |
| 3689.1 | 6.797 | 5.6577 |
| 3689.3 | 6.6474 | 5.6577 |
| 3689.4 | 6.2229 | 5.6577 |
| 3689.5 | 5.7131 | 5.6577 |
| 3689.6 | 5.6577 | 5.6577 |
| 3689.8 | 5.7929 | 5.6577 |
| 3689.9 | 5.8502 | 5.6577 |
| 3690 | 5.8555 | 5.6577 |
| 3690.1 | 6.0175 | 5.6577 |
| 3690.3 | 6.4498 | 7.1753 |
| 3690.4 | 7.0182 | 7.1753 |
| 3690.5 | 7.1753 | 7.1753 |
| 3690.6 | 6.8634 | 7.1753 |
| 3690.8 | 6.4006 | 7.1753 |
| 3690.9 | 6.0525 | 7.1753 |
| 3691 | 5.56 | 7.1753 |
| 3691.1 | 5.0802 | 4.8683 |
| 3691.3 | 4.788 | 4.8683 |
| 3691.4 | 4.7137 | 4.8683 |
| 3691.5 | 4.673 | 4.8683 |
| 3691.6 | 4.6471 | 4.8683 |
| 3691.8 | 4.6737 | 4.8683 |
| 3691.9 | 4.7557 | 4.8683 |
| 3692 | 4.8126 | 4.8683 |
| 3692.1 | 4.7825 | 4.8683 |
| 3692.3 | 4.693 | 4.8683 |
| 3692.4 | 4.5636 | 4.8683 |
| 3692.5 | 4.4232 | 4.8683 |
| 3692.6 | 4.339 | 4.8683 |
| 3692.8 | 4.3605 | 4.8683 |
| 3692.9 | 4.4416 | 4.8683 |
| 3693 | 4.5329 | 4.8683 |
| 3693.1 | 4.6541 | 4.8683 |
| 3693.3 | 4.8599 | 4.8683 |
| 3693.4 | 5.1617 | 4.8683 |
| 3693.5 | 5.4145 | 4.8683 |
| 3693.6 | 5.4862 | 4.8683 |
| 3693.8 | 5.3848 | 4.8683 |
| 3693.9 | 5.2895 | 4.8683 |
| 3694 | 5.1928 | 4.8683 |
| 3694.1 | 5.0999 | 4.8683 |
| 3694.3 | 5.0348 | 4.8683 |
| 3694.4 | 5.0566 | 4.8683 |
| 3694.5 | 5.0797 | 4.8683 |
| 3694.6 | 5.0826 | 4.8683 |
| 3694.8 | 5.0411 | 4.8683 |
| 3694.9 | 4.958 | 4.8683 |
| 3695 | 4.8631 | 4.8683 |
| 3695.1 | 4.7916 | 4.8683 |
| 3695.3 | 4.7355 | 4.8683 |
| 3695.4 | 4.7039 | 4.8683 |
| 3695.5 | 4.7126 | 4.8683 |
| 3695.6 | 4.7683 | 4.8683 |
| 3695.8 | 4.8865 | 4.8683 |
| 3695.9 | 5.1473 | 4.8683 |
| 3696 | 5.5875 | 4.8683 |
| 3696.1 | 6.3849 | 8.9609 |
| 3696.3 | 7.4399 | 8.9609 |
| 3696.4 | 8.4399 | 8.9609 |
| 3696.5 | 8.9609 | 8.9609 |
| 3696.6 | 8.892 | 8.9609 |
| 3696.8 | 7.8835 | 8.9609 |
| 3696.9 | 6.4913 | 8.9609 |
| 3697 | 5.7279 | 8.9609 |
| 3697.1 | 6.2548 | 8.4568 |
| 3697.3 | 7.4329 | 8.4568 |
| 3697.4 | 8.2603 | 8.4568 |
| 3697.5 | 8.4568 | 8.4568 |
| 3697.6 | 8.321 | 8.4568 |
| 3697.8 | 8.251 | 8.4568 |
| 3697.9 | 8.0797 | 8.4568 |
| 3698 | 7.7053 | 8.4568 |
| 3698.1 | 7.2963 | 8.4568 |
| 3698.3 | 6.9905 | 8.4568 |
| 3698.4 | 6.9899 | 8.4568 |
| 3698.5 | 7.3783 | 9.0727 |
| 3698.6 | 8.2599 | 9.0727 |
| 3698.8 | 9.0727 | 9.0727 |
| 3698.9 | 9.0304 | 9.0304 |
| 3699 | 8.3402 | 9.0304 |
| 3699.1 | 7.7292 | 9.0304 |
| 3699.3 | 7.4928 | 9.0304 |
| 3699.4 | 7.3134 | 7.0577 |
| 3699.5 | 7.2266 | 7.0577 |
| 3699.6 | 7.1505 | 7.0577 |
| 3699.8 | 7.0577 | 7.0577 |
| 3699.9 | 7.0703 | 7.0577 |
| 3700 | 7.3362 | 7.0577 |
| 3700.1 | 7.355 | 7.0577 |
| 3700.3 | 6.6829 | 6.2653 |
| 3700.4 | 5.9568 | 6.2653 |
| 3700.5 | 5.793 | 6.2653 |
| 3700.6 | 5.9716 | 6.2653 |
| 3700.8 | 6.1615 | 6.2653 |
| 3700.9 | 6.4478 | 6.2653 |
| 3701 | 6.6503 | 6.2653 |
| 3701.1 | 6.5674 | 6.2653 |
| 3701.3 | 6.6676 | 6.2653 |
| 3701.4 | 7.495 | 6.2653 |
| 3701.5 | 9.3761 | 13.473 |
| 3701.6 | 11.639 | 13.473 |
| 3701.8 | 13.266 | 13.473 |
| 3701.9 | 13.473 | 13.473 |
| 3702 | 12.399 | 13.473 |
| 3702.1 | 10.668 | 13.473 |
| 3702.3 | 9.3555 | 8.8203 |
| 3702.4 | 8.8539 | 8.8203 |
| 3702.5 | 8.8203 | 8.8203 |
| 3702.6 | 8.8481 | 8.8203 |
| 3702.8 | 9.0367 | 8.8203 |
| 3702.9 | 9.7264 | 8.8203 |
| 3703 | 11.091 | 12.17 |
| 3703.1 | 12.17 | 12.17 |
| 3703.3 | 12.16 | 12.17 |
| 3703.4 | 11.293 | 12.17 |
| 3703.5 | 10.683 | 12.17 |
| 3703.6 | 10.513 | 12.17 |
| 3703.8 | 9.979 | 12.17 |
| 3703.9 | 8.8085 | 7.0736 |
| 3704 | 7.5731 | 7.0736 |
| 3704.1 | 7.0736 | 7.0736 |
| 3704.3 | 7.1194 | 7.0736 |
| 3704.4 | 7.3565 | 7.0736 |
| 3704.5 | 7.6615 | 7.0736 |
| 3704.6 | 8.179 | 7.0736 |
| 3704.8 | 8.8797 | 7.0736 |
| 3704.9 | 9.3181 | 7.0736 |
| 3705 | 9.7497 | 7.0736 |
| 3705.1 | 10.154 | 7.0736 |
| 3705.3 | 10.572 | 7.0736 |
| 3705.4 | 10.864 | 13.883 |
| 3705.5 | 11.504 | 13.883 |
| 3705.6 | 12.64 | 13.883 |
| 3705.8 | 13.883 | 13.883 |
| 3705.9 | 13.808 | 13.883 |
| 3706 | 12.27 | 13.883 |
| 3706.1 | 10.485 | 8.9726 |
| 3706.3 | 9.4394 | 8.9726 |
| 3706.4 | 8.9726 | 8.9726 |
| 3706.5 | 9.3285 | 8.9726 |
| 3706.6 | 10.372 | 11.029 |
| 3706.8 | 11.029 | 11.029 |
| 3706.9 | 10.333 | 11.029 |
| 3707 | 8.4488 | 11.029 |
| 3707.1 | 6.7774 | 6.8146 |
| 3707.3 | 6.4186 | 6.8146 |
| 3707.4 | 6.9041 | 6.8146 |
| 3707.5 | 7.1725 | 6.8146 |
| 3707.6 | 7.0589 | 6.8146 |
| 3707.8 | 6.9663 | 6.8146 |
| 3707.9 | 6.9302 | 6.8146 |
| 3708 | 6.8293 | 6.8146 |
| 3708.1 | 6.6656 | 6.8146 |
| 3708.3 | 6.4829 | 6.8146 |
| 3708.4 | 6.4109 | 6.8146 |
| 3708.5 | 6.5298 | 6.8146 |
| 3708.6 | 6.6997 | 6.8146 |
| 3708.8 | 6.8423 | 6.8146 |
| 3708.9 | 7.0359 | 6.8146 |
| 3709 | 7.4298 | 6.8146 |
| 3709.1 | 8.0464 | 6.8146 |
| 3709.3 | 8.7052 | 11.533 |
| 3709.4 | 9.4395 | 11.533 |
| 3709.5 | 10.327 | 11.533 |
| 3709.6 | 11.25 | 11.533 |
| 3709.8 | 11.533 | 11.533 |
| 3709.9 | 10.852 | 11.533 |
| 3710 | 9.9446 | 11.533 |
| 3710.1 | 9.2301 | 11.533 |
| 3710.3 | 8.5617 | 7.5842 |
| 3710.4 | 7.7242 | 7.5842 |
| 3710.5 | 7.5842 | 7.5842 |
| 3710.6 | 8.0605 | 7.5842 |
| 3710.8 | 8.7713 | 7.5842 |
| 3710.9 | 9.3404 | 9.9231 |
| 3711 | 9.8557 | 9.9231 |
| 3711.1 | 9.9231 | 9.9231 |
| 3711.3 | 9.1323 | 9.9231 |
| 3711.4 | 7.8218 | 9.9231 |
| 3711.5 | 6.7412 | 9.9231 |
| 3711.6 | 6.3343 | 6.263 |
| 3711.8 | 6.263 | 6.263 |
| 3711.9 | 6.2954 | 6.263 |
| 3712 | 6.4746 | 6.263 |
| 3712.1 | 6.7736 | 6.263 |
| 3712.3 | 7.0827 | 6.263 |
| 3712.4 | 7.3306 | 7.9601 |
| 3712.5 | 7.5937 | 7.9601 |
| 3712.6 | 7.8409 | 7.9601 |
| 3712.8 | 7.9601 | 7.9601 |
| 3712.9 | 7.8721 | 7.9601 |
| 3713 | 7.6131 | 7.9601 |
| 3713.1 | 7.285 | 7.2823 |
| 3713.3 | 7.0219 | 7.2823 |
| 3713.4 | 6.9148 | 6.9486 |
| 3713.5 | 6.9019 | 6.9486 |
| 3713.6 | 6.9151 | 6.9486 |
| 3713.8 | 6.9432 | 6.9486 |
| 3713.9 | 6.997 | 6.9486 |
| 3714 | 6.9746 | 6.9486 |
| 3714.1 | 6.9594 | 6.9486 |
| 3714.3 | 7.0365 | 6.9486 |
| 3714.4 | 7.1794 | 6.9486 |
| 3714.5 | 7.2647 | 6.9486 |
| 3714.6 | 7.3851 | 6.9486 |
| 3714.8 | 7.7587 | 8.0977 |
| 3714.9 | 8.0977 | 8.0977 |
| 3715 | 8.0815 | 8.0977 |
| 3715.1 | 7.6719 | 8.0977 |
| 3715.3 | 7.2521 | 8.0977 |
| 3715.4 | 6.9663 | 8.0977 |
| 3715.5 | 6.78 | 8.0977 |
| 3715.6 | 6.6301 | 6.1097 |
| 3715.8 | 6.5141 | 6.1097 |
| 3715.9 | 6.4708 | 6.1097 |
| 3716 | 6.3894 | 6.1097 |
| 3716.1 | 6.2325 | 6.1097 |
| 3716.3 | 6.1097 | 6.1097 |
| 3716.4 | 6.1317 | 6.1097 |
| 3716.5 | 6.2894 | 6.1097 |
| 3716.6 | 6.5113 | 6.1097 |
| 3716.8 | 6.6904 | 6.1097 |
| 3716.9 | 6.7884 | 6.1097 |
| 3717 | 7.0129 | 7.3793 |
| 3717.1 | 7.3896 | 7.3793 |
| 3717.3 | 7.7161 | 7.3793 |
| 3717.4 | 7.812 | 7.3793 |
| 3717.5 | 7.8217 | 7.3793 |
| 3717.6 | 7.8041 | 7.3793 |
| 3717.8 | 7.7395 | 7.3793 |
| 3717.9 | 7.6411 | 7.3793 |
| 3718 | 7.6007 | 7.3793 |
| 3718.1 | 7.6879 | 7.3793 |
| 3718.3 | 7.8096 | 7.3793 |
| 3718.4 | 7.707 | 7.3793 |
| 3718.5 | 7.3426 | 7.3793 |
| 3718.6 | 6.9336 | 7.3793 |
| 3718.8 | 6.7151 | 7.3793 |
| 3718.9 | 6.7395 | 7.3793 |
| 3719 | 6.9351 | 7.3793 |
| 3719.1 | 7.1054 | 7.3793 |
| 3719.3 | 7.1071 | 7.3793 |
| 3719.4 | 6.9718 | 7.3793 |
| 3719.5 | 6.7367 | 7.3793 |
| 3719.6 | 6.3391 | 6.3391 |
| 3719.8 | 5.8933 | 6.3391 |
| 3719.9 | 5.5271 | 6.3391 |
| 3720 | 5.3497 | 6.3391 |
| 3720.1 | 5.3593 | 5.2389 |
| 3720.3 | 5.496 | 5.2389 |
| 3720.4 | 5.6787 | 5.2389 |
| 3720.5 | 5.7457 | 5.2389 |
| 3720.6 | 5.6311 | 5.2389 |
| 3720.8 | 5.3305 | 5.2389 |
| 3720.9 | 5.049 | 5.2389 |
| 3721 | 4.8718 | 5.2389 |
| 3721.1 | 4.8055 | 5.2389 |
| 3721.3 | 4.883 | 5.2389 |
| 3721.4 | 5.1656 | 5.2389 |
| 3721.5 | 5.6024 | 5.2389 |
| 3721.6 | 5.942 | 5.7464 |
| 3721.8 | 6.0303 | 5.7464 |
| 3721.9 | 5.8762 | 5.7464 |
| 3722 | 5.6562 | 5.7464 |
| 3722.1 | 5.5095 | 5.7464 |
| 3722.3 | 5.5528 | 5.7464 |
| 3722.4 | 5.6825 | 5.7464 |
| 3722.5 | 5.8086 | 5.7464 |
| 3722.6 | 5.8547 | 5.7464 |
| 3722.8 | 5.7796 | 5.7464 |
| 3722.9 | 5.519 | 5.7464 |
| 3723 | 5.2363 | 5.7464 |
| 3723.1 | 5.0255 | 5.7464 |
| 3723.3 | 4.9205 | 4.6352 |
| 3723.4 | 4.8616 | 4.6352 |
| 3723.5 | 4.8288 | 4.6352 |
| 3723.6 | 4.7357 | 4.6352 |
| 3723.8 | 4.6352 | 4.6352 |
| 3723.9 | 4.6553 | 4.6352 |
| 3724 | 4.9173 | 4.6352 |
| 3724.1 | 5.4445 | 6.5492 |
| 3724.3 | 6.1029 | 6.5492 |
| 3724.4 | 6.5492 | 6.5492 |
| 3724.5 | 6.5007 | 6.5492 |
| 3724.6 | 6.0911 | 6.5492 |
| 3724.8 | 5.6635 | 6.5492 |
| 3724.9 | 5.4133 | 6.5492 |
| 3725 | 5.3357 | 6.5492 |
| 3725.1 | 5.4816 | 5.5816 |
| 3725.3 | 5.7791 | 5.5816 |
| 3725.4 | 5.8717 | 5.5816 |
| 3725.5 | 5.655 | 5.5816 |
| 3725.6 | 5.4053 | 5.5816 |
| 3725.8 | 5.4168 | 5.5816 |
| 3725.9 | 5.5141 | 5.5816 |
| 3726 | 5.5593 | 5.5816 |
| 3726.1 | 5.5495 | 5.5816 |
| 3726.3 | 5.5734 | 5.5816 |
| 3726.4 | 5.6647 | 5.5816 |
| 3726.5 | 5.8847 | 5.5816 |
| 3726.6 | 6.1487 | 5.5816 |
| 3726.8 | 6.2579 | 5.5816 |
| 3726.9 | 5.9804 | 5.5816 |
| 3727 | 5.3707 | 5.5816 |
| 3727.1 | 4.7955 | 5.5816 |
4)对所有岩性地层的特征值进行数值归约，将特征值相同的相邻岩性层归为一类，生成测井数据岩性层数值归约数据表，得到测井数据岩性层数值归约数据表如下(如表2所示)，其中up表示岩性层顶界深度，down表示岩性层底界深度，val表示岩性层测井特征值。
表2测井数据岩性层数值归约表示数据表
| up | down | val |
| 3680 | 3681.5 | 11.336 |
| 3681.6 | 3682.6 | 7.4447 |
| 3682.8 | 3685.5 | 11.334 |
| 3685.6 | 3686.3 | 7.923 |
| 3686.4 | 3687.3 | 12.735 |
| 3687.4 | 3688 | 6.1855 |
| 3688.1 | 3689 | 10.317 |
| 3689.1 | 3690.1 | 5.6577 |
| 3690.3 | 3691 | 7.1753 |
| 3691.1 | 3696 | 4.8683 |
| 3696.1 | 3697 | 8.9609 |
| 3697.1 | 3698.4 | 8.4568 |
| 3698.5 | 3698.8 | 9.0727 |
| 3698.9 | 3699.3 | 9.0304 |
| 3699.4 | 3700.1 | 7.0577 |
| 3700.3 | 3701.4 | 6.2653 |
| 3701.5 | 3702.1 | 13.473 |
| 3702.8 | 3702.9 | 8.8203 |
| 3703 | 3703.8 | 12.17 |
| 3703.9 | 3705.3 | 7.0736 |
| 3705.4 | 3706 | 13.883 |
| 3706.1 | 3706.5 | 8.9726 |
| 3706.6 | 3707 | 11.029 |
| 3707.1 | 3709.1 | 6.8146 |
| 3709.3 | 3710.1 | 11.533 |
| 3710.3 | 3710.8 | 7.5842 |
| 3710.9 | 3711.5 | 9.9231 |
| 3711.6 | 3712.3 | 6.263 |
| 3712.4 | 3713 | 7.9601 |
| 3713.1 | 3713.3 | 7.2823 |
| 3713.4 | 3714.6 | 6.9486 |
| 3714.8 | 3715.5 | 8.0977 |
| 3715.6 | 3716.9 | 6.1097 |
| 3717 | 3719.5 | 7.3793 |
| 3719.6 | 3720 | 6.3391 |
| 3720.1 | 3721.5 | 5.2389 |
| 3721.6 | 3723.1 | 5.7464 |
| 3723.3 | 3724 | 4.6352 |
| 3724.1 | 3725 | 6.5492 |
| 3725.1 | 3727.1 | 5.5816 |
上述各实施例仅用于说明本发明，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。
Claims (2)
1.一种测井数据岩性层数值归约方法，包括以下步骤：
1)采用五点二次法对原始测井数据进行平滑处理，得到滤波后的平滑测井数据；
2)采用活度分层法对平滑测井数据进行岩性分层，计算所述平滑测井数据的活度值，并根据活度截止值得到各个地层的分界点，相邻分界点形成一个岩性地层，进而得到岩性地层分层数据；
3)根据岩性地层分层数据对原始测井数据进行分段，在每一段原始测井数据上，根据测井数据形态和测井仪器的纵向分辨率，结合岩性分层数据，提取各岩性地层的特征值数据；
提取各岩性地层特征值数据的规则为：
①当h<＝d时，取该岩性地层测井数据的极值点作为特征值，其中h表示所述步骤2)中得到的各岩性地层的厚度，d表示测井仪器的纵向分辨率；
②当h>d时，根据该岩性地层测井数据的形态采用不同的取值规则，测井数据形态分为三类：
a、第一类：该岩性地层测井数据的形态特征为其只有一个极值点，此时取该极值点的值作为该岩性地层的特征值；
b、第二类：该岩性地层测井数据的形态特征为其第一个极值点及最后一个极值点同为极大或极小值点，且第一个极值点与最后一个极值点差值大于该层测井数据的最大值与最小值的差值的1/3，此时取第一个极值点与最后一个极值点中较大或较小的值作为该岩性地层的特征值；
c、第三类：该岩性地层测井数据的形态特征为第一个极值点与最后一个极值点同为极大或极小值，或者呈现相反的极性，且各极值点在一个值上下摆动，也就是各个极值点的两两差值都小于该段测井数据的最大值与最小值差值的1/3；则取第一个极值点与最后一个极值点之间所有极值点的平均值作为该岩性地层的特征值；
4)对所有岩性地层的特征值进行数值归约，将特征值相同的相邻岩性层归为一类，生成测井数据岩性层数值归约数据表。
2.如权利要求1所述的一种测井数据岩性层数值归约方法，其特征在于：所述步骤2)中，所述平滑测井数据中第i个滑动平均值的活度值Ei为：
其中，N为窗长L内平滑测井数据的存取样点数，L为给定活度分层窗长，且L∈N+，为滑动平均值yk前、后各N/2范围内所有滑动平均值的平均值，即为：
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410682792.5A CN104462792B (zh) | 2014-11-24 | 2014-11-24 | 一种测井数据岩性层数值归约方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410682792.5A CN104462792B (zh) | 2014-11-24 | 2014-11-24 | 一种测井数据岩性层数值归约方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104462792A CN104462792A (zh) | 2015-03-25 |
| CN104462792B true CN104462792B (zh) | 2017-04-12 |
Family
ID=52908821
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410682792.5A Active CN104462792B (zh) | 2014-11-24 | 2014-11-24 | 一种测井数据岩性层数值归约方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104462792B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111980685B (zh) * | 2019-05-22 | 2023-08-22 | 中国石油天然气股份有限公司 | 测井曲线处理方法及装置 |
| CN112230303A (zh) * | 2020-09-28 | 2021-01-15 | 中国地质科学院地球物理地球化学勘查研究所 | 一种地球物理测井数据的处理方法和系统 |
-
2014
- 2014-11-24 CN CN201410682792.5A patent/CN104462792B/zh active Active
Non-Patent Citations (6)
| Title |
|---|
| A Heuristic Search Method for Optimal Zonation of Well Logs;Hui Chuan Chen et al;《Mathematical Geology》;19861231;第18卷(第5期);第489-500页 * |
| 常用数学物理方法在测井解释中的应用;张磊;《中国优秀硕士学位论文全文数据库基础科学辑》;20130115;第2013年卷(第1期);A011-241 * |
| 提高测井曲线分辨率方法;王涛;《化工管理》;20140531;第112页 * |
| 测井曲线自动分层方法回顾与展望;肖波等;《地球物理学进展》;20101031;第25卷(第5期);第1802-1810页 * |
| 测井曲线自动分层解释方法研究;彭智等;《石油仪器》;20110630;第25卷(第3期);第43-44、48页 * |
| 测井自动分层方法研究及其在岩性识别中的应用;韩科龙;《中国优秀硕士学位论文全文数据库基础科学辑》;20110815;第2011年卷(第8期);正文第2章第2.2节、第3章第3.1节、第3.3节 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104462792A (zh) | 2015-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104533400B (zh) | 一种重构测井曲线的方法 | |
| CN104747185B (zh) | 非均质油藏储层综合分类评价方法 | |
| CN112343587A (zh) | 一种特低渗透油藏优势渗流通道识别表征方法 | |
| CN107728227B (zh) | 一种井网成熟区快速判别地下河道砂体方法 | |
| CN111173507B (zh) | 高含水油田剩余油预测方法 | |
| CN107701179B (zh) | 一种基于常规测井资料的页岩气储层的可压性评价方法 | |
| CN106353820A (zh) | 一种基于知识模型的测井解释方法 | |
| CN104133250A (zh) | 膏盐地层地质卡层方法 | |
| CN104564006A (zh) | 一种低渗气井压裂产水能力判断方法 | |
| CN107885893B (zh) | 描述二氧化碳驱储层非均质性的构建方法 | |
| CN108150160A (zh) | 一种地层中欠压实作用超压的求取方法 | |
| CN112145165B (zh) | 一种微裂缝-孔隙型储层动静态渗透率转换方法 | |
| CN104462792B (zh) | 一种测井数据岩性层数值归约方法 | |
| CN113361111B (zh) | 一种低渗透油藏储层特征模型的划分方法 | |
| CN104459790A (zh) | 含油气性盆地有效储层的分析方法和装置 | |
| CN110390154A (zh) | 一种提高复杂断块油气田油藏数值模拟效率的方法 | |
| CN106481315A (zh) | 陆上砂岩油藏单井可采储量快速确定模型及建立方法 | |
| CN106468171A (zh) | 多信息耦合的递进式砂砾岩岩相判别方法 | |
| CN104834934B (zh) | 一种用于识别储层流体的核体俘获法 | |
| CN111155980B (zh) | 一种水流优势通道识别方法及装置 | |
| CN107369101A (zh) | 滚评目标快速经济可动评价方法 | |
| CN107060744A (zh) | 一种测井地质评价系统 | |
| CN108875109B (zh) | 一种异常地层压力的预测方法及系统 | |
| CN112859172A (zh) | 纵横波微测井资料处理方法及装置 | |
| CN105386757A (zh) | 碎屑岩储层孔隙度正常压实趋势线厘定方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Co-patentee after: China University of Petroleum (Beijing) Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC Research Institute Patentee before: China National Offshore Oil Corporation Co-patentee before: China University of Petroleum (Beijing) |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20191209 Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC research institute limited liability company Patentee before: China Offshore Oil Group Co., Ltd. Co-patentee before: China University of Petroleum (Beijing) |