CN107942400A - 一种预测钻井风险的方法及装置 - Google Patents

一种预测钻井风险的方法及装置 Download PDF

Info

Publication number
CN107942400A
CN107942400A CN201711103418.5A CN201711103418A CN107942400A CN 107942400 A CN107942400 A CN 107942400A CN 201711103418 A CN201711103418 A CN 201711103418A CN 107942400 A CN107942400 A CN 107942400A
Authority
CN
China
Prior art keywords
mtd
mrow
msub
mtr
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.)
Granted
Application number
CN201711103418.5A
Other languages
English (en)
Other versions
CN107942400B (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.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas 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 China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201711103418.5A priority Critical patent/CN107942400B/zh
Publication of CN107942400A publication Critical patent/CN107942400A/zh
Application granted granted Critical
Publication of CN107942400B publication Critical patent/CN107942400B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00—Prospecting or detecting by the use of nuclear radiation, e.g. of natural or induced radioactivity
    • G01V5/04—Prospecting or detecting by the use of nuclear radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
    • G01V5/08—Prospecting or detecting by the use of nuclear radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays
    • G01V5/12—Prospecting or detecting by the use of nuclear radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using gamma or X-ray sources

Abstract

本发明提供了一种预测钻井风险的方法及装置，该方法包括：获取待评价的储层内所有已钻井在水平面内的中心点的二维坐标值，并获取已钻井在多个预设深度的多个第一伽马值，相邻预设深度的深度差相等；将每一个已钻井的第一伽马值进行归一化处理，得到与已钻井的第一伽马值一一对应的第二伽马值；根据已钻井的第二伽马值、第二伽马值对应的预设深度以及已钻井的二维坐标值、储层内未钻井的单元区域的二维坐标值及预设深度，得到储层内未钻井的单元区域在多个预设深度的第二伽马值；计算单元区域的砂岩发育概率；根据砂岩发育概率，确定该单元区域的钻井风险大小。本发明提供的方法及装置，预测的储层内未钻井的区域的钻井风险准确率较高。

Description

一种预测钻井风险的方法及装置
技术领域
本发明涉及油气开采技术领域，特别涉及一种预测钻井风险的方法及装置。
背景技术
目前，在我国已投入开发的油田主要以陆相沉积为主。陆相沉积普遍存在储层厚度薄、横向变化快的特点。为了提高油田的开采效率，避免钻好的油井开采不出来油气，一般在钻井之前先进行风险评价。
目前主要通过地震资料进行风险评价，将地震资料的振幅平面图中地震波的振幅值较高的区域对应的区域确定为砂岩，振幅值较低的区域对应的位置确定为泥岩，确定砂岩为钻井风险小的区域，泥岩为钻井风险大的区域。
在实现本发明的过程中，本发明人发现现有技术中至少存在以下问题：
采用现有的方法进行钻井之前的风险评价时，仅根据地震资料的振幅平面图确定钻井的风险大小，预测的结果与实际钻井结果吻合程度不高，准确率较低。
发明内容
有鉴于此，本发明提供一种预测钻井风险的方法及装置，通过预测待钻井区域的砂岩发育概率的大小来评价钻井的风险。
具体而言，包括以下的技术方案：
第一方面，本发明提供了一种预测钻井风险的方法，包括：
获取待评价的储层内所有已钻井在水平面内的中心点的二维坐标值，并获取每一个所述已钻井在多个预设深度的多个第一伽马值，相邻的每两个预设深度之间的深度差相同；
将每一个所述已钻井的第一伽马值进行归一化处理，得到与所述已钻井的第一伽马值一一对应的第二伽马值；
根据所述已钻井的第二伽马值、所述已钻井的第二伽马值对应的预设深度以及所述已钻井的二维坐标值、所述储层内未钻井的单元区域的二维坐标值及所述多个预设深度，得到所述储层内未钻井的单元区域在所述多个预设深度的第二伽马值；
根据每一所述单元区域在所述多个预设深度的第二伽马值，计算得到该单元区域的砂岩发育概率；
根据所述单元区域的砂岩发育概率，确定该单元区域的钻井风险大小。
可选择地，所述将每一个所述已钻井的第一伽马值进行归一化处理，得到与所述已钻井的第一伽马值一一对应的所述已钻井的第二伽马值，采用的公式为：
式中，
GR*——已钻井在某一预设深度的第二伽马值；
GR——该已钻井在该预设深度的第一伽马值；
GRmax——该已钻井在多个预设深度的第一伽马值中的最大值；
GRmin——该已钻井在多个预设深度的第一伽马值中的最小值。
可选择地，所述根据所述已钻井的第二伽马值、所述已钻井的第二伽马值对应的预设深度以及所述已钻井的二维坐标值、所述储层内未钻井的单元区域的二维坐标值及所述多个预设深度，得到所述储层内未钻井的单元区域在所述多个预设深度的第二伽马值，包括以下步骤：
步骤1：将每一个所述已钻井在多个预设深度值处的多个位置作为第一位置，每一个所述第一位置的二维坐标值为其对应的已钻井的二维坐标值，深度值为对应的预设深度值；
步骤2：将所述储层内未钻井的单元区域在多个预设深度值处的多个位置作为第二位置，所述第二位置的二维坐标值为其对应的单元区域的二维坐标值，深度值与对应的预设深度值；
步骤3：根据第一位置的二维坐标值及深度值、所述第二位置的二维坐标值及深度值，计算所述第二位置与每一个所述第一位置之间的距离；
步骤4：将与所述第二位置之间的距离小于或等于预设距离的第一位置作为该第二位置的相关位置；
步骤5：根据所述第二位置的相关位置的第二伽马值计算该第二位置的第二伽马值；
步骤6：将计算得到第二伽马值的第二位置作为新的第一位置，重复步骤3至步骤5，计算其他第二位置的第二伽马值；
直到计算出所述储层内所有第二位置的第二伽马值。
可选择地，所述根据所述第二位置的相关位置的第二伽马值计算该第二位置的第二伽马值，采用的计算公式为：
式中，
GR*(a0)——第二位置a0的第二伽马值；
GR*(ai)——第二位置a0的相关位置ai的第二伽马值；
λi——权重系数；
n——第二位置a0的相关位置的个数；
其中，权重系数λi的计算公式为：
式中，
cij——第二位置a0的相关位置ai与aj的相关值，i＝1,2,…,n，j＝1,2,…,n；
ci0——第二位置a0与其相关位置ai的相关值，i＝1,2,…,n；
μ——拉格朗日乘数因子；
cij的计算公式如下：
当i＝j时，cij＝c11＝c22＝…＝c0+c
当i≠j时，cij＝c0+c-γ(|ai-aj|)
式中，c0和c为常数；
γ(|ai-aj|)的计算公式如下：
式中，
xi——位置ai的横坐标的值；
yi——位置ai的纵坐标的值；
zi——位置ai的深度值；
xj——位置aj的横坐标的值；
yj——位置aj的纵坐标的值；
zj——位置aj的深度值；
a——常数。
可选择地，所述根据每一所述单元区域在所述多个预设深度的第二伽马值，计算得到该单元区域的砂岩发育概率，采用的计算公式为：
式中，
η——某一单元区域内的砂岩发育概率；
m——该单元区域在多个预设深度的多个第二伽马值中小于0.4的第二伽马值的个数；
n——该单元区域在多个预设深度的多个第二伽马值的个数。
可选择地，所述根据所述单元区域内的砂岩发育概率，确定所述该单元区域的钻井风险大小，包括：
将砂岩发育概率大于80％的所述单元区域的钻井风险确定为低风险；
将砂岩发育概率在60％到80％之间的所述单元区域的钻井风险确定为较低风险；
将砂岩发育概率在40％到60％之间的所述单元区域的钻井风险确定为中度风险；
将砂岩发育概率在20％到40％之间的所述单元区域的钻井风险确定为较高风险；
将砂岩发育概率小于20％的所述单元区域的钻井风险确定为高风险。
第二方面，本发明还提供了一种预测钻井风险的装置，包括：
获取模块，用于获取待评价的储层内所有已钻井在水平面内的中心点的二维坐标值，并获取每一个所述已钻井在多个预设深度的多个第一伽马值，相邻的每两个预设深度之间的深度差相同；
处理模块，用于将每一个所述已钻井的第一伽马值进行归一化处理，得到与所述已钻井的第一伽马值一一对应的第二伽马值；
第一计算模块，用于根据所述已钻井的第二伽马值、所述已钻井的第二伽马值对应的预设深度以及所述已钻井的二维坐标值、所述储层内未钻井的单元区域的二维坐标值及所述多个预设深度，得到所述储层内未钻井的单元区域在所述多个预设深度的第二伽马值；
第二计算模块，用于根据每一所述单元区域在所述多个预设深度的第二伽马值，计算得到该单元区域的砂岩发育概率；
确定模块，用于根据所述单元区域的砂岩发育概率，确定该单元区域的钻井风险大小。
可选择地，所述处理模块采用的公式为：
式中，
GR*——已钻井在某一预设深度的第二伽马值；
GR——该已钻井在该预设深度的第一伽马值；
GRmax——该已钻井在多个预设深度的第一伽马值中的最大值；
GRmin——该已钻井在多个预设深度的第一伽马值中的最小值。
可选择地，所述第一计算模块的计算流程包括：
步骤1：将每一个所述已钻井在多个预设深度值处的多个位置作为第一位置，每一个所述第一位置的二维坐标值为其对应的已钻井的二维坐标值，深度值为对应的预设深度值；
步骤2：将所述储层内未钻井的单元区域在多个预设深度值处的多个位置作为第二位置，所述第二位置的二维坐标值为其对应的单元区域的二维坐标值，深度值与对应的预设深度值；
步骤3：根据第一位置的二维坐标值及深度值、所述第二位置的二维坐标值及深度值，计算所述第二位置与每一个所述第一位置之间的距离；
步骤4：将与所述第二位置之间的距离小于或等于预设距离的第一位置作为该第二位置的相关位置；
步骤5：根据所述第二位置的相关位置的第二伽马值计算该第二位置的第二伽马值；
步骤6：将计算得到第二伽马值的第二位置作为新的第一位置，重复步骤3至步骤5，计算其他第二位置的第二伽马值；
直到计算出所述储层内所有第二位置的第二伽马值。
可选择地，所述根据所述第二位置的相关位置的第二伽马值计算该第二位置的第二伽马值，采用的计算公式为：
式中，
GR*(a0)——第二位置a0的第二伽马值；
GR*(ai)——第二位置a0的相关位置ai的第二伽马值；
λi——权重系数；
n——第二位置a0的相关位置的个数；
其中，权重系数λi的计算公式为：
式中，
cij——第二位置a0的相关位置ai与aj的相关值，i＝1,2,…,n，j＝1,2,…,n；
ci0——第二位置a0与其相关位置ai的相关值，i＝1,2,…,n；
μ——拉格朗日乘数因子；
cij的计算公式如下：
当i＝j时，cij＝c11＝c22＝…＝c0+c
当i≠j时，cij＝c0+c-γ(|ai-aj|)
式中，c0和c为常数；
γ(|ai-aj|)的计算公式如下：
式中，
xi——位置ai的横坐标的值；
yi——位置ai的纵坐标的值；
zi——位置ai的深度值；
xj——位置aj的横坐标的值；
yj——位置aj的纵坐标的值；
zj——位置aj的深度值；
a——常数。
可选择地，所述第二计算模块采用的计算公式为：
式中，
η——某一单元区域内的砂岩发育概率；
m——该单元区域在多个预设深度的多个第二伽马值中小于0.4的第二伽马值的个数；
n——该单元区域在多个预设深度的多个第二伽马值的个数。
可选择地，所述确定模块用于：
将砂岩发育概率大于80％的所述单元区域的钻井风险确定为低风险；
将砂岩发育概率在60％到80％之间的所述单元区域的钻井风险确定为较低风险；
将砂岩发育概率在40％到60％之间的所述单元区域的钻井风险确定为中度风险；
将砂岩发育概率在20％到40％之间的所述单元区域的钻井风险确定为较高风险；
将砂岩发育概率小于20％的所述单元区域的钻井风险确定为高风险。
本发明实施例提供的技术方案的有益效果：
本发明实施例提供的预测钻井风险的方法及装置，通过对储层内的所有已钻井的预设深度处的第一伽玛值进行归一化处理，得到第二伽马值，从而可解决测井过中的人为和仪器等因素带来的误差；根据所有已钻井的第二伽马值，进行计算得到储层内未钻井的单元区域的砂岩发育概率；根据不同单元区域的砂岩发育概率，确定不同的单元区域的钻井风险大小。因此，本发明提供的方法及装置，是根据所有已钻井的第二伽马值对未钻井的区域的钻井风险进行预测，准确率较高。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种预测钻井风险的方法的流程图；
图2为本发明实施例提供的一种预测钻井风险的装置的框图。
具体实施方式
为使本发明的技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
本发明实施例提供了一种预测钻井风险的方法，如图1所示，包括步骤101、102、103、104和105。下面将对各步骤进行具体说明。
步骤101：获取待评价的储层内的已钻井在水平面内的中心点的二维坐标值，并获取每一个已钻井在多个预设深度的多个第一伽马值，相邻的每两个预设深度之间的深度差相同。
在储层内钻井时，获取每一个已钻井的自然伽玛测井曲线。自然伽玛测井曲线是伽玛值关于深度的曲线。根据自然伽玛曲线及多个预设深度，可以确定与多个预设深度一一对应的多个伽玛值，即为第一伽马值。
获取已钻井在水平面内的中心点的二维坐标值，用来后续计算储层内其他未钻井的区域的第二伽马值。
为了保证对第一伽马值的优化过程中的准确性，在确定预设深度时，确定相邻的两个预设深度之间的深度差相同。
具体地，预设深度及相邻的两个预设深度之间的深度差可根据储层的已钻井的自然伽玛测井曲线确定。
步骤102：将每一个已钻井的第一伽马值进行归一化处理，得到与已钻井的第一伽马值一一对应的第二伽马值。
由于已钻井条件的差异、测井所有的仪器刻度、测量的时间及测量的人员不同等因素，都会对已钻井的自然伽玛测井曲线造成影响。为了减小这些因素对后续计算储层内其他区域的砂岩发育概率的影响，先对每一个已钻井的第一伽马值进行归一化处理。归一化处理采用的公式如下：
式中，
GR*——已钻井在某一预设深度的第二伽马值；
GR——该已钻井在该预设深度的第一伽马值；
GRmax——该已钻井在多个预设深度的第一伽马值中的最大值；
GRmin——该已钻井在多个预设深度的第一伽马值中的最小值。
步骤103：根据已钻井的第二伽马值、已钻井的第二伽马值对应的预设深度以及已钻井的二维坐标值、储层内未钻井的单元区域的二维坐标值及多个预设深度，得到储层内未钻井的单元区域在多个预设深度的第二伽马值。
该步骤可通过子步骤1-6来完成。
步骤1：将每一个已钻井在多个预设深度值处的多个位置作为第一位置，每一个第一位置的二维坐标值为其对应的已钻井的二维坐标值，深度值为对应的预设深度值。
步骤2：将储层内未钻井的单元区域在多个预设深度值处的多个位置作为第二位置，第二位置的二维坐标值为其对应的单元区域的二维坐标值，深度值为对应的预设深度值。
步骤3：根据第一位置的二维坐标值及深度值、第二位置的二维坐标值及深度值，计算第二位置与每一个第一位置之间的距离。
步骤4：将与第二位置之间的距离小于或等于预设距离的第一位置作为该第二位置的相关位置。
步骤5：根据第二位置的相关位置的第二伽马值计算该第二位置的第二伽马值。
在该步骤中采用的计算公式为：
式中，
GR*(a0)——第二位置a0的第二伽马值；
GR*(ai)——第二位置a0的相关位置ai的第二伽马值；
λi——权重系数；
n——第二位置a0的相关位置的个数；
其中，权重系数λi的计算公式为：
式中，
cij——第二位置a0的相关位置ai与aj的相关值，i＝1,2,…,n，j＝1,2,…,n；
ci0——第二位置a0与其相关位置ai的相关值，i＝1,2,…,n；
μ——拉格朗日乘数因子；
cij的计算公式如下：
当i＝j时，cij＝c11＝c22＝…＝c0+c
当i≠j时，cij＝c0+c-γ(|ai-aj|)
式中，c0和c为常数；
γ(|ai-aj|)的计算公式如下：
式中，
xi——位置ai的横坐标的值；
yi——位置ai的纵坐标的值；
zi——位置ai的深度值；
xj——位置aj的横坐标的值；
yj——位置aj的纵坐标的值；
zj——位置aj的深度值；
a——常数。
步骤6：将计算得到第二伽马值的第二位置作为新的第一位置，重复步骤3至步骤5，计算其他第二位置的第二伽马值。
直到计算出所述储层内所有第二位置的第二伽马值。
其中，单元区域的大小及预设距离的值可根据储层的实际情况确定。一般，单元区域的面积越小，得到的钻井风险就越精确。比如，可取单元区域的面积为1m2。
步骤104：根据每一单元区域在多个预设深度的第二伽马值，计算得到该单元区域的砂岩发育概率。
计算单元区域的发育概率的计算公式如下：
采用的计算公式为：
式中，
η——某一单元区域内的砂岩发育概率；
m——该单元区域在多个预设深度的多个第二伽马值中小于0.4的第二伽马值的个数；
n——该单元区域在多个预设深度的多个第二伽马值的个数。
根据自然伽玛测井对岩性的划分原理，在沉积岩中含泥质越多，其放射性就越强，归一化后砂岩的第二伽马值就越大。通过对已钻井的自然伽玛测井曲线分析可以确定，第二伽玛值小于0.4的沉积岩为砂岩，第二伽马值大于0.4的沉积岩为泥岩。
根据单元区域的不同深度处的第二伽马值来计算得到单元区域的砂岩发育概率。
步骤105：根据单元区域的砂岩发育概率，确定该单元区域的钻井风险大小。
具体地，将砂岩发育概率大于80％的单元区域的钻井风险确定为低风险；将砂岩发育概率在60％到80％之间的单元区域的钻井风险确定为较低风险；将砂岩发育概率在40％到60％之间的单元区域的钻井风险确定为中度风险；将砂岩发育概率在20％到40％之间的单元区域的钻井风险确定为较高风险；将砂岩发育概率小于20％的单元区域的钻井风险确定为高风险。
得到不同单元区域的钻井风险的情况后，根据不同单元区域的钻井风险大小在储层内确定钻井的位置。
本发明实施例提供的预测钻井风险的方法，通过对储层内的所有已钻井的预设深度处的第一伽玛值进行归一化处理，得到第二伽马值，从而可解决测井过中的人为和仪器等因素带来的误差；根据所有已钻井的第二伽马值，进行计算得到储层内未钻井的单元区域的砂岩发育概率；根据不同单元区域的砂岩发育概率，确定不同的单元区域的钻井风险大小。因此，本发明提供的方法，根据所有已钻井的第二伽马值对未钻井的区域的钻井风险进行预测，准确率较高。
本发明实施例还提供概率一种预测钻井风险的装置，如图2所示，包括：
获取模块201，用于获取待评价的储层内所有已钻井在水平面内的中心点的二维坐标值，并获取每一个已钻井在多个预设深度的多个第一伽马值，相邻的每两个预设深度之间的深度差相同。
处理模块202，用于将每一个已钻井的第一伽马值进行归一化处理，得到与已钻井的第一伽马值一一对应的第二伽马值。
第一计算模块203，用于根据已钻井的第二伽马值、已钻井的第二伽马值对应的预设深度以及已钻井的二维坐标值、储层内未钻井的单元区域的二维坐标值及多个预设深度，得到储层内未钻井的单元区域在多个预设深度的第二伽马值。
第二计算模块204，用于根据每一单元区域在多个预设深度的第二伽马值，计算得到该单元区域的砂岩发育概率。
确定模块205，用于根据单元区域的砂岩发育概率，确定该单元区域的钻井风险大小。
可选择地，处理模块202采用的公式为：
式中，
GR*——已钻井在某一预设深度的第二伽马值；
GR——该已钻井在该预设深度的第一伽马值；
GRmax——该已钻井在多个预设深度的第一伽马值中的最大值；
GRmin——该已钻井在多个预设深度的第一伽马值中的最小值。
可选择地，第一计算模块203的计算流程包括：
步骤1：将每一个已钻井在多个预设深度值处的多个位置作为第一位置，每一个第一位置的二维坐标值为其对应的已钻井的二维坐标值，深度值为对应的预设深度值；
步骤2：将储层内未钻井的单元区域在多个预设深度值处的多个位置作为第二位置，第二位置的二维坐标值为其对应的单元区域的二维坐标值，深度值与对应的预设深度值；
步骤3：根据第一位置的二维坐标值及深度值、第二位置的二维坐标值及深度值，计算第二位置与每一个第一位置之间的距离；
步骤4：将与第二位置之间的距离小于或等于预设距离的第一位置作为该第二位置的相关位置；
步骤5：根据第二位置的相关位置的第二伽马值计算该第二位置的第二伽马值；
步骤6：将计算得到第二伽马值的第二位置作为新的第一位置，重复步骤3至步骤5，计算其他第二位置的第二伽马值；
直到计算出储层内所有第二位置的第二伽马值。
可选择地，根据所述第二位置的相关位置的第二伽马值计算该第二位置的第二伽马值，采用的计算公式为：
式中，
GR*(a0)——第二位置a0的第二伽马值；
GR*(ai)——第二位置a0的相关位置ai的第二伽马值；
λi——权重系数；
n——第二位置a0的相关位置的个数；
其中，权重系数λi的计算公式为：
式中，
cij——第二位置a0的相关位置ai与aj的相关值，i＝1,2,…,n，j＝1,2,…,n；
ci0——第二位置a0与其相关位置ai的相关值，i＝1,2,…,n；
μ——拉格朗日乘数因子；
cij的计算公式如下：
当i＝j时，cij＝c11＝c22＝…＝c0+c
当i≠j时，cij＝c0+c-γ(|ai-aj|)
式中，c0和c为常数；
γ(|ai-aj|)的计算公式如下：
式中，
xi——位置ai的横坐标的值；
yi——位置ai的纵坐标的值；
zi——位置ai的深度值；
xj——位置aj的横坐标的值；
yj——位置aj的纵坐标的值；
zj——位置aj的深度值；
a——常数。
可选择地，第二计算模块204采用的计算公式为：
式中，
η——某一单元区域内的砂岩发育概率；
m——该单元区域在多个预设深度的多个第二伽马值中小于0.4的第二伽马值的个数；
n——该单元区域在多个预设深度的多个第二伽马值的个数。
可选择地，确定模块205用于：
将砂岩发育概率大于80％的单元区域的钻井风险确定为低风险；
将砂岩发育概率在60％到80％之间的单元区域的钻井风险确定为较低风险；
将砂岩发育概率在40％到60％之间的单元区域的钻井风险确定为中度风险；
将砂岩发育概率在20％到40％之间的单元区域的钻井风险确定为较高风险；
将砂岩发育概率小于20％的单元区域的钻井风险确定为高风险。
本发明实施例提供的预测钻井风险的装置，通过对储层内的已钻井的预设深度处的第一伽玛值进行归一化处理，得到第二伽马值，从而可解决测井过中的人为和仪器等因素带来的误差；根据已钻井的第二伽马值，进行计算得到储层内未钻井的单元区域的砂岩发育概率；根据不同单元区域的砂岩发育概率，确定不同的单元区域的钻井风险大小。因此，本发明提供的方法，根据所有已钻井的第二伽马值对未钻井的区域的钻井风险进行预测，准确率较高。
以上所述仅是为了便于本领域的技术人员理解本发明的技术方案，并不用以限制本发明。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (12)

1.一种预测钻井风险的方法，其特征在于，包括：
获取待评价的储层内所有已钻井在水平面内的中心点的二维坐标值，并获取每一个所述已钻井在多个预设深度的多个第一伽马值，相邻的每两个预设深度之间的深度差相同；
将每一个所述已钻井的第一伽马值进行归一化处理，得到与所述已钻井的第一伽马值一一对应的第二伽马值；
根据所述已钻井的第二伽马值、所述已钻井的第二伽马值对应的预设深度以及所述已钻井的二维坐标值、所述储层内未钻井的单元区域的二维坐标值及所述多个预设深度，得到所述储层内未钻井的单元区域在所述多个预设深度的第二伽马值；
根据每一所述单元区域在所述多个预设深度的第二伽马值，计算得到该单元区域的砂岩发育概率；
根据所述单元区域的砂岩发育概率，确定该单元区域的钻井风险大小。
2.根据权利要求1所述的预测钻井风险的方法，其特征在于，所述将每一个所述已钻井的第一伽马值进行归一化处理，得到与所述已钻井的第一伽马值一一对应的所述已钻井的第二伽马值，采用的公式为：
<mrow> <msup> <mi>GR</mi> <mo>*</mo> </msup> <mo>=</mo> <mfrac> <mrow> <mi>G</mi> <mi>R</mi> <mo>-</mo> <msub> <mi>GR</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>GR</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>GR</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> </mfrac> </mrow>
式中，
GR*——已钻井在某一预设深度的第二伽马值；
GR——该已钻井在该预设深度的第一伽马值；
GRmax——该已钻井在多个预设深度的第一伽马值中的最大值；
GRmin——该已钻井在多个预设深度的第一伽马值中的最小值。
3.根据权利要求1所述的预测钻井风险的方法，其特征在于，所述根据所述已钻井的第二伽马值、所述已钻井的第二伽马值对应的预设深度以及所述已钻井的二维坐标值、所述储层内未钻井的单元区域的二维坐标值及所述多个预设深度，得到所述储层内未钻井的单元区域在所述多个预设深度的第二伽马值，包括以下步骤：
步骤1：将每一个所述已钻井在多个预设深度值处的多个位置作为第一位置，每一个所述第一位置的二维坐标值为其对应的已钻井的二维坐标值，深度值为对应的预设深度值；
步骤2：将所述储层内未钻井的单元区域在多个预设深度值处的多个位置作为第二位置，所述第二位置的二维坐标值为其对应的单元区域的二维坐标值，深度值与对应的预设深度值；
步骤3：根据第一位置的二维坐标值及深度值、所述第二位置的二维坐标值及深度值，计算所述第二位置与每一个所述第一位置之间的距离；
步骤4：将与所述第二位置之间的距离小于或等于预设距离的第一位置作为该第二位置的相关位置；
步骤5：根据所述第二位置的相关位置的第二伽马值计算该第二位置的第二伽马值；
步骤6：将计算得到第二伽马值的第二位置作为新的第一位置，重复步骤3至步骤5，计算其他第二位置的第二伽马值；
直到计算出所述储层内所有第二位置的第二伽马值。
4.根据权利要求3所述的预测钻井风险的方法，其特征在于，所述根据所述第二位置的相关位置的第二伽马值计算该第二位置的第二伽马值，采用的计算公式为：
<mrow> <msup> <mi>GR</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>&amp;lambda;</mi> <mi>i</mi> </msub> <msup> <mi>GR</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mrow> <mo>(</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow>
式中，
GR*(a0)——第二位置a0的第二伽马值；
GR*(ai)——第二位置a0的相关位置ai的第二伽马值；
λi——权重系数；
n——第二位置a0的相关位置的个数；
其中，权重系数λi的计算公式为：
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>&amp;lambda;</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;lambda;</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;lambda;</mi> <mi>n</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>&amp;mu;</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>*</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>c</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>12</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>c</mi> <mrow> <mn>1</mn> <mi>n</mi> </mrow> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>22</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>c</mi> <mrow> <mn>2</mn> <mi>n</mi> </mrow> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mrow> <mi>n</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mrow> <mi>n</mi> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>c</mi> <mrow> <mi>n</mi> <mi>n</mi> </mrow> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>c</mi> <mn>10</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mn>20</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mrow> <mi>n</mi> <mn>0</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中，
cij——第二位置a0的相关位置ai与aj的相关值，i＝1,2,…,n，j＝1,2,…,n；
ci0——第二位置a0与其相关位置ai的相关值，i＝1,2,…,n；
μ——拉格朗日乘数因子；
cij的计算公式如下：
当i＝j时，cij＝c11＝c22＝…＝c0+c
当i≠j时，cij＝c0+c-γ(|ai-aj|)
式中，c0和c为常数；
γ(|ai-aj|)的计算公式如下：
<mrow> <mi>h</mi> <mo>=</mo> <mo>|</mo> <msub> <mi>a</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>a</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>z</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow>
<mrow> <mi>&amp;gamma;</mi> <mrow> <mo>(</mo> <mi>h</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>h</mi> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>c</mi> <mn>0</mn> </msub> <mo>+</mo> <mi>c</mi> <mo>*</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>3</mn> <mi>h</mi> </mrow> <mrow> <mn>2</mn> <mi>a</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <msup> <mi>h</mi> <mn>3</mn> </msup> <mrow> <mn>2</mn> <msup> <mi>a</mi> <mn>3</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mn>0</mn> <mo>&lt;</mo> <mi>h</mi> <mo>&amp;le;</mo> <mi>a</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>c</mi> <mn>0</mn> </msub> <mo>+</mo> <mi>c</mi> </mrow> </mtd> <mtd> <mrow> <mi>h</mi> <mo>&gt;</mo> <mi>a</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中，
xi——位置ai的横坐标的值；
yi——位置ai的纵坐标的值；
zi——位置ai的深度值；
xj——位置aj的横坐标的值；
yj——位置aj的纵坐标的值；
zj——位置aj的深度值；
a——常数。
5.根据权利要求1所述的预测钻井风险的方法，其特征在于，所述根据每一所述单元区域在所述多个预设深度的第二伽马值，计算得到该单元区域的砂岩发育概率，采用的计算公式为：
<mrow> <mi>&amp;eta;</mi> <mo>=</mo> <mfrac> <mi>m</mi> <mi>n</mi> </mfrac> </mrow>
式中，
η——某一单元区域内的砂岩发育概率；
m——该单元区域在多个预设深度的多个第二伽马值中小于0.4的第二伽马值的个数；
n——该单元区域在多个预设深度的多个第二伽马值的个数。
6.根据权利要求1-5任一项所述的预测钻井风险的方法，其特征在于，所述根据所述单元区域内的砂岩发育概率，确定所述该单元区域的钻井风险大小，包括：
将砂岩发育概率大于80％的所述单元区域的钻井风险确定为低风险；
将砂岩发育概率在60％到80％之间的所述单元区域的钻井风险确定为较低风险；
将砂岩发育概率在40％到60％之间的所述单元区域的钻井风险确定为中度风险；
将砂岩发育概率在20％到40％之间的所述单元区域的钻井风险确定为较高风险；
将砂岩发育概率小于20％的所述单元区域的钻井风险确定为高风险。
7.一种预测钻井风险的装置，其特征在于，包括：
获取模块，用于获取待评价的储层内所有已钻井在水平面内的中心点的二维坐标值，并获取每一个所述已钻井在多个预设深度的多个第一伽马值，相邻的每两个预设深度之间的深度差相同；
处理模块，用于将每一个所述已钻井的第一伽马值进行归一化处理，得到与所述已钻井的第一伽马值一一对应的第二伽马值；
第一计算模块，用于根据所述已钻井的第二伽马值、所述已钻井的第二伽马值对应的预设深度以及所述已钻井的二维坐标值、所述储层内未钻井的单元区域的二维坐标值及所述多个预设深度，得到所述储层内未钻井的单元区域在所述多个预设深度的第二伽马值；
第二计算模块，用于根据每一所述单元区域在所述多个预设深度的第二伽马值，计算得到该单元区域的砂岩发育概率；
确定模块，用于根据所述单元区域的砂岩发育概率，确定该单元区域的钻井风险大小。
8.根据权利要求7所述的预测钻井风险的装置，其特征在于，所述处理模块采用的公式为：
<mrow> <msup> <mi>GR</mi> <mo>*</mo> </msup> <mo>=</mo> <mfrac> <mrow> <mi>G</mi> <mi>R</mi> <mo>-</mo> <msub> <mi>GR</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>GR</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>GR</mi> <mrow> <mi>m</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> </mfrac> </mrow>
式中，
GR*——已钻井在某一预设深度的第二伽马值；
GR——该已钻井在该预设深度的第一伽马值；
GRmax——该已钻井在多个预设深度的第一伽马值中的最大值；
GRmin——该已钻井在多个预设深度的第一伽马值中的最小值。
9.根据权利要求7所述的预测钻井风险的装置，其特征在于，所述第一计算模块的计算流程包括：
步骤1：将每一个所述已钻井在多个预设深度值处的多个位置作为第一位置，每一个所述第一位置的二维坐标值为其对应的已钻井的二维坐标值，深度值为对应的预设深度值；
步骤2：将所述储层内未钻井的单元区域在多个预设深度值处的多个位置作为第二位置，所述第二位置的二维坐标值为其对应的单元区域的二维坐标值，深度值与对应的预设深度值；
步骤3：根据第一位置的二维坐标值及深度值、所述第二位置的二维坐标值及深度值，计算所述第二位置与每一个所述第一位置之间的距离；
步骤4：将与所述第二位置之间的距离小于或等于预设距离的第一位置作为该第二位置的相关位置；
步骤5：根据所述第二位置的相关位置的第二伽马值计算该第二位置的第二伽马值；
步骤6：将计算得到第二伽马值的第二位置作为新的第一位置，重复步骤3至步骤5，计算其他第二位置的第二伽马值；
直到计算出所述储层内所有第二位置的第二伽马值。
10.根据权利要求9所述的预测钻井风险的装置，其特征在于，所述根据所述第二位置的相关位置的第二伽马值计算该第二位置的第二伽马值，采用的计算公式为：
<mrow> <msup> <mi>GR</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <msub> <mi>&amp;lambda;</mi> <mi>i</mi> </msub> <msup> <mi>GR</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <msub> <mi>a</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> <mo>,</mo> <mrow> <mo>(</mo> <mi>i</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow>
式中，
GR*(a0)——第二位置a0的第二伽马值；
GR*(ai)——第二位置a0的相关位置ai的第二伽马值；
λi——权重系数；
n——第二位置a0的相关位置的个数；
其中，权重系数λi的计算公式为：
<mrow> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>&amp;lambda;</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;lambda;</mi> <mn>2</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&amp;lambda;</mi> <mi>n</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>&amp;mu;</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>*</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>c</mi> <mn>11</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>12</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>c</mi> <mrow> <mn>1</mn> <mi>n</mi> </mrow> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mn>21</mn> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mn>22</mn> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>c</mi> <mrow> <mn>2</mn> <mi>n</mi> </mrow> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mrow> <mi>n</mi> <mn>1</mn> </mrow> </msub> </mtd> <mtd> <msub> <mi>c</mi> <mrow> <mi>n</mi> <mn>2</mn> </mrow> </msub> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <msub> <mi>c</mi> <mrow> <mi>n</mi> <mi>n</mi> </mrow> </msub> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>...</mn> </mtd> <mtd> <mn>1</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>c</mi> <mn>10</mn> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mn>20</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <msub> <mi>c</mi> <mrow> <mi>n</mi> <mn>0</mn> </mrow> </msub> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中，
cij——第二位置a0的相关位置ai与aj的相关值，i＝1,2,…,n，j＝1,2,…,n；
ci0——第二位置a0与其相关位置ai的相关值，i＝1,2,…,n；
μ——拉格朗日乘数因子；
cij的计算公式如下：
当i＝j时，cij＝c11＝c22＝…＝c0+c
当i≠j时，cij＝c0+c-γ(|ai-aj|)
式中，c0和c为常数；
γ(|ai-aj|)的计算公式如下：
<mrow> <mi>h</mi> <mo>=</mo> <mo>|</mo> <msub> <mi>a</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>a</mi> <mi>j</mi> </msub> <mo>|</mo> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>x</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>y</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>z</mi> <mi>j</mi> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> </mrow>
<mrow> <mi>&amp;gamma;</mi> <mrow> <mo>(</mo> <mi>h</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mrow> <mi>h</mi> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>c</mi> <mn>0</mn> </msub> <mo>+</mo> <mi>c</mi> <mo>*</mo> <mrow> <mo>(</mo> <mfrac> <mrow> <mn>3</mn> <mi>h</mi> </mrow> <mrow> <mn>2</mn> <mi>a</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <msup> <mi>h</mi> <mn>3</mn> </msup> <mrow> <mn>2</mn> <msup> <mi>a</mi> <mn>3</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow> </mtd> <mtd> <mrow> <mn>0</mn> <mo>&lt;</mo> <mi>h</mi> <mo>&amp;le;</mo> <mi>a</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>c</mi> <mn>0</mn> </msub> <mo>+</mo> <mi>c</mi> </mrow> </mtd> <mtd> <mrow> <mi>h</mi> <mo>&gt;</mo> <mi>a</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中，
xi——位置ai的横坐标的值；
yi——位置ai的纵坐标的值；
zi——位置ai的深度值；
xj——位置aj的横坐标的值；
yj——位置aj的纵坐标的值；
zj——位置aj的深度值；
a——常数。
11.根据权利要求7所述的预测钻井风险的装置，其特征在于，所述第二计算模块采用的计算公式为：
<mrow> <mi>&amp;eta;</mi> <mo>=</mo> <mfrac> <mi>m</mi> <mi>n</mi> </mfrac> </mrow>
式中，
η——某一单元区域内的砂岩发育概率；
m——该单元区域在多个预设深度的多个第二伽马值中小于0.4的第二伽马值的个数；
n——该单元区域在多个预设深度的多个第二伽马值的个数。
12.根据权利要求7-11任一项所述的预测钻井风险的装置，其特征在于，所述确定模块用于：
将砂岩发育概率大于80％的所述单元区域的钻井风险确定为低风险；
将砂岩发育概率在60％到80％之间的所述单元区域的钻井风险确定为较低风险；
将砂岩发育概率在40％到60％之间的所述单元区域的钻井风险确定为中度风险；
将砂岩发育概率在20％到40％之间的所述单元区域的钻井风险确定为较高风险；
将砂岩发育概率小于20％的所述单元区域的钻井风险确定为高风险。
CN201711103418.5A 2017-11-10 2017-11-10 一种预测钻井风险的方法及装置 Active CN107942400B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711103418.5A CN107942400B (zh) 2017-11-10 2017-11-10 一种预测钻井风险的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711103418.5A CN107942400B (zh) 2017-11-10 2017-11-10 一种预测钻井风险的方法及装置

Publications (2)

Publication Number Publication Date
CN107942400A true CN107942400A (zh) 2018-04-20
CN107942400B CN107942400B (zh) 2019-06-11

Family

ID=61933718

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711103418.5A Active CN107942400B (zh) 2017-11-10 2017-11-10 一种预测钻井风险的方法及装置

Country Status (1)

Country Link
CN (1) CN107942400B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112302640A (zh) * 2020-10-28 2021-02-02 中国石油天然气集团有限公司 水平井压裂工程风险预测方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104297787A (zh) * 2014-10-17 2015-01-21 中国石油天然气股份有限公司 河流相低渗致密砂岩储层的三维岩相数据处理方法及装置
CN106094052A (zh) * 2016-06-01 2016-11-09 中国地质大学（武汉） 一种致密白云岩储层的裂缝发育程度识别方法
US9556728B2 (en) * 2014-01-13 2017-01-31 Varel Europe S.A.S. Methods and systems of analyzing wellbore drilling operations
CN106991279A (zh) * 2017-03-29 2017-07-28 中国石油大学(北京) 一种致密砂岩气成藏概率的确定方法和装置
CN107045154A (zh) * 2017-02-08 2017-08-15 中国海洋石油总公司 一种水平井环境中的识别地层产状的方法和装置
CN107329171A (zh) * 2017-06-07 2017-11-07 中国石油天然气股份有限公司 深度域储层地震反演方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9556728B2 (en) * 2014-01-13 2017-01-31 Varel Europe S.A.S. Methods and systems of analyzing wellbore drilling operations
CN104297787A (zh) * 2014-10-17 2015-01-21 中国石油天然气股份有限公司 河流相低渗致密砂岩储层的三维岩相数据处理方法及装置
CN106094052A (zh) * 2016-06-01 2016-11-09 中国地质大学（武汉） 一种致密白云岩储层的裂缝发育程度识别方法
CN107045154A (zh) * 2017-02-08 2017-08-15 中国海洋石油总公司 一种水平井环境中的识别地层产状的方法和装置
CN106991279A (zh) * 2017-03-29 2017-07-28 中国石油大学(北京) 一种致密砂岩气成藏概率的确定方法和装置
CN107329171A (zh) * 2017-06-07 2017-11-07 中国石油天然气股份有限公司 深度域储层地震反演方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
方立亭: "伊拉克鲁迈拉油田岩性特征及钻井风险评价", 《西部探矿工程》 *
郭智 等: "辫状河相致密砂岩气藏三维地质建模", 《石油勘探与开发》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112302640A (zh) * 2020-10-28 2021-02-02 中国石油天然气集团有限公司 水平井压裂工程风险预测方法及装置

Also Published As

Publication number Publication date
CN107942400B (zh) 2019-06-11

Similar Documents

Publication Publication Date Title
US20210332690A1 (en) Method and system of combined support for a well drilling process
RU2617877C2 (ru) Способ и устройство для оптимизации глубинных измерений удельного сопротивления с использованием многокомпонентных антенн
US20120312090A1 (en) Method for evaluation of hydrocarbon content of shale
CN103114840A (zh) 一种高-过高成熟页岩有机碳含量计算方法及装置
Santos et al. Characterization of natural fracture systems: Analysis of uncertainty effects in linear scanline results
CN103774988A (zh) 一种实时随钻伽马正演地质导向钻井方法
CN104213899A (zh) 一种地层岩石骨架的测井识别方法
CN107784159A (zh) 一种储层电阻率各向异性系数的确定方法
CN106154322A (zh) 测井曲线校正方法和装置
CN107807410A (zh) 储层地球物理敏感参数优选方法及系统
US11008835B2 (en) Horizontal reservoir description systems
CN103615230B (zh) 一种双泥质指示因子含水饱和度模型的建立方法
US10317555B2 (en) Method of minimizing tool response for downhole logging operations
Liu et al. An improved probability combination scheme based on principal component analysis and permanence of ratios model-An application to a fractured reservoir modeling, Ordos Basin
CN103527172B (zh) 可变岩电耦合指数含水饱和度计算方法
CN104459800A (zh) 一种预测砂体尖灭的方法和装置
US11249218B2 (en) Horizontal reservoir description systems
CN107942400B (zh) 一种预测钻井风险的方法及装置
CN106522919A (zh) 三探头密度测井的围岩校正方法及装置
CN106226217A (zh) 一种确定地下裂缝张开度的方法、装置及其应用
CN110261467A (zh) 一种识别碳酸盐岩古岩溶储层垂向分带性的方法
CN116930023A (zh) 致密砂岩相控分类孔隙度测井精细解释方法和装置
CN105242317B (zh) 一种纵波速度的确定方法及装置
Freeman et al. Fault seal mapping–incorporating geometric and property uncertainty
Lemieux et al. Fracture network analysis of the St-Eustache quarry, Quebec, Canada, for groundwater resources management

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