CN111123378B - 确定划分岩性类型的伽马射线强度临界值的方法及装置 - Google Patents

确定划分岩性类型的伽马射线强度临界值的方法及装置 Download PDF

Info

Publication number
CN111123378B
CN111123378B CN201911354826.7A CN201911354826A CN111123378B CN 111123378 B CN111123378 B CN 111123378B CN 201911354826 A CN201911354826 A CN 201911354826A CN 111123378 B CN111123378 B CN 111123378B
Authority
CN
China
Prior art keywords
value
gamma ray
ray intensity
target
oil
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
Application number
CN201911354826.7A
Other languages
English (en)
Other versions
CN111123378A (zh
Inventor
黄平
黄科
安辉
邹定永
张福宏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201911354826.7A priority Critical patent/CN111123378B/zh
Publication of CN111123378A publication Critical patent/CN111123378A/zh
Application granted granted Critical
Publication of CN111123378B publication Critical patent/CN111123378B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
    • G01V5/04Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种确定划分岩性类型的伽马射线强度临界值的方法及装置,属于石油地质勘探技术领域。方法包括:获取目标区域中每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线,确定每个油气井的伽马射线强度值与泥质含量值的关系,获取目标区域的检测泥质含量值,确定目标区域中的目标油气井在检测泥质含量值下的目标伽马射线强度值,将目标伽马射线强度值作为目标区域的划分岩性类型的伽马射线强度临界值。通过上述方法确定划分岩性类型的伽马射线强度临界值,考虑到了目标区域的检测泥质含量值因素,能够更好的反映目标区域的岩性类型分布情况,所以基于该伽马射线强度临界值划分岩性类型结果更加准确。

Description

确定划分岩性类型的伽马射线强度临界值的方法及装置
技术领域
本发明涉及石油地质勘探技术领域,特别涉及一种确定划分岩性类型的伽马射线强度临界值的方法及装置。
背景技术
为了对油气藏进行高效钻探,作业人员需要对油气藏所在储层的岩性类型进行划分。储层岩石中所含放射性元素衰变时会放出伽马射线,不同的储层岩石所产生的伽马射线强度值不同,因此,根据伽马射线强度值可以划分岩性类型。
相关技术中,统计每一个油气井不同深度的储层岩石所对应的伽马射线强度值,利用离散型频率直方图法绘制每一个油气井的伽马射线强度值概率分布直方图。根据该图绘制出关于伽马射线强度值概率分布的曲线,确定伽马射线强度值大概率分布的中心位置,该中心位置处所对应的伽马射线强度值作为该油气井的伽马射线强度值。然后在这些伽马射线强度值中,选取一个作为确定划分岩性类型的伽马射线强度临界值。
通过上述方法获得的伽马射线强度临界值,可能使各油气井的伽马射线强度临界值所对应的测井解释岩性所采用的泥质含量值不一致,导致在目标区域内油气井与油气井之间岩性解释标准与测井解释不一致,岩性类型划分结果不够准确。
发明内容
本发明实施例提供了一种确定划分岩性类型的伽马射线强度临界值的方法及装置,可以解决相关技术中岩性类型划分结果不够准确的问题。所述技术方案如下:
一方面,提供了一种确定划分岩性类型的伽马射线强度临界值的方法,所述方法包括:
获取目标区域中每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线;
根据所述每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线,确定每个油气井的伽马射线强度值与泥质含量值的关系;
获取所述目标区域的检测泥质含量值;
基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,将所述目标伽马射线强度值作为所述目标区域的划分岩性类型的伽马射线强度临界值。
在一种可能的实现方式中,所述基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,包括:
基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量的关系,确定每个油气井检测泥质含量值对应的伽马射线强度值;
基于所述每个油气井检测泥质含量值对应的伽马射线强度值,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,所述基于所述每个油气井检测泥质含量值对应的伽马射线强度值,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,包括:
如果确定出的伽马射线强度值的个数为奇数个,将奇数个伽马射线强度值按数值大小顺序排列,选取排在中间一个位置的伽马射线强度值,作为所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,所述基于所述每个油气井检测泥质含量值对应的伽马射线强度值,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,还包括:
如果确定出的伽马射线强度值的个数为偶数个,将偶数个伽马射线强度值按数值大小顺序排列,在排在中间两个位置的伽马射线强度值中,随机选取一个伽马射线强度值,作为所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,所述基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,包括:
将所述目标区域中的评价井确定为目标油气井;
基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,所述基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,将所述目标伽马射线强度值作为所述目标区域的划分岩性类型的伽马射线强度临界值,包括:
获取所述目标区域中每个油气井在所述检测泥质含量值下的伽马射线强度值;
根据所述目标油气井的目标伽马射线强度值分别与每个油气井的伽马射线强度值的差值,平移每个油气井的伽马自然测井曲线;
将所述目标油气井的伽马自然测井曲线和平移得到的每个油气井的伽马自然测井曲线,输入到地震反演模型中,得到用于体现所述目标区域不同位置的岩性类型的地质模型。
另一方面,提供了确定划分岩性类型的伽马射线强度临界值的装置,所述装置包括:
第一获取模块,用于获取目标区域中每个油气井的自然伽马测井曲线、以及泥质含量与深度的关系曲线;
第一确定模块,用于根据所述每个油气井的自然伽马测井曲线、以及泥质含量与深度的关系曲线,确定每个油气井的伽马射线强度值与泥质含量的关系;
第二获取模块,用于获取所述目标区域的检测泥质含量值;
第二确定模块,用于基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,将所述目标伽马射线强度值作为所述目标区域的划分岩性类型的伽马射线强度临界值。
在一种可能的实现方式中,所述第二确定模块,包括:
第一确定子单元,用于基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量的关系,确定每个油气井检测泥质含量值对应的伽马射线强度值;
第二确定子单元,用于基于所述每个油气井检测泥质含量值对应的伽马射线强度值,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,所述第二确定子单元用于:
如果确定出的伽马射线强度值的个数为奇数个,将奇数个伽马射线强度值按数值大小顺序排列,选取排在中间一个位置的伽马射线强度值,作为所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,所述第二确定子单元还用于:
如果确定出的伽马射线强度值的个数为偶数个,将偶数个伽马射线强度值按数值大小顺序排列,在排在中间两个位置的伽马射线强度值中,随机选取一个伽马射线强度值,作为所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,所述第二确定模块用于:
将所述目标区域中的评价井确定为目标油气井;
基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,所述基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,将所述目标伽马射线强度值作为所述目标区域的划分岩性类型的伽马射线强度临界值,包括:
获取所述目标区域中每个油气井在所述检测泥质含量值下的伽马射线强度值;
根据所述目标油气井的目标伽马射线强度值分别与每个油气井的伽马射线强度值的差值,平移每个油气井的伽马自然测井曲线;
将所述目标油气井的伽马自然测井曲线和平移得到的每个油气井的伽马自然测井曲线,输入到地震反演模型中,得到用于体现所述目标区域不同位置的岩性类型的地质模型。
本发明的实施例提供的技术方案至少带来以下有益效果:
本发明实施例通过获取目标区域中每个油气井的自然伽马测井曲线、以及泥质含量与深度的关系曲线,根据每个油气井的自然伽马测井曲线、以及泥质含量与深度的关系曲线,确定每个油气井的伽马射线强度值与泥质含量的关系,获取目标区域的检测泥质含量值,基于检测泥质含量值、以及每个油气井的伽马射线强度值与泥质含量的关系,确定目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值,将该目标伽马射线强度值作为目标区域的划分岩性类型的伽马射线强度临界值。通过上述方法确定划分岩性类型的伽马射线强度临界值,考虑到了目标区域的检测泥质含量值因素,能够更好的反映目标区域的岩性类型分布情况,所以基于该伽马射线强度临界值划分岩性类型结果更加准确。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种确定划分岩性类型的伽马射线强度临界值的方法流程图;
图2是本发明实施例提供的一种确定划分岩性类型的伽马射线强度临界值的装置结构示意图;
图3是本发明实施例提供的一种电子设备的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
本发明实施例提供的方法,可以应用在石油地质勘探技术领域。具体的用于划分岩性类型。当地震技术人员想要预测目标区域内的砂岩和泥岩分布范围,即需要划分岩性类型时,可以预先收集目标区域中每个油气井的泥质含量与深度关系曲线及自然伽马测井曲线,通过分析得到每个油气井的伽马射线强度值与泥质含量的关系。根据目标区域的地质情况,确定该目标区域的检测泥质含量值,进而确定每个油气井在该检测泥质含量值下的伽马射线强度值。在目标区域中,选取一个目标油气井,例如,可以直接将资料比较齐全、井壁没有垮塌的评价井,作为目标油气井。也可以通过计算伽马射线强度值的概率分布范围来确定目标油气井。最后,根据检测泥质含量值、以及每个油气井的伽马射线强度值与泥质含量的关系,确定目标区域中的目标油气井在检测泥质含量值下的目标伽马射线强度值,将该目标伽马射线强度值作为该目标区域的划分岩性类型的伽马射线强度临界值。
图1是本发明实施例提供的一种确定划分岩性类型的伽马射线强度临界值的方法流程图,参见图1,该实施例包括:
步骤101、获取目标区域中每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线。
其中,泥质是指岩石中颗粒含量小于0.01mm(单位:毫米)的碎屑物质,而泥质含量值指的是泥质的体积含量占岩石总体积含量的比。
自然伽马测井曲线是通过将井下仪器放入油气井中,测量并记录每一深度点上接收的自然伽马射线总强度,从而获得的锯齿状不规则的曲线。在石油地质勘探和开发中,由于储层岩石中所含放射性元素衰变时会放出伽马射线,不同岩性的岩层含有的放射性元素的种类和数量不同,所产生的伽马射线强度也不同,利用自然伽马测井曲线可得到伽马射线强度值与油气井深度的关系。
泥质含量值与深度的关系曲线由测井工程师根据钻井时提取的岩样作为约束条件,用计算自然伽马射线强度值的方式获得。具体的是,可以通过测量和记录井剖面上岩层的自然伽马性质,获得油气井的每一深度点上的自然伽马值,用钻井时提取的岩石样品对自然伽马值进行标定,计算泥质在整个岩石地层中含量的百分比,该百分比即为泥质含量值,可以反映上、下岩石岩性的差异性。其中,自然伽马性质包括:岩石所含自然伽马射线强度值等。由于岩石中的矿物成分即岩性不同,造成不同的岩石处的自然伽马存在差异。利用自然伽马测井可以获得自然伽马射线强度值,通过伽马射线强度值可以获得大致的岩性类型,例如,在砂泥岩剖面中,砂岩伽马值较小,泥岩伽马值较大。
在实施中,自然伽马测井曲线的获取过程为:主要使用地面仪器和井下仪器相结合的方式去探测自然伽马射线强度。其中,地面仪器主要包括电源、记录仪以及将来自井下的一连串电脉冲转换成连续电流的一整套电路等。井下仪器主要有伽马射线探测器、供给该伽马射线探测器所需的高压电源以及将输出的电脉冲进行放大的放大器等。当井下仪器在井内由下而上提升时,来自岩层的自然伽马射线穿过井内泥浆和仪器外壳进入伽马射线探测器,该伽马射线探测器接收到的一连串伽马射线转换成一个个的电脉冲,然后经过井下放大器加以放大,由电缆线送至地面仪器,地面仪器把每分钟接收到的电脉冲数转变成与其成比例的点位差由记录仪进行记录。由此,井下仪器在井内自下而上移动测量,就连续记录出井剖面的自然伽马测井曲线了。
在实施中,利用自然伽马测井获得泥质含量值与深度的关系曲线的过程为:由测井工程师根据钻井时提取的岩样作为约束条件,用计算自然伽马射线强度值的方式获得。具体的是,可以通过测量和记录井剖面上岩层的自然伽马射线强度性质,获得油气井的每一深度点上的自然伽马射线强度值,用钻井时提取的岩石样品对自然伽马射线强度值进行标定,计算井剖面每一深度点泥质含量的百分比,可得到井剖面每一深度点的泥质含量。根据该泥质含量值可以划分每口油气井的井剖面上的岩性类型。
步骤102、根据该每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线,确定每个油气井的伽马射线强度值与泥质含量值的关系。
在实施中,根据自然伽马测井曲线可以获得油气井的每一深度点上的伽马射线强度值,,根据泥质含量值与深度的关系曲线,可以获得油气井的每一深度点上的泥质含量值。根据获得的自然伽马测井曲线、以及泥质含量值与深度的关系曲线,可以确定每个油气井的伽马射线强度值与泥质含量值的关系。
步骤103、获取该目标区域的检测泥质含量值。
其中,检测泥质含量值是指该目标区域泥质含量值的临界值。
在实施中,测井工程师可以收集该目标区域内已钻井的测井数据,直接得到目标区域的检测泥质含量值。
例如,测井数据可以包括岩心分析数据、是否获得工业气流等数据。根据这些测井数据以确定岩石的性质,从而确定出该目标区域的检测泥质含量值。其中,岩心分析是指按一定要求在岩石样本的岩心上钻取直径为25mm的标准岩心柱,用来测定岩石的孔隙度、渗透率等参数的标准样品。岩心分析主要包括岩心渗透率、孔隙度及含油水饱和度等数据分析。工业气流指工业气井通过酸化、压裂等增产措施后,稳定产量达到或超过规定下限的工业气流。
步骤104、基于该检测泥质含量值、以及该每个油气井的伽马射线强度值与泥质含量值的关系,确定该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值,将该目标伽马射线强度值作为该目标区域的划分岩性类型的伽马射线强度临界值。
其中,岩性是指反映储层中岩石特征的一些属性,例如,岩石的颜色、岩石的成分、岩石的结构、胶结物或特殊矿物等。岩性类型主要包括两大类:砂岩和泥岩,砂岩主要包括:砂质砂岩、粉砂岩和细砂岩等,泥岩主要包括:粉砂质泥岩和页岩等。
为了解整个目标区域的岩性分布,需要用地震伽马反演获取目标伽马射线强度值,并用该目标伽马射线强度值对目标区域进行岩性类型划分。而地震伽马反演需要用每口油气井的自然伽马曲线放入构造模型中,建立岩石地质模型。为适应多个油气井对比,需要在该目标区域内对每个油气井的测井响应特征进行统一标定,尽可能还原储层特征,准确的划分该目标区域的岩性类型。此时,需要选取一个目标油气井作为标准井,对每个油气井进行测井响应特征的统一标定。测井响应特征可以有很多,例如,自然伽马强度值、井径、岩性密度和声波时差等测井响应特征。
在实施中,目标油气井的选取原则可以包括:目标油气井必须能在整体上体现出目标区域储层地质特征的变化趋势。自然伽马测井曲线、以及每个油气井的伽马射线强度值与泥质含量值的关系曲线的测井响应特征与储层岩性相符合,且测井响应特征明显。若存在异常测井响应特征,例如,井壁垮塌等,引起的误差应尽可能要小。该油气井所在储层在目标区域内沉积具有稳定性,以便进行统一的对比分析。
在实施中,根据检测泥质含量值、以及该每个油气井的伽马射线强度值与泥质含量值的关系,确定每个油气井检测泥质含量值对应的伽马射线强度值,基于该每个油气井检测泥质含量值对应的伽马射线强度值,确定该目标区域中的目标油气井在该检测泥质含量值值下的目标伽马射线强度值。
例如,根据获得的检测泥质含量值、以及该每个油气井的伽马射线强度值与泥质含量值的对应关系,确定每个油气井在该检测泥质含量值下对应的伽马射线强度值。
在实施中,基于该每个油气井检测泥质含量值对应的伽马射线强度值,确定该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值的过程可以为:
在一种可能的实现方式中,如果确定出的伽马射线强度值的个数为奇数个,将奇数个伽马射线强度值按数值大小顺序排列,选取排在中间一个位置的伽马射线强度值,作为该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值。
需要说明的是,伽马射线强度值按数值大小顺序排列,具体的,伽马射线强度值按数值由大到小的顺序排列,也可以按数值由小到大的顺序排列。本发明实施例对此不做限制。
例如,该目标区域内有3个油气井,这3个油气井在检测泥质含量值下的伽马射线强度值分别是11API(单位:刻度单位,自然伽马测井曲线的1/200定义为API)、14API和20API,这些伽马射线强度值的个数为3个,将这3个伽马射线强度值按数值由小到大的顺序排列:11API、14API、20API,选取排在中间一个位置的伽马射线强度值14API,则根据该伽马射线强度值与油气井的对应关系,将伽马射线强度值为14API对应的油气井,确定为目标油气井。并将14API作为该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,如果确定出的伽马射线强度值的个数为偶数个,将偶数个伽马射线强度值按数值大小顺序排列,在排在中间两个位置的伽马射线强度值中,随机选取一个伽马射线强度值,作为该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值。
例如,该目标区域内有4个油气井,这4个油气井在检测泥质含量值下的伽马射线强度值分别是10API、16API、15API和20API,这些伽马射线强度值的个数为4个,将这4个伽马射线强度值按数值由小到大的顺序排列:10API、15API、16API、20API,在排在中间两个位置的伽马射线强度值15API和16API中,随机选取一个伽马射线强度值。如果选择伽马射线强度值为15API,则根据该伽马射线强度值与油气井的对应关系,将伽马射线强度值为15API对应的油气井,确定为目标油气井。并将15API作为该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值。
如果选择伽马射线强度值为16API,则根据该伽马射线强度值与油气井的对应关系,将伽马射线强度值为16API对应的油气井,确定为目标油气井。并将16API作为该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值。
在实施中,根据该每个油气井对应的伽马射线强度值,确定该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值的过程也可以为:将该目标区域中的评价井确定为目标油气井,基于该检测泥质含量值、以及该每个油气井的伽马射线强度值与泥质含量值的关系曲线,确定该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值。
其中,评价井是指在三维地震评价的基础上,且在已获得工业油气流的圈闭上,为获得油气藏类型、油气藏构造形态、油气藏厚度及物性变化,评价油气藏区域的规模、产能及经济价值,以建立探明储量为目的而钻的探井。评价井的井壁没有垮塌且测井数据资料比较齐全。
例如,将该目标区域中的评价井确定为目标油气井,从评价井的测井数据资料中可以确定该评价井的伽马射线强度值与泥质含量值的关系,在获得的每个油气井的伽马射线强度值与泥质含量值的关系中,找到在目标油气井的伽马射线强度值与泥质含量值的关系,根据检测泥质含量值,在目标油气井的伽马射线强度值与泥质含量值的关系曲线上,找到该目标油气井在该检测泥质含量值下对应的伽马射线强度值,即确定出目标油气井对应的目标伽马射线强度值。
在实施中,该基于该检测泥质含量值、以及该每个油气井的伽马射线强度值与泥质含量值的关系曲线,确定该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值,将该目标伽马射线强度值作为该目标区域的划分岩性类型的伽马射线强度临界值之后,可以获取该区域中每个油气井在该检测泥质含量值下的伽马射线强度值,根据该目标油气井的目标伽马射线强度值分别与每个油气井的伽马射线强度值的差值,平移每个油气井的伽马自然测井曲线,将该目标油气井的伽马自然测井曲线和平移得到的每个油气井的伽马自然测井曲线,输入到地震反演模型中,得到用于体现该目标区域不同位置的岩性类型的地质模型。
其中,地震反演模型是指利用地表观测地震资料,以已知地质规律和钻井、测井数据资料为约束,对地下岩层空间结构和物理性质进行成像的过程。为了使地震资料能与测井数据资料直接对比,即把界面型的反射剖面转换成岩层型的测井剖面。
由于测井仪器的差异及测井工作环境的不同,有可能使相同的储层体现出不同的测井响应特征,由此产生系统误差。由于井径、井壁粗糙度、泥浆密度与矿化度或测井仪器外径等非地层因素下的测井环境,会引起自然伽马测井曲线存在一些假象,例如,在油气井的井眼及泥浆质量不好的情况下,自然伽马测井曲线会发生严重的变形。因此,对自然伽马测井曲线标准化之前,必须进行测井环境的校正。并且,目标区域中的油气井是很少的,为了消除伽马自然测井曲线的系统误差且将这个目标区域的岩性了解清楚,需要在目标区域内利用目标油气井对每个油气井的测井响应特征进行统一标定,即对每个油气井的伽马自然测井曲线进行标准化。最后,可以将该标准化的伽马自然测井曲线作为地震反演模型的输入参量。
在实施中,技术人员可以获取到目标区域中的目标油气井在检测泥质含量值下的目标伽马射线强度值,以及中每个油气井在所述检测泥质含量值下的伽马射线强度值,然后根据该目标油气井的目标伽马射线强度值和该每个油气井的伽马射线强度值的差值,以伽马射线强度值为横坐标轴方向上平移该每个油气井的伽马自然测井曲线,即在该检测泥质含量值下,可以在目标油气井的伽马自然测井曲线和每个油气井的伽马自然测井曲线上,分别找到该泥质含量值对应的目标油气井的目标伽马射线强度值和每个油气井的伽马射线强度值,使得目标油气井的目标伽马射线强度值和每个油气井的伽马射线强度值对应的两个点能够重合。再将该目标油气井的伽马自然测井曲线和平移动得到的该每个油气井的伽马自然测井曲线,输入到地震反演模型中,进行反演计算,得到用于体现该目标区域不同位置的岩性类型的地质模型,根据该地质模型可以确定该目标区域的岩性类型分布。
例如,如果目标区域的检测泥质含量值为25%,h1在检测泥质含量值为25%下的伽马射线强度值为90API,h5油气井在检测泥质含量值为25%下的伽马射线强度值为74API。如果以h1油气井为目标油气井,即h1油气井作为基准井,将h5油气井的伽马射线强度值校正到h1油气井的伽马射线强度值上,则h5油气井的伽马自然测井曲线上的每一个伽马射线强度值均需要加上16API,即可以把h5油气井的伽马自然测井曲线平行向右移动16个API单位长度。获得h5油气井平行移动后的伽马自然测井曲线,即获得h5油气井标准化的伽马自然测井曲线。将h1油气井的伽马自然测井曲线和h5油气井标准化的伽马自然测井曲线,输入到地震反演模型中,进行反演计算,得到用于体现该目标区域不同位置的岩性类型的地质模型,根据该地质模型可以确定该目标区域的岩性类型分布。通过地震反演模型中得到的地质模型的伽马射线强度值与本发明实施例提供的方法确定出的伽马射线强度值保持一致,因此,基于该伽马射线强度临界值划分岩性类型结果更加准确。
本发明实施例提供的方法,通过获取目标区域中每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线,根据该每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线,确定每个油气井的伽马射线强度值与泥质含量值的关系,获取该目标区域的检测泥质含量值,基于该检测泥质含量值、以及该每个油气井的伽马射线强度值与泥质含量值的关系曲线,确定该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值,将该目标伽马射线强度值作为该目标区域的划分岩性类型的伽马射线强度临界值。通过上述方法确定划分岩性类型的伽马射线强度临界值,考虑到了目标区域的检测泥质含量值因素,能够更好的反映目标区域的岩性类型分布情况,所以基于该伽马射线强度临界值划分岩性类型结果更加准确。
上述所有可选技术方案,可以采用任意结合形成本发明的可选实施例,在此不再一一赘述。
图2是本发明实施例提供的一种确定划分岩性类型的伽马射线强度临界值的装置结构示意图,该方法可以应用于电子设备,参见图2,该装置包括:
第一获取模块201,用于获取目标区域中每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线;
第一确定模块202,用于根据该每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线,确定每个油气井的伽马射线强度值与泥质含量值的关系;
第二获取模块203,用于获取该目标区域的检测泥质含量值;
第二确定模块204,用于基于该检测泥质含量值、以及该每个油气井的伽马射线强度值与泥质含量值的关系曲线,确定该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值,将该目标伽马射线强度值作为该目标区域的划分岩性类型的伽马射线强度临界值。
在一种可能的实现方式中,该第二确定模块204,包括:
第一确定子单元,用于基于该检测泥质含量值、以及该每个油气井的伽马射线强度值与泥质含量值的关系曲线,确定每个油气井对应的伽马射线强度值;
第二确定子单元,用于基于该每个油气井对应的伽马射线强度值,确定该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,该第二确定子单元用于:
如果确定出的伽马射线强度值的个数为奇数个,将奇数个伽马射线强度值按数值大小顺序排列,选取排在中间一个位置的伽马射线强度值,作为该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,该第二确定子单元还用于:
如果确定出的伽马射线强度值的个数为偶数个,将偶数个伽马射线强度值按数值大小顺序排列,在排在中间两个位置的伽马射线强度值中,随机选取一个伽马射线强度值,作为该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,该第二确定模块204用于:
将该目标区域中的评价井确定为目标油气井;
基于该检测泥质含量值、以及该每个油气井的伽马射线强度值与泥质含量值的关系,确定该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值。
在一种可能的实现方式中,该基于该检测泥质含量值、以及该每个油气井的伽马射线强度值与泥质含量值的关系曲线,确定该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值,将该目标伽马射线强度值作为该目标区域的划分岩性类型的伽马射线强度临界值之后,包括:
获取该目标区域中每个油气井在该检测泥质含量值下的伽马射线强度值;
根据该目标油气井的目标伽马射线强度值分别与每个油气井的伽马射线强度值的差值,平移每个油气井的伽马自然测井曲线;
将该目标油气井的伽马自然测井曲线和平移得到的每个油气井的伽马自然测井曲线,输入到地震反演模型中,得到用于体现该目标区域不同位置的岩性类型的地质模型。
本发明实施例提供的装置,通过获取目标区域中每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线,根据该每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线,确定每个油气井的伽马射线强度值与泥质含量值的关系,获取该目标区域的检测泥质含量值,基于该检测泥质含量值、以及该每个油气井的伽马射线强度值与泥质含量值的关系,确定该目标区域中的目标油气井在该检测泥质含量值下的目标伽马射线强度值,将该目标伽马射线强度值作为该目标区域的划分岩性类型的伽马射线强度临界值。通过上述方法确定划分岩性类型的伽马射线强度临界值,考虑到了目标区域的检测泥质含量值因素,能够更好的反映目标区域的岩性类型分布情况,所以基于该伽马射线强度临界值划分岩性类型结果更加准确。
需要说明的是:上述实施例提供的确定划分岩性类型的伽马射线强度临界值的表征装置在进行确定划分岩性类型的伽马射线强度临界值时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的确定划分岩性类型的伽马射线强度临界值的装置与确定划分岩性类型的伽马射线强度临界值的方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
图3本发明实施例提供的一种电子设备的结构示意图,该电子设备300可因配置或性能不同而产生比较大的差异,可以包括一个或一个以上处理器(central processingunits,CPU)301和一个或一个以上的存储器302,其中,上述存储器302中存储有至少一条指令,上述至少一条指令由上述处理器301加载并执行以实现上述各个实施例提供的确定划分岩性类型的伽马射线强度临界值的方法。当然,该电子设备还可以具有有线或无线网络接口、键盘以及输入输出接口等部件,以便进行输入输出,该电子设备还可以包括其他用于实现设备功能的部件,在此不做赘述。
在示例性实施例中,还提供了一种计算机可读存储介质,例如包括指令的存储器,上述指令可由处理器执行以完成上述实施例中的确定划分岩性类型的伽马射线强度临界值的方法。例如,该计算机可读存储介质可以是只读存储器(Read-Only Memory,ROM)、随机存取存储器(Random Access Memory,RAM)、只读光盘(Compact Disc Read-Only Memory,CD-ROM)、磁带、软盘和光数据存储设备等。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,该程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
上述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种确定划分岩性类型的伽马射线强度临界值的方法,其特征在于,所述方法包括:
获取目标区域中每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线;
根据所述每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线,确定每个油气井的伽马射线强度值与泥质含量值的关系;
获取所述目标区域的检测泥质含量值,所述检测泥质含量值为所述目标区域的泥质含量值的临界值;
基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,将所述目标伽马射线强度值作为所述目标区域的划分岩性类型的伽马射线强度临界值。
2.根据权利要求1所述的方法,其特征在于,所述基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量值的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,包括:
基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量值的关系,确定每个油气井检测泥质含量值对应的伽马射线强度值;
基于所述每个油气井检测泥质含量值对应的伽马射线强度值,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
3.根据权利要求2所述的方法,其特征在于,所述基于所述每个油气井检测泥质含量值对应的伽马射线强度值,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,包括:
如果确定出的伽马射线强度值的个数为奇数个,将奇数个伽马射线强度值按数值大小顺序排列,选取排在中间一个位置的伽马射线强度值,作为所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
4.根据权利要求2所述的方法,其特征在于,所述基于所述每个油气井检测泥质含量值对应的伽马射线强度值,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,还包括:
如果确定出的伽马射线强度值的个数为偶数个,将偶数个伽马射线强度值按数值大小顺序排列,在排在中间两个位置的伽马射线强度值中,随机选取一个伽马射线强度值,作为所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
5.根据权利要求1所述的方法,其特征在于,所述基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量值的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,包括:
将所述目标区域中的评价井确定为目标油气井;
基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量值的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
6.根据权利要求1所述的方法,其特征在于,所述基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量值的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,将所述目标伽马射线强度值作为所述目标区域的划分岩性类型的伽马射线强度临界值之后,包括:
获取所述目标区域中每个油气井在所述检测泥质含量值下的伽马射线强度值;
根据所述目标油气井的目标伽马射线强度值分别与每个油气井的伽马射线强度值的差值,平移每个油气井的伽马自然测井曲线;
将所述目标油气井的伽马自然测井曲线和平移得到的每个油气井的伽马自然测井曲线,输入到地震反演模型中,得到用于体现所述目标区域不同位置的岩性类型的地质模型。
7.一种确定划分岩性类型的伽马射线强度临界值的装置,其特征在于,所述装置包括:
第一获取模块,用于获取目标区域中每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线;
第一确定模块,用于根据所述每个油气井的自然伽马测井曲线、以及泥质含量值与深度的关系曲线,确定每个油气井的伽马射线强度值与泥质含量值的关系;
第二获取模块,用于获取所述目标区域的检测泥质含量值,所述检测泥质含量值为所述目标区域的泥质含量值的临界值;
第二确定模块,用于基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量值的关系,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值,将所述目标伽马射线强度值作为所述目标区域的划分岩性类型的伽马射线强度临界值。
8.根据权利要求7所述的装置,其特征在于,所述第二确定模块,包括:
第一确定子单元,用于基于所述检测泥质含量值、以及所述每个油气井的伽马射线强度值与泥质含量值的关系,确定每个油气井对应的伽马射线强度值;
第二确定子单元,用于基于所述每个油气井检测泥质含量值对应的伽马射线强度值,确定所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
9.根据权利要求8所述的装置,其特征在于,所述第二确定子单元用于:
如果确定出的伽马射线强度值的个数为奇数个,将奇数个伽马射线强度值按数值大小顺序排列,选取排在中间一个位置的伽马射线强度值,作为所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
10.根据权利要求8所述的装置,其特征在于,所述第二确定子单元还用于:
如果确定出的伽马射线强度值的个数为偶数个,将偶数个伽马射线强度值按数值大小顺序排列,在排在中间两个位置的伽马射线强度值中,随机选取一个伽马射线强度值,作为所述目标区域中的目标油气井在所述检测泥质含量值下的目标伽马射线强度值。
CN201911354826.7A 2019-12-25 2019-12-25 确定划分岩性类型的伽马射线强度临界值的方法及装置 Active CN111123378B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911354826.7A CN111123378B (zh) 2019-12-25 2019-12-25 确定划分岩性类型的伽马射线强度临界值的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911354826.7A CN111123378B (zh) 2019-12-25 2019-12-25 确定划分岩性类型的伽马射线强度临界值的方法及装置

Publications (2)

Publication Number Publication Date
CN111123378A CN111123378A (zh) 2020-05-08
CN111123378B true CN111123378B (zh) 2022-06-03

Family

ID=70503260

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911354826.7A Active CN111123378B (zh) 2019-12-25 2019-12-25 确定划分岩性类型的伽马射线强度临界值的方法及装置

Country Status (1)

Country Link
CN (1) CN111123378B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103809217A (zh) * 2014-01-20 2014-05-21 西安石油大学 一种计算高、低伽马交互储层泥质含量的方法
CN104636588A (zh) * 2014-08-25 2015-05-20 中国石油天然气股份有限公司 计算烃源岩中有机碳含量的方法及装置
CN104898161A (zh) * 2014-03-05 2015-09-09 中国石油集团东方地球物理勘探有限责任公司 一种基于测井响应模拟体的有效砂岩预测方法
CN109814174A (zh) * 2019-01-24 2019-05-28 中国石油大学(华东) 一种碎屑岩不整合结构体测井综合识别方法

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA1145066A (en) * 1980-03-12 1983-04-19 Jeffrey S. Schweitzer Lithology determination from the calcium and magnesium activation lines
GB2438121B (en) * 2005-02-19 2010-11-17 Baker Hughes Inc Use of the dynamic downhole measurements as lithology indicators
CN101501531A (zh) * 2006-04-19 2009-08-05 贝克休斯公司 用于次表地层的定量岩性和矿物性评估的方法
US20080162050A1 (en) * 2006-12-28 2008-07-03 Peter Harris Method for interpreting seismic data and controlled source electromagnetic data to estimate subsurface reservoir properties
US8899348B2 (en) * 2009-10-16 2014-12-02 Weatherford/Lamb, Inc. Surface gas evaluation during controlled pressure drilling
CN102279417A (zh) * 2011-04-11 2011-12-14 中国石油化工集团公司 一种近地表岩性探测仪及探测方法
EP2749731A1 (en) * 2012-12-28 2014-07-02 Services Pétroliers Schlumberger Integrated wellbore caliper
CN105093304A (zh) * 2014-05-15 2015-11-25 中国石油化工股份有限公司 一种地球物理勘探利用测井曲线自动求取岩性曲线的方法
WO2017024530A1 (zh) * 2015-08-11 2017-02-16 深圳朝伟达科技有限公司 一种计算烃源岩中有机碳含量的方法
CN108829980B (zh) * 2018-06-20 2022-06-07 西南石油大学 利用pnn测井资料建立碳氧比和碳氢比解释模型的方法
CN109521487B (zh) * 2018-11-08 2020-03-06 中国石油大学(华东) 一种利用元素伽马能谱测井识别气层的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103809217A (zh) * 2014-01-20 2014-05-21 西安石油大学 一种计算高、低伽马交互储层泥质含量的方法
CN104898161A (zh) * 2014-03-05 2015-09-09 中国石油集团东方地球物理勘探有限责任公司 一种基于测井响应模拟体的有效砂岩预测方法
CN104636588A (zh) * 2014-08-25 2015-05-20 中国石油天然气股份有限公司 计算烃源岩中有机碳含量的方法及装置
CN109814174A (zh) * 2019-01-24 2019-05-28 中国石油大学(华东) 一种碎屑岩不整合结构体测井综合识别方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
自然伽马能谱测井在春光区块的应用研究;刘洋;《工程地球物理学报》;20180531;第15卷(第3期);第299-307页 *

Also Published As

Publication number Publication date
CN111123378A (zh) 2020-05-08

Similar Documents

Publication Publication Date Title
US9939548B2 (en) Systems, methods, and computer medium to produce efficient, consistent, and high-confidence image-based electrofacies analysis in stratigraphic interpretations across multiple wells
US6816787B2 (en) Generating and displaying a virtual core and a virtual plug associated with a selected piece of the virtual core
US8244473B2 (en) System and method for automated data analysis and parameter selection
US9097821B2 (en) Integrated workflow or method for petrophysical rock typing in carbonates
US11802985B2 (en) Method and system for analyzing filling for karst reservoir based on spectrum decomposition and machine learning
US20160349389A1 (en) Method for developing a geomechanical model based on seismic data, well logs and sem analysis of horizontal and vertical drill cuttings
US4071755A (en) Method for in situ evaluation of the source rock potential of earth formations
US9613253B2 (en) Pore-space measurement method apparatus and system
AU2011302599B2 (en) System and method for hydrocarbon gas pay zone characterization in a subterranean reservoir
US20170023689A1 (en) Predicting mechanical and elastic rock properties of the subsurface
CN104755960A (zh) 基于盆地建模对用于处理地震数据的速度模型的改进
WO2021167980A1 (en) Determination of calibrated minimum horizontal stress magnitude using fracture closure pressure and multiple mechanical earth model realizations
CA2695638A1 (en) Automated borehole image interpretation
Campbell et al. Natural and hydraulic fracture density prediction and identification of controllers
Abdideh et al. Cluster analysis of petrophysical and geological parameters for separating the electrofacies of a gas carbonate reservoir sequence
US20130292111A1 (en) Method of constructing a well log of a quantitative property from sample measurements and log data
Jorden et al. Well logging: rock properties, borehole environment, mud and temperature logging
US10151849B2 (en) Well placement using closure stress based landing map
US11215732B2 (en) Geological constraint using probability functions in stochastic mineralogy modeling
Duan et al. Architectural characterization of Ordovician fault-controlled paleokarst carbonate reservoirs, Tahe oilfield, China
Lai et al. Typical misinterpretations and scientific concepts in well-logging geologic studies
CN111123378B (zh) 确定划分岩性类型的伽马射线强度临界值的方法及装置
CN112528106A (zh) 火山岩岩性识别方法
Anxionnaz et al. Computer-generated corelike descriptions from open-hole logs (1)
US11905809B2 (en) Determining reservoir heterogeneity for optimized drilling location

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
GR01 Patent grant
GR01 Patent grant