CN107218033A - 一种识别地层产状的方法和地层参数的反演方法 - Google Patents
一种识别地层产状的方法和地层参数的反演方法 Download PDFInfo
- Publication number
- CN107218033A CN107218033A CN201710337945.6A CN201710337945A CN107218033A CN 107218033 A CN107218033 A CN 107218033A CN 201710337945 A CN201710337945 A CN 201710337945A CN 107218033 A CN107218033 A CN 107218033A
- Authority
- CN
- China
- Prior art keywords
- log
- depth
- bed boundary
- composite log
- composite
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 55
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 34
- 239000002131 composite material Substances 0.000 claims abstract description 83
- 238000003384 imaging method Methods 0.000 claims abstract description 33
- 230000000052 comparative effect Effects 0.000 claims abstract description 17
- 230000009977 dual effect Effects 0.000 claims description 17
- 230000004044 response Effects 0.000 claims description 12
- 230000000694 effects Effects 0.000 claims description 11
- 230000000875 corresponding effect Effects 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 6
- 238000013329 compounding Methods 0.000 claims description 6
- 239000002002 slurry Substances 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 5
- 230000002195 synergetic effect Effects 0.000 claims description 4
- 230000001174 ascending effect Effects 0.000 claims description 3
- 230000009545 invasion Effects 0.000 abstract description 9
- 238000004364 calculation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000003786 synthesis reaction Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000011835 investigation Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012887 quadratic function Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000012952 Resampling Methods 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000004087 circulation Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000013519 translation Methods 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
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 And Detection Of Objects (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种识别地层产状的方法和地层参数的反演方法，包括如下步骤：获取电成像仪器每个极板上的合成测井曲线；对每条合成测井曲线进行预分层，得到每条合成测井曲线的预分层界面位置；依次选取包含预分层界面的深度段作为参考段，对其它极板上的合成测井曲线进行相关对比，得到各个合成测井曲线的有效层界面；对所有极板上的合成测井曲线的有效层界面按深度进行排序，获得地层界面并计算产状信息。进一步地，在地层分层的基础上，结合双侧向资料建立正演模型，通过反演得到泥浆侵入深度、侵入带与原状地层电阻率。本发明提取的地层界面与测井图像吻合较好，产状信息更准确，反演出的地层参数可以定量描述泥浆的侵入特性。
Description
技术领域
本发明涉及油气资源地球物理勘查技术领域，特别涉及一种识别地层产状的方法和地层参数的反演方法。
背景技术
为了有效探测油气层，在油气田勘探开发的过程中，测井资料数字处理成为广大测井解释人员的主要任务之一，现有的商业软件中的地层产状计算都是基于固定窗长的相关对比法，这种方法的优点是计算速度快，但是缺点是无法识别出层界面，计算出的地层产状是窗长范围内所有层界面倾角的综合均值，这种产状信息可以反映出窗长范围内层界面的大致趋势，但是无法反映出每个层界面的精确倾角与倾向，无法应用于精细沉积分析。为此，测井解释人员通常的做法是选择手工方式，通过软件的人机交互功能逐一拾取层界面并计算地层产状，这种做法需要付出很大的工作量，且效率不高。
此外，目前电成像测井与双侧向测井是油气勘探的两项重要技术手段，电成像测井的探测深度约为2厘米，因此，电成像测井仅能获取井壁很浅的地层视电阻率信息，对井壁地层特征进行评价；双侧向测井的深探测深度和浅探测深度分别约为110厘米和20厘米，因此，双侧向测井能够获取距离井壁20厘米之外的视电阻率信息，对渗透性较好的地层进行判别。
但是，通过电成像测井与双侧向测井，均无法获取介于探测深度2至20厘米之间区域的地层电阻率信息，因而很难对该区域内的泥浆侵入特性进行精细评价。因此考虑如何通过电成像与双侧向测井信息来进行联合反演，进而精确获取全井段地层的泥浆侵入信息，对地层的渗透性评价与储层判别具有重要意义。
发明内容
为了解决上述技术问题，本发明提供了一种识别地层产状的方法和地层参数的反演方法，能够精确获取全井段地层的泥浆侵入信息。
为了达到本发明目的，本发明实施例的技术方案是这样实现的：
本发明实施例提供了一种识别地层产状的方法，包括：
获取电成像仪器每个极板上的合成测井曲线；
对每条合成测井曲线进行预分层，得到每条合成测井曲线的预分层界面位置；
依次选取包含预分层界面位置的深度段作为参考段，对其它极板上的合成测井曲线进行相关对比，得到各个合成测井曲线的有效层界面；
对所有极板上的合成测井曲线的有效层界面按深度进行排序，获得地层界面并计算地层产状信息。
进一步地，所述获取电成像仪器每块极板上的合成测井曲线，具体包括：
对所述每个极板内若干个纽扣电极的测井曲线进行深度对齐；
对所述每个极板内每个采样点上若干个钮扣电极的测井数据求取平均值，将所述平均值作为所述采样点的测量值，从而得到所述每个极板上的合成测井曲线。
进一步地，所述获取电成像仪器每块极板上的合成测井曲线之前，所述方法还包括：
对所述每个极板内每个采样点上所有钮扣电极的测井数据进行升序或降序排序，按照预设的第一剔除比例值剔除较大测井数据，以及第二剔除比例值剔除较小测井数据。
进一步地，所述对每条合成测井曲线进行预分层，得到每条合成测井曲线的预分层界面位置的步骤，具体包括：
对每条所述合成测井曲线进行平滑滤波处理；
获取每条所述合成测井曲线对应的活度曲线和对比度曲线；
求取每条活度曲线的极大值点，并当极大值点在预设的窗长阈值范围内与其余采样点的对比度值大于预设的对比度阈值时，将所述极大值点的深度作为所述合成测井曲线的预分层界面位置。
进一步地，所述依次选取包含预分层界面位置的深度段作为参考段，对其它的合成测井曲线进行相关对比，得到各个合成测井曲线的有效层界面的步骤，具体包括：
对所述每个极板上的合成测井曲线执行如下步骤：
选取每个极板上的所述合成测井曲线作为参考曲线，将其它极板上的所述合成测井曲线作为对比曲线；
在参考曲线上，按照深度顺序依次选取三个所述预分层界面的深度段数据作为参考段，第一个深度点与第三个深度点之间的深度差为本次对比窗长；
在各对比曲线上，使用相关对比法在预设的探索长度内依次计算出最大相关系数时的对比段中点深度；
计算满足最大相关系数大于预设的相关系数阈值的对比曲线条数，检测所述对比曲线条数是否大于或等于预设的对比曲线条数阈值，如果大于或等于预设的对比曲线条数阈值，将参考曲线中的所述预分层界面作为有效层界面，将所述获得最大相关系数时的对比段中点深度，作为与参考曲线中目标层界面相同层的层界面深度；
对参考曲线上的所有预分层界面进行相关对比，得到该极板上的合成测井曲线的所有有效层界面。
进一步地，所述在各对比曲线上，使用相关对比法在预设的探索长度内依次计算出最大相关系数时的对比段中点深度的步骤，具体包括：
在各对比曲线上，在预设的探索长度范围内依次选取与本次对比窗长长度相等的对比段，计算各对比段与参考段的合成测井曲线的相关系数，记录每次探索长度范围内的最大相关系数和获得最大相关系数时的对比段中点深度。
进一步地：所述方法还包括：
当相邻两个所述有效层界面的深度差小于预设的深度阈值时，或者当所述相邻两个所述有效层界面存在交叉时，在所述获得的地层界面中，去除所述相邻两个有效层界面中相关系数小的有效层界面。
本发明实施例还提供了一种地层参数的反演方法，包括：
获取电成像仪器每个极板上的合成测井曲线；
对每条合成测井曲线进行预分层，得到每条合成测井曲线的预分层界面位置；
依次选取包含预分层界面位置的深度段作为参考段，对其它极板上的合成测井曲线进行相关对比，得到各个合成测井曲线的有效层界面；
对所有极板上的合成测井曲线的有效层界面按深度进行排序，获得地层界面；
基于电成像测井数据，提取各地层界面对应的地层倾角、地层倾向，提取冲洗带视电阻率信息；
根据获得的地层界面、提取的地层倾角与冲洗带视电阻率信息，创建初始正演模型，进行地层仿真，将仿真获取的测井响应值与实际测井数据进行对比，判断是否满足终止迭代条件，如果不满足终止迭代条件，修正地层参数，重新进行仿真；如果满足终止迭代条件，提取全井段地层的泥浆侵入信息。
进一步地，在所述根据获得的地层界面、提取的地层倾角与冲洗带视电阻率信息，创建初始正演模型，进行地层仿真的步骤中，当所述地层倾角小于预设的角度阈值时，通过二维模型进行仿真；当所述地层倾角大于预设的角度阈值时，通过三维模型进行仿真。
进一步地，所述终止迭代条件包括以下任意一个：
将所述双侧向测井响应值与双侧向实测值进行比较，响应值与实测值之间的差值小于预先设置的第一差值阈值；
迭代步数大于预先设置的最大迭代步数；
将相邻两次仿真得到的双侧向测井响应值与双侧向实测值的差值进行比较，所述差值的变化小于预先设置的第二差值阈值充电器测量自身的振动数据。
本发明的技术方案，具有如下有益效果：
由于电成像图像上的层界面本质是不同岩性的地层电阻率不一致引起的图像颜色跳变，在电扣曲线上表现为极值点或者突变值点，本发明提供的识别地层产状的方法是在检测这些极值点或突变点的基础上，识别出地层界面，然后计算出层界面产状信息，所以提取的层界面与图像特征吻合较好，计算出的地层产状与真实状态很接近，通过本发明计算地层产状快速高效，准确度高，基本可以取代目前的人工拾取方法；
进一步地，本发明提供的地层参数的反演方法，通过电成像与双侧向测井信息进行联合反演，可以计算出泥浆侵入深度、侵入带地层电阻率与原状地层电阻率，进而获取到全井段地层的泥浆侵入信息，反演结果更加接近地层的真实状况。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本发明实施例的一种识别地层产状的方法的流程示意图；
图2为本发明实施例的一种地层参数的反演方法的流程示意图；
图3为基于电成像测井数据，获取的地层界面和地层倾角示意图；
图4为根据本发明的反演方法，获取的侵入带与原状地层电阻率、泥浆侵入深度结果示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下文中将结合附图对本发明的实施例进行详细说明。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
参照图1，根据本发明实施例的一种识别地层产状的方法，包括如下步骤：
步骤101：获取电成像仪器每个极板上的合成测井曲线；
需要说明的是，由于电成像仪器有多个极板，每个极板上有多个纽扣电极，因此，需要对每个极板上多个钮扣电极的数据进行合成。
进一步地，，所述步骤101具体包括：
步骤1011：对每个极板内若干个纽扣电极的测井曲线进行深度对齐。
具体地，由于每个极板上的纽扣电极分上下两排排列，不同排的纽扣电极在极板上的垂直位置不同，使电极响应存在深度差，数据处理时应进行深度对齐。在本发明一实施例中，对上排纽扣电极的测井曲线进行深度平移，使其深度与下排纽扣电极的测井曲线保持一致。
步骤1012：通过求取每个极板内每个采样点的多个钮扣电极的测井数据的平均值，得到所述每个极板上的合成测井曲线。
进一步地，在获取电成像仪器每块极板上的合成测井曲线之前，所述方法还包括：对每个采样点上所有钮扣电极的测井数据进行升序或降序排序，按照预设的第一剔除比例值剔除较大测井数据和/或第二剔除比例值剔除较小测井数据。
具体地，在电成像测井过程中，某一个或某几个电极可能临时性工作不正常，或者当出现裂缝、孔洞、砾岩等地质现象时，电极的测量数据不反映地层电导率的变化，因此，通过预先设置第一剔除比例值和/或第二剔除比例值剔除异常测量数据。例如，预先设置第一剔除比例值为10％，第二剔除比例值为10％，将按大小排序后的测井数据中的较大值的10％和较小值的10％全部剔除，保留中间主体数据。
步骤102：对每条合成测井曲线进行预分层，得到每条合成测井曲线的预分层界面位置；
在本发明一实施例中，通过活度分层法对每条合成测井曲线进行预分层。与现有的其它分层方法相比，活度分层法具有运算简单、计算量小、分层速度快等特点。
进一步地，所述步骤102包括如下步骤：
步骤1021：对每条合成测井曲线进行平滑滤波处理；
进一步地，通过钟形函数平滑滤波或均值滤波对每条合成测井曲线进行平滑滤波处理。
在本发明一实施例中，通过钟形函数平滑滤波处理时，假设相邻采样点之间满足二次函数变化关系，则可以用如下的二次函数曲线Zt对采样点进行拟合：
Zt＝a0+a1t+a2t2 (1)
采用相邻5个采样点做平滑计算，Ti+t为测井曲线第i+t个采样点值，t为[-2，2]之间的整数，残差Q的计算公式为：
当残差Q取最小值时，Ti+t与Zt的偏差最小，此时，通过对残差Q分别对参数a0，a1与a2求偏导数得到参数a0，a1与a2，即通过和三个方程解出。
通过最小二乘法导出滑动平均值的计算公式为：
步骤1022：获取每条合成测井曲线对应的活度曲线和对比度曲线；
每条合成测井曲线的活度曲线计算公式为：
式中，E(d)表示d点的活度函数值；x(t)表示合成测井曲线的测井值；表示x(t)在区间[d-Δ/2,d+Δ/2]内的平均值；Δ表示计算活度函数值时所取合成测井曲线的长度(下称窗长)。
d点对应的对比度值定义为：
Contrast_Ratio＝Max/Min (5)
式中Contrast_Ratio表示对比度值，Max与Min分别表示计算活度函数值时所取窗长内测井值的最大与最小值；
步骤1023：求取每条活度曲线的极大值点，并当极大值点在预设的窗长阈值范围内与其余采样点的对比度值大于预设的对比度阈值时，将所述极大值点的深度作为所述合成测井曲线的预分层界面位置；
需要说明的是，对于测井曲线的地层界面，在活度曲线上出现峰值，这表明在界面处测井曲线的变化最为激烈。如果能在测井曲线对应的活度曲线上找到活度峰值点，那么该点所对应的位置就为可能的地层界面。
步骤103：依次选取包含预分层界面位置的深度段作为参考段，对其它极板上的合成测井曲线进行相关对比，得到各个合成测井曲线的有效层界面；
具体地，以某个极板为例，所述步骤103包括如下步骤：
步骤1031：选取所述极板上的合成测井曲线作为参考曲线，将其它极板上的合成测井曲线作为对比曲线；
步骤1032：在参考曲线上，按照深度顺序依次选取三个预分层界面的深度段数据作为参考段，第一个深度点与第三个深度点之间的深度差为本次对比窗长；
步骤1033：在各对比曲线上，使用相关对比法在预设的探索长度内依次计算出最大相关系数时的对比段中点深度；
具体地，所述步骤1033具体包括：在各对比曲线上，在预设的探索长度范围内依次选取与本次对比窗长长度相等的对比段，计算各对比段与参考段的合成测井曲线的相关系数，记录每次探索长度范围内的最大相关系数和获得最大相关系数时的对比段中点深度。
步骤1034：计算满足最大相关系数大于预设的相关系数阈值的对比曲线条数，检测所述对比曲线条数是否大于或等于预设的对比曲线条数阈值，如果大于或等于预设的对比曲线条数阈值，将参考曲线中的所述预分层界面作为有效层界面，将所述获得最大相关系数时的对比段中点深度，作为与参考曲线中目标层界面相同层的层界面深度；
步骤1035：对参考曲线上的所有预分层界面循环执行步骤1032-1034，得到该极板上的合成测井曲线的所有有效层界面。
对其它极板上的合成测井曲线，依次执行步骤1031-1035，得到所有极板上的合成测井曲线的所有有效层界面。
步骤104：对所有极板上的合成测井曲线的有效层界面按深度进行排序，获得地层界面并计算地层产状信息。
在本发明一实施例中，所述步骤104还包括：在获得的地层界面中，去除不合理的有效层界面。
进一步地，所述去除不合理的层界面，具体为：当相邻两个有效层界面的深度差小于预设的深度阈值时，在获得的地层界面中，去除所述相邻两个有效层界面中相关系数小的有效层界面。
进一步地，所述去除不合理的层界面，具体为：当所述相邻两个有效层界面存在交叉时，在获得的地层界面中，去除所述相邻两个有效层界面中相关系数小的有效层界面。
进一步地，所述识别地层产状的方法还包含：根据双侧向测井数据提取地层界面，对双侧向测井数据提取的地层界面，和基于电成像测井数据获得的地层界面进行合并。
由于电成像图像上的层界面本质是不同岩性的地层电阻率不一致引起的图像颜色跳变，在电扣曲线上表现为极值点或者突变值点，本发明提供的识别地层产状的方法是在检测这些极值点或突变点的基础上，识别出地层界面，然后计算出层界面产状信息，所以提取的层界面与产状信息与图像特征吻合较好，与地层的真实状态很接近，且通过本发明计算地层产状快速高效，准确度高，基本可以取代目前的人工拾取方法。
参照图2，根据本发明实施例的一种地层参数的反演方法，包括如下步骤：
步骤201：获取电成像仪器每个极板上的合成测井曲线；
具体地，步骤201中的获取电成像仪器每个极板上的合成测井曲线的方法同前述步骤101，此处不再赘述。
步骤202：对每条合成测井曲线进行预分层，得到每条合成测井曲线的预分层界面位置；
具体地，步骤202中的对每条合成测井曲线进行预分层，得到每条合成测井曲线的预分层界面位置的方法同前述步骤102，此处不再赘述。
步骤203：依次选取包含预分层界面位置的深度段作为参考段，对其它极板上的合成测井曲线进行相关对比，得到各个合成测井曲线的有效层界面；
具体地，步骤203中的对其它极板上的合成测井曲线进行相关对比，得到各个合成测井曲线的有效层界面的方法同前述步骤103，此处不再赘述。
步骤204：对所有极板上的合成测井曲线的有效层界面按深度进行排序，获得地层界面；
具体地，步骤204中的对所有极板上的合成测井曲线的有效层界面按深度进行排序，获得地层界面的方法同前述步骤104，此处不再赘述。
步骤205：基于电成像测井数据，提取各地层界面对应的地层倾角、地层倾向；
进一步地，步骤205包含以下步骤：
步骤2051：求取地层视倾角和视倾向；
进一步地，步骤2051包括：通过最小二乘算法求取各个地层界面的幅度A、初始相位以及基线深度y0；根据地层界面的幅度A、初始相位计算地层视倾角和视倾向。
具体地，同一地层界面在不同极板合成曲线上的深度满足以下正弦关系：
式中，i表示极板序号，yi是第i个极板对应的深度，xi是第i个极板所在方位，A为幅度值也即高程差值的一半，为初始相位值，与地层视倾向有关，y0为地层界面的基线深度值。
地层视倾角和视倾向计算公式如下：
Dip＝Tan-1(2*A/Del) (7)
式中Dip为地层视倾角，Azi为地层视倾向，Del为预先设置的电直径值。
步骤2052：将地层视倾角和视倾向转化为地层真倾角和真倾向；
具体地，地层真倾角和真倾向可以由井眼倾角、井斜方位角和视倾角、视倾向经过坐标旋转变换得到。图3为基于电成像测井数据，获取的地层界面和地层倾角示意图，参照图3，从左起第一道为电成像原始测量信号成像和提取到的层界面指示曲线道，其中“N-E-S-W-N”指“北-东-南-西-北”，第二道为深度道，第三道为地层倾角蝌蚪图。
步骤206：基于电成像测井数据，提取冲洗带视电阻率信息；
进一步地，步骤206包含以下步骤：
步骤2061：求取冲洗带视电阻率值；
具体地，将多条合成测井曲线上对应采样深度的测井值加权平均后作为对应深度点的冲洗带电导率值，使用仪器发射电压除以该冲洗带电导率并乘以预先设置的仪器系数，得到冲洗带视电阻率值；
步骤2062：对冲洗带视电阻率进行重采样，获得与双侧向测井采样深度相匹配的冲洗带视电阻率。
步骤207：根据电成像、双侧向测井数据进行联合反演，得到泥浆侵入深度、侵入带地层电阻率与原状地层电阻率。
进一步地，步骤207包含以下步骤：
步骤2071：根据预先设置的泥浆侵入深度初值、原状地层电阻率初值，以及所述提取的地层界面、地层倾角、冲洗带视电阻率等地层信息创建初始正演模型；
步骤2072：进行地层仿真，得到深浅两种探测模式下的双侧向测井响应值；
进一步地，步骤2072中，当地层倾角小于预设的角度阈值时，通过二维模型进行仿真；当地层倾角大于预设的角度阈值时，通过三维模型进行仿真。
步骤2073：判断是否满足终止迭代条件，当满足终止迭代条件时，终止迭代并提取模型中的泥浆侵入深度、侵入带电阻率与原状地层电阻率信息；当不满足终止迭代条件时，修正泥浆侵入深度与原状地层电阻率参数，重复步骤2071至2073，直至满足终止迭代条件为止。
进一步地，步骤2073中，所述终止迭代条件为以下任意一个：将所述双侧向测井响应值与双侧向实测值进行比较，响应值与实测值之间的差值小于预先设置的第一差值阈值；迭代步数大于预先设置的最大迭代步数；将相邻两次仿真得到的双侧向测井响应值与双侧向实测值的差值进行比较，所述差值的变化小于预先设置的第二差值阈值。
图4为根据本发明的反演方法，获取的侵入带与原状地层电阻率、泥浆侵入深度结果示意图，参照图4，从左起第一道为电成像图像道，第二道为双侧向电阻率曲线道，其中“N-E-S-W-N”指“北-东-南-西-北”，第三道为深度道，第四道为反演出的侵入带与原状地层电阻率道，第五道为反演出的泥浆侵入深度道。从图4可以看出，本发明提供的地层参数的反演方法，所提取的地层界面及产状信息与测井图像吻合较好，解决了复杂背景干扰下地层产状难以提取的问题；本发明分辨率高，充分利用地层产状信息，反演结果更加接近地层的真实状况，计算结果更为准确。
本领域普通技术人员可以理解上述方法中的全部或部分步骤可通过程序来指令相关硬件完成，所述程序可以存储于计算机可读存储介质中，如只读存储器、磁盘或光盘等。可选地，上述实施例的全部或部分步骤也可以使用一个或多个集成电路来实现，相应地，上述实施例中的各模块/单元可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。本发明不限制于任何特定形式的硬件和软件的结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种识别地层产状的方法，其特征在于，包括：
获取电成像仪器每个极板上的合成测井曲线；
对每条合成测井曲线进行预分层，得到每条合成测井曲线的预分层界面位置；
依次选取包含预分层界面位置的深度段作为参考段，对其它极板上的合成测井曲线进行相关对比，得到各个合成测井曲线的有效层界面；
对所有极板上的合成测井曲线的有效层界面按深度进行排序，获得地层界面并计算地层产状信息。
2.根据权利要求1所述的识别地层产状的方法，其特征在于，所述获取电成像仪器每块极板上的合成测井曲线，具体包括：
对所述每个极板内若干个纽扣电极的测井曲线进行深度对齐；
对所述每个极板内每个采样点上若干个钮扣电极的测井数据求取平均值，将所述平均值作为所述采样点的测量值，从而得到所述每个极板上的合成测井曲线。
3.根据权利要求1所述的识别地层产状的方法，其特征在于，所述获取电成像仪器每块极板上的合成测井曲线之前，所述方法还包括：
对所述每个极板内每个采样点上所有钮扣电极的测井数据进行升序或降序排序，按照预设的第一剔除比例值剔除较大测井数据，以及第二剔除比例值剔除较小测井数据。
4.根据权利要求1所述的识别地层产状的方法，其特征在于，所述对每条合成测井曲线进行预分层，得到每条合成测井曲线的预分层界面位置的步骤，具体包括：
对每条所述合成测井曲线进行平滑滤波处理；
获取每条所述合成测井曲线对应的活度曲线和对比度曲线；
求取每条活度曲线的极大值点，并当极大值点在预设的窗长阈值范围内与其余采样点的对比度值大于预设的对比度阈值时，将所述极大值点的深度作为所述合成测井曲线的预分层界面位置。
5.根据权利要求1所述的识别地层产状的方法，其特征在于，所述依次选取包含预分层界面位置的深度段作为参考段，对其它的合成测井曲线进行相关对比，得到各个合成测井曲线的有效层界面的步骤，具体包括：
对所述每个极板上的合成测井曲线执行如下步骤：
选取每个极板上的所述合成测井曲线作为参考曲线，将其它极板上的所述合成测井曲线作为对比曲线；
在参考曲线上，按照深度顺序依次选取三个所述预分层界面的深度段数据作为参考段，第一个深度点与第三个深度点之间的深度差为本次对比窗长；
在各对比曲线上，使用相关对比法在预设的探索长度内依次计算出最大相关系数时的对比段中点深度；
计算满足最大相关系数大于预设的相关系数阈值的对比曲线条数，检测所述对比曲线条数是否大于或等于预设的对比曲线条数阈值，如果大于或等于预设的对比曲线条数阈值，将参考曲线中的所述预分层界面作为有效层界面，将所述获得最大相关系数时的对比段中点深度，作为与参考曲线中目标层界面相同层的层界面深度；
对参考曲线上的所有预分层界面进行相关对比，得到该极板上的合成测井曲线的所有有效层界面。
6.根据权利要求5所述的识别地层产状的方法，其特征在于，所述在各对比曲线上，使用相关对比法在预设的探索长度内依次计算出最大相关系数时的对比段中点深度的步骤，具体包括：
在各对比曲线上，在预设的探索长度范围内依次选取与本次对比窗长长度相等的对比段，计算各对比段与参考段的合成测井曲线的相关系数，记录每次探索长度范围内的最大相关系数和获得最大相关系数时的对比段中点深度。
7.根据权利要求1所述的识别地层产状的方法，其特征在于：所述方法还包括：
当相邻两个所述有效层界面的深度差小于预设的深度阈值时，或者当所述相邻两个所述有效层界面存在交叉时，在所述获得的地层界面中，去除所述相邻两个有效层界面中相关系数小的有效层界面。
8.一种地层参数的反演方法，其特征在于，包括：
获取电成像仪器每个极板上的合成测井曲线；
对每条合成测井曲线进行预分层，得到每条合成测井曲线的预分层界面位置；
依次选取包含预分层界面位置的深度段作为参考段，对其它极板上的合成测井曲线进行相关对比，得到各个合成测井曲线的有效层界面；
对所有极板上的合成测井曲线的有效层界面按深度进行排序，获得地层界面；
基于电成像测井数据，提取各地层界面对应的地层倾角、地层倾向，提取冲洗带视电阻率信息；
根据获得的地层界面、提取的地层倾角与冲洗带视电阻率信息，创建初始正演模型，进行地层仿真，将仿真获取的测井响应值与实际测井数据进行对比，判断是否满足终止迭代条件，如果不满足终止迭代条件，修正地层参数，重新进行仿真；如果满足终止迭代条件，提取全井段地层的泥浆侵入信息。
9.根据权利要求8所述的地层参数的反演方法，其特征在于，在所述根据获得的地层界面、提取的地层倾角与冲洗带视电阻率信息，创建初始正演模型，进行地层仿真的步骤中，当所述地层倾角小于预设的角度阈值时，通过二维模型进行仿真；当所述地层倾角大于预设的角度阈值时，通过三维模型进行仿真。
10.根据权利要求8所述的地层参数的反演方法，其特征在于，所述终止迭代条件包括以下任意一个：
将所述双侧向测井响应值与双侧向实测值进行比较，响应值与实测值之间的差值小于预先设置的第一差值阈值；
迭代步数大于预先设置的最大迭代步数；
将相邻两次仿真得到的双侧向测井响应值与双侧向实测值的差值进行比较，所述差值的变化小于预先设置的第二差值阈值充电器测量自身的振动数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710337945.6A CN107218033B (zh) | 2017-05-15 | 2017-05-15 | 一种识别地层产状的方法和地层参数的反演方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710337945.6A CN107218033B (zh) | 2017-05-15 | 2017-05-15 | 一种识别地层产状的方法和地层参数的反演方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107218033A true CN107218033A (zh) | 2017-09-29 |
| CN107218033B CN107218033B (zh) | 2020-10-30 |
Family
ID=59944410
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710337945.6A Active CN107218033B (zh) | 2017-05-15 | 2017-05-15 | 一种识别地层产状的方法和地层参数的反演方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107218033B (zh) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108915676A (zh) * | 2018-07-20 | 2018-11-30 | 陕西延长石油(集团)有限责任公司研究院 | 一种致密储层孔隙可动流体侵入剖面成像方法 |
| CN108952690A (zh) * | 2018-08-01 | 2018-12-07 | 中国石油大学(华东) | 基于随钻方位电磁波测井资料的地层界面实时提取方法 |
| CN110852144A (zh) * | 2019-08-23 | 2020-02-28 | 中国石油大学(北京) | 一种基于dtw的智能地层对比方法和系统 |
| CN111472765A (zh) * | 2020-04-24 | 2020-07-31 | 北京月新时代科技股份有限公司 | 目标井的地层划分方法和装置 |
| CN112253090A (zh) * | 2020-10-14 | 2021-01-22 | 中海油田服务股份有限公司 | 一种多频电成像的数据参数反演方法和装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5870690A (en) * | 1997-02-05 | 1999-02-09 | Western Atlas International, Inc. | Joint inversion processing method for resistivity and acoustic well log data |
| CN1696742A (zh) * | 2005-06-17 | 2005-11-16 | 中油测井技术服务有限责任公司 | 一种计算井筒中地层产状的测井数据处理方法 |
| CN104808248A (zh) * | 2015-04-22 | 2015-07-29 | 长江大学 | 基于电成像测井全井眼图像的地层产状检测方法 |
| CN104849766A (zh) * | 2015-05-19 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司 | 电成像测井图像的地层产状识别方法 |
| CN105114063A (zh) * | 2015-08-04 | 2015-12-02 | 中国石油大学(华东) | 一种斜井双侧向测井曲线校正方法 |
| CN106324689A (zh) * | 2016-06-24 | 2017-01-11 | 杭州迅美科技有限公司 | 一种水平井地层环境下电阻率各向异性识别方法 |
-
2017
- 2017-05-15 CN CN201710337945.6A patent/CN107218033B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5870690A (en) * | 1997-02-05 | 1999-02-09 | Western Atlas International, Inc. | Joint inversion processing method for resistivity and acoustic well log data |
| CN1696742A (zh) * | 2005-06-17 | 2005-11-16 | 中油测井技术服务有限责任公司 | 一种计算井筒中地层产状的测井数据处理方法 |
| CN104808248A (zh) * | 2015-04-22 | 2015-07-29 | 长江大学 | 基于电成像测井全井眼图像的地层产状检测方法 |
| CN104849766A (zh) * | 2015-05-19 | 2015-08-19 | 中国石油集团川庆钻探工程有限公司 | 电成像测井图像的地层产状识别方法 |
| CN105114063A (zh) * | 2015-08-04 | 2015-12-02 | 中国石油大学(华东) | 一种斜井双侧向测井曲线校正方法 |
| CN106324689A (zh) * | 2016-06-24 | 2017-01-11 | 杭州迅美科技有限公司 | 一种水平井地层环境下电阻率各向异性识别方法 |
Non-Patent Citations (1)
| Title |
|---|
| 陈岩: "新疆准噶尔盆地乌尔禾油田低渗难采油藏精细描述方法研究", 《中国优秀博硕士学位论文全文数据库（博士） 工程科技Ⅰ辑》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108915676A (zh) * | 2018-07-20 | 2018-11-30 | 陕西延长石油(集团)有限责任公司研究院 | 一种致密储层孔隙可动流体侵入剖面成像方法 |
| CN108915676B (zh) * | 2018-07-20 | 2021-09-14 | 陕西延长石油(集团)有限责任公司研究院 | 一种致密储层孔隙可动流体侵入剖面成像方法 |
| CN108952690A (zh) * | 2018-08-01 | 2018-12-07 | 中国石油大学(华东) | 基于随钻方位电磁波测井资料的地层界面实时提取方法 |
| CN108952690B (zh) * | 2018-08-01 | 2022-01-25 | 中国石油大学(华东) | 基于随钻方位电磁波测井资料的地层界面实时提取方法 |
| CN110852144A (zh) * | 2019-08-23 | 2020-02-28 | 中国石油大学(北京) | 一种基于dtw的智能地层对比方法和系统 |
| CN111472765A (zh) * | 2020-04-24 | 2020-07-31 | 北京月新时代科技股份有限公司 | 目标井的地层划分方法和装置 |
| CN111472765B (zh) * | 2020-04-24 | 2023-07-14 | 北京月新时代科技股份有限公司 | 目标井的地层划分方法和装置 |
| CN112253090A (zh) * | 2020-10-14 | 2021-01-22 | 中海油田服务股份有限公司 | 一种多频电成像的数据参数反演方法和装置 |
| CN112253090B (zh) * | 2020-10-14 | 2022-10-18 | 中海油田服务股份有限公司 | 一种多频电成像的数据参数反演方法和装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107218033B (zh) | 2020-10-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113759425B (zh) | 井震联合评价深层古岩溶储层充填特征的方法与系统 | |
| CN107218033A (zh) | 一种识别地层产状的方法和地层参数的反演方法 | |
| CN106526693B (zh) | 裂缝识别方法和装置 | |
| CN105221133B (zh) | 一种基于测井多参数确定烃源岩有机碳含量的方法和装置 | |
| US5862100A (en) | Method and system for detecting hydrocarbon reservoirs using statistical normalization of amplitude-versus-offset indicators based upon seismic signals | |
| US7062072B2 (en) | Methods of producing images of underground formations surrounding a borehole | |
| US6574566B2 (en) | Automated feature identification in data displays | |
| CN111580181B (zh) | 一种基于多场多特征信息融合的导水陷落柱识别方法 | |
| CN108802812A (zh) | 一种井震融合的地层岩性反演方法 | |
| EA019763B1 (ru) | Способ и компьютерный носитель информации для анализа характеристик подземной формации | |
| CN105938503B (zh) | 一种方向信号多层界面识别方法 | |
| CN109763814B (zh) | 基于多维测井数据的地层匹配可视分析方法 | |
| CN111596364A (zh) | 基于高精度层序地层格架的地震沉积微相组合分析方法 | |
| CN105005080A (zh) | 一种利用振幅比率属性识别地层圈闭尖灭线的方法 | |
| CN104011566A (zh) | 用于分析地质构造的特性的基于小波变换的系统和方法 | |
| CN109425900A (zh) | 一种地震储层预测方法 | |
| CN109522578A (zh) | 基于贝叶斯模糊判别的复杂岩相预测方法及系统 | |
| CN117274520A (zh) | 一种地质勘探的测绘方法 | |
| CN111025394A (zh) | 一种基于深度域地震资料精细小断层探测方法及装置 | |
| Panara et al. | Fracture intensity and associated variability: A new methodology for 3D digital outcrop model analysis of carbonate reservoirs | |
| CN112835124A (zh) | 基于成像测井和阵列声波测井资料的裂缝有效性评价方法 | |
| CN108957554A (zh) | 一种地球物理勘探中地震反演方法 | |
| CN108035710A (zh) | 基于数据挖掘划分深层岩石地质相的方法 | |
| CN115857047B (zh) | 一种地震储层综合预测方法 | |
| CN113075748B (zh) | 基于成像测井和声波远探测测井资料的裂缝有效性评价方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB02 | Change of applicant information | ||
| CB02 | Change of applicant information |
Address after: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing Applicant after: China Offshore Oil Group Co., Ltd. Applicant after: China Oilfield Services Limited Address before: 100010 Chaoyangmen North Street, Dongcheng District, Dongcheng District, Beijing Applicant before: China National Offshore Oil Corporation Applicant before: China Oilfield Services Limited |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |