CN108961352A - 一种随钻测井曲线的绘制方法 - Google Patents

一种随钻测井曲线的绘制方法 Download PDF

Info

Publication number
CN108961352A
CN108961352A CN201710356097.3A CN201710356097A CN108961352A CN 108961352 A CN108961352 A CN 108961352A CN 201710356097 A CN201710356097 A CN 201710356097A CN 108961352 A CN108961352 A CN 108961352A
Authority
CN
China
Prior art keywords
point
vertical depth
drafting
well logging
sequence
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
CN201710356097.3A
Other languages
English (en)
Other versions
CN108961352B (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 Chemical Corp
Sinopec Research Institute of Petroleum Engineering
Original Assignee
China Petroleum and Chemical Corp
Sinopec Research Institute of Petroleum Engineering
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 Chemical Corp, Sinopec Research Institute of Petroleum Engineering filed Critical China Petroleum and Chemical Corp
Priority to CN201710356097.3A priority Critical patent/CN108961352B/zh
Publication of CN108961352A publication Critical patent/CN108961352A/zh
Application granted granted Critical
Publication of CN108961352B publication Critical patent/CN108961352B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/20Drawing from basic elements, e.g. lines or circles
    • G06T11/203Drawing of straight lines or curves

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种随钻测井曲线的绘制方法,属于油气开发与勘探领域。该方法包括:第一步,建立纵向地层剖面;第二步,将井眼轨迹数据投影到纵向地层剖面上;第三步,沿井眼轨迹构建二维的垂深‑测井响应量数据对序列;第四步,沿井眼轨迹构建二维的垂深‑水平位移对序列;第五步,建立曲线绘制坐标绘制道,并绘制二维的垂深‑测井响应量数据对序列中的第一点N0;第六步,沿井眼轨迹,以当前绘制点之后连续M个二维的垂深‑水平位移对序列中的点为基础,采用三次样条算法,形成回归方程;第七步,对回归方程进行二阶求导,然后搜索拐点,从拐点处分道绘制测井曲线;第八步,结束。利用本发明方法实现了更加清晰的垂向的地层特征显示。

Description

一种随钻测井曲线的绘制方法
技术领域
本发明属于油气开发与勘探领域,具体涉及一种随钻测井曲线的绘制方法,实现更加清晰的垂向的地层特征显示。
背景技术
石油测井曲线是直观分析地层岩性、物性和流体性质的主要手段。测井资料一直是按照垂直深度以曲线形式显示测井数值;即规定井口位置为零点、以垂直深度和某深度点的测井物理量幅度值为直角坐标,通过纵横比例变换后确定其在显示介质上的位置、并逐点连线形成测井曲线。在直井中这种显示方法可根据测井曲线的幅度变化真实反映地层某个物理参数的变化。
目前,在石油、天然气等钻井勘探开发技术领域,在大斜度井、水平井作业应用的越来越广泛。在大斜度井和水平井普遍使用随钻测井方式,测井的曲线不适合再使用传统直井的测井曲线的绘制方式。若采用沿斜深绘制测井曲线的传统方式显示大斜度井、水平井测井资料,使得斜井中的地层视厚度远大于真厚度、不能确定地层真实厚度,也不能反映井眼轨迹在水平方向的位移大小、直观表达井眼轨迹与地层界面的空间接触关系。常规的做法是将水平井的测量曲线投影到垂向线上。然而,对于测井轨迹是在上下地层中来回反复穿行情况下,则投影后的曲线是重复的,人们无法分辨出地层的特征,垂向投影的方法失去了作用。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种随钻测井曲线的绘制方法,克服对于测井轨迹是在上下地层中来回反复穿行情况下,则投影后的曲线是重复的,无法分辨出地层的特征,垂向投影的方法失去了作用的缺点,从而指导钻井人员快速准确调整钻井轨迹倾角,预防打穿目标层,确保轨迹在目标层内穿行。并为钻井轨迹、钻井参数调整提提供决策依据,从而提高水平井钻井效率、优质储层钻遇率。
本发明是通过以下技术方案实现的:
一种随钻测井曲线的绘制方法,包括:
第一步,建立纵向地层剖面;
第二步,将井眼轨迹数据投影到纵向地层剖面上;
第三步,沿井眼轨迹构建二维的垂深-测井响应量数据对序列;
第四步,沿井眼轨迹构建二维的垂深-水平位移对序列;
第五步,建立曲线绘制坐标绘制道,并绘制二维的垂深-测井响应量数据对序列中的第一点N0;
第六步,沿井眼轨迹,以当前绘制点之后连续M个二维的垂深-水平位移对序列中的点为基础,采用三次样条算法,形成回归方程;
第七步,对回归方程进行二阶求导,然后搜索拐点,从拐点处分道绘制测井曲线;
第八步,结束。
所述第一步是这样实现的:
以造斜点A和钻头位置C连线AC,过A点作垂线,其中B点是该垂线上的一点,然后做过连线ABC的平面,该平面为纵向地层剖面,标记为ABCD;作过AB线段的且垂直于ABCD剖面的剖面FGHI,FGHI为垂向地层剖面。纵向地层剖面ABCD与垂向地层剖面FGHI相交于线段AB。
所述第二步是这样实现的:
沿井眼轨迹,计算轨迹每一点相对于A点的沿AB方向的垂向位移zi(i=1、2、3、……);计算轨迹每一点相对于A点的沿AD方向的水平位移xi(i=1、2、3、……)。
所述第三步是这样实现的:
构建(z1,val1)、(z2,val2)、(z3,val3)……,形成二维的垂深-测井响应量数据对序列:(zi,vali)(i=1,2,3,……),其中,zi为垂直深度,vali为对应于垂深的测井物理量幅度值。
所述第四步是这样实现的:
构建(z1,x1)、(z2,x2)、(z3,x3)……,形成沿井眼轨迹构建二维的垂深-水平位移对序列:(zi,xi)(i=1,2,3,……),其中,x为水平位移量。
所述第五步是这样实现的:
在二维剖面中,以(zi,vali)为样本点,以垂深zi为垂向坐标、以value1为横向左刻度、以value2为横向右刻度、以vali数据为横向坐标,作为一个曲线绘制坐标绘制道,以(z1,val1)值为坐标数据,绘制第一点,并标记该位置为N0。
所述第六步是这样实现的:
取M为10,沿井眼轨迹,以当前绘制点之后连续10个点(zi,xi)(i=1、2、3、……、10)为基础,采用三次样条算法,形成回归方程:Z=aX3+bX2+cX+d。
所述第七步是这样实现的:
对回归方程Z=aX3+bX2+cX+d按x为自变量进行二阶求导,判断是否存在拐点;
如果存在拐点,则记录下该拐点N,采用第五步同样的刻度信息,自动创建一个新的曲线道,并以该拐点N为起始曲线绘制点,并重新标记起始点N0,此时,N0=N,然后返回第六步;
如果不存在拐点,则置N=M,自第五步绘制的数据点N0开始到拐点N为止,以垂深-测井响应量数据对序列中的下一个点为样本点,在第五步形成的绘制道中绘制数据的坐标点,并逐点连接相邻的坐标点,然后判断二维的垂深-测井响应量数据对序列是否全部完成,如果完成,则转入第八步,如果未完成,则返回第六步。
与现有技术相比,本发明的有益效果是:运用本发明的随钻测井曲线绘制方法,特别是通过计算机软件实现自动随钻测井曲线绘制,有效地解决了对于测井轨迹是在上下地层中来回反复穿行情况下,则投影后的曲线是重复的,人们无法分辨出地层的特征,垂向投影的方法失去了作用的问题。这样在地质导向应用过程,计算机软件可以快速直观的为地质导向工程师提供相关曲线展示,便于他们正确的做出决策。
附图说明
图1井眼轨迹投影剖面示意图
图2垂向曲线示意图
图3分道显示曲线示意图。
具体实施方式
下面结合附图对本发明作进一步详细描述:
本发明提出了一种基于二阶导数自动搜寻拐点自动分道显示的方法。方法可以在模型建立完成以后,由计算软件自动搜索拐点,自动分道显示,实现在随钻条件下的测井曲线的绘制。
本发明是通过采用下述技术方案实现的:
第一步,选择纵向地层投影剖面
以造斜点A和钻头位置C连线AC,过A点作垂线,其中B点是该垂线上的一点,然后做过连线ABC的平面,这里称该平面为纵向地层剖面,标记为ABCD;作过AB线段的且垂直于ABCD剖面的剖面FGHI,这里称FGHI为垂向地层剖面。纵向地层剖面ABCD与垂向地层剖面相交线段于AB。
图1为本发明定义的概念与标记符号。
第二步,将井眼轨迹数据投影到平面ABCD上,沿井眼轨迹,计算轨迹每一点相对于A点的沿AB方向的垂向位移zi(i=1、2、3、……);计算轨迹每一点相对于A点的沿AD方向的水平位移xi(i=1、2、3、……);
第三步,沿井眼轨迹构建二维的垂深-测井响应量数据对序列:(z1,value1)、(z2,value2)、(z3,value3)……。
第四步,沿井眼轨迹构建二维的垂深-水平位移对数据序列:(z1,x1)、(z2,x2)、(z3,x3)……。
第五步,在二维剖面中,以(zi,vali)为样本点,以垂深zi为垂向坐标;以value1为横向左刻度,以value2为横向右刻度,以vali数据为横向坐标,作为一个曲线绘制坐标绘制道。以(z1,val1)值为坐标数据,绘制第一点,并标记该位置为N0。(如图2所示)。
第六步,沿井眼轨迹,以当前绘制点之后连续10点(zi,xi)(i=1、2、3、……、10)序列点为基础,采用三次样条算法,形成回归方程:Z=aX3+bX2+cX+d;
第七步,对Z-X进行二阶求导,查找是否存在拐点,如果存在,则记录下该拐点位置N,采用第五步同样的刻度信息,自动创建一个新的曲线道,并以该拐点为起始曲线绘制点,绘制垂深-测井响应量数据对曲线。重复上述第五步后的过程,直到所有曲线绘制完成。其参考如图3所示。重复上述第五步后的过程,直到所有曲线绘制完成;如果不存在拐点,则置N=10。自第五步绘制的数据点N0开始到N数据位置为止,以(z2,val2)为样本点,在第五步形成的道中绘制数据的坐标点,并逐点连接相邻的坐标点。以前面的N对应点为开始点,并重新标记N0,(此时,N0=N),继续重复第六步后操纵作,直至(zi,vali)全部数据序列完成;
第八步,结束。
本发明涉及油气开发与勘探领域,具体用于随钻测井曲线的绘制,实现更加清晰的垂向的地层特征显示。此项发明的成果应用到随钻地质导向,将非常有效的提高地质导向效率,具有非常好的市场需求。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。

Claims (8)

1.一种随钻测井曲线的绘制方法,其特征在于:所述方法包括:
第一步,建立纵向地层剖面;
第二步,将井眼轨迹数据投影到纵向地层剖面上;
第三步,沿井眼轨迹构建二维的垂深-测井响应量数据对序列;
第四步,沿井眼轨迹构建二维的垂深-水平位移对序列;
第五步,建立曲线绘制坐标绘制道,并绘制二维的垂深-测井响应量数据对序列中的第一点N0;
第六步,沿井眼轨迹,以当前绘制点之后连续M个二维的垂深-水平位移对序列中的点为基础,采用三次样条算法,形成回归方程;
第七步,对回归方程进行二阶求导,然后搜索拐点,从拐点处分道绘制测井曲线;
第八步,结束。
2.根据权利要求1所述的随钻测井曲线的绘制方法,其特征在于:所述第一步是这样实现的:
以造斜点A和钻头位置C连线AC,过A点作垂线,其中B点是该垂线上的一点,然后做过连线ABC的平面,该平面为纵向地层剖面,标记为ABCD;作过AB线段的且垂直于ABCD剖面的剖面FGHI,FGHI为垂向地层剖面。纵向地层剖面ABCD与垂向地层剖面FGHI相交于线段AB。
3.根据权利要求2所述的随钻测井曲线的绘制方法,其特征在于:所述第二步是这样实现的:
沿井眼轨迹,计算轨迹每一点相对于A点的沿AB方向的垂向位移zi(i=1、2、3、……);计算轨迹每一点相对于A点的沿AD方向的水平位移xi(i=1、2、3、……)。
4.根据权利要求3所述的随钻测井曲线的绘制方法,其特征在于:所述第三步是这样实现的:
构建(z1,val1)、(z2,val2)、(z3,val3)……,形成二维的垂深-测井响应量数据对序列:(zi,vali)(i=1,2,3,……),其中,zi为垂直深度,vali为对应于垂深的测井物理量幅度值。
5.根据权利要求4所述的随钻测井曲线的绘制方法,其特征在于:所述第四步是这样实现的:
构建(z1,x1)、(z2,x2)、(z3,x3)……,形成沿井眼轨迹构建二维的垂深-水平位移对序列:(zi,xi)(i=1,2,3,……),其中,x为水平位移量。
6.根据权利要求5所述的随钻测井曲线的绘制方法,其特征在于:所述第五步是这样实现的:
在二维剖面中,以(zi,vali)为样本点,以垂深zi为垂向坐标、以value1为横向左刻度、以value2为横向右刻度、以vali数据为横向坐标,作为一个曲线绘制坐标绘制道,以(z1,val1)值为坐标数据,绘制第一点,并标记该位置为N0。
7.根据权利要求6所述的随钻测井曲线的绘制方法,其特征在于:所述第六步是这样实现的:
取M为10,沿井眼轨迹,以当前绘制点之后连续10个点(zi,xi)(i=1、2、3、……、10)为基础,采用三次样条算法,形成回归方程:Z=aX3+bX2+cX+d。
8.根据权利要求7所述的随钻测井曲线的绘制方法,其特征在于:所述第七步是这样实现的:
对回归方程Z=aX3+bX2+cX+d按x为自变量进行二阶求导,判断是否存在拐点;
如果存在拐点,则记录下该拐点N,采用第五步同样的刻度信息,自动创建一个新的曲线道,并以该拐点N为起始曲线绘制点,并重新标记起始点N0,此时,N0=N,然后返回第六步;
如果不存在拐点,则置N=M,自第五步绘制的数据点N0开始到拐点N为止,以垂深-测井响应量数据对序列中的下一个点为样本点,在第五步形成的绘制道中绘制数据的坐标点,并逐点连接相邻的坐标点,然后判断二维的垂深-测井响应量数据对序列是否全部完成,如果完成,则转入第八步,如果未完成,则返回第六步。
CN201710356097.3A 2017-05-19 2017-05-19 一种随钻测井曲线的绘制方法 Active CN108961352B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710356097.3A CN108961352B (zh) 2017-05-19 2017-05-19 一种随钻测井曲线的绘制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710356097.3A CN108961352B (zh) 2017-05-19 2017-05-19 一种随钻测井曲线的绘制方法

Publications (2)

Publication Number Publication Date
CN108961352A true CN108961352A (zh) 2018-12-07
CN108961352B CN108961352B (zh) 2021-07-20

Family

ID=64462080

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710356097.3A Active CN108961352B (zh) 2017-05-19 2017-05-19 一种随钻测井曲线的绘制方法

Country Status (1)

Country Link
CN (1) CN108961352B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112145156A (zh) * 2020-07-16 2020-12-29 中国石油大学(华东) 一种井眼轨迹自适应测斜计算方法
CN112346129A (zh) * 2020-10-29 2021-02-09 中国石油天然气集团有限公司 测井曲线提取及合成地震记录制作方法及装置
CN117523035A (zh) * 2024-01-05 2024-02-06 中国石油集团东方地球物理勘探有限责任公司 水平井剖面中曲线变密填充方法、装置、设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2446278C1 (ru) * 2010-11-02 2012-03-27 Открытое акционерное общество "Всероссийский нефтегазовый научно-исследовательский институт имени академика А.П. Крылова" (ОАО "ВНИИнефть") Способ разработки месторождений с высоковязкими нефтями и битумами системой горизонтально-наклонных скважин
CN102662671A (zh) * 2012-04-06 2012-09-12 中国石油大学(华东) 大斜度、水平井测井资料二维分解显示方法
CN105551083A (zh) * 2015-12-08 2016-05-04 中国石油天然气集团公司 一种刻画水平井井眼轨迹与地层关系的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2446278C1 (ru) * 2010-11-02 2012-03-27 Открытое акционерное общество "Всероссийский нефтегазовый научно-исследовательский институт имени академика А.П. Крылова" (ОАО "ВНИИнефть") Способ разработки месторождений с высоковязкими нефтями и битумами системой горизонтально-наклонных скважин
CN102662671A (zh) * 2012-04-06 2012-09-12 中国石油大学(华东) 大斜度、水平井测井资料二维分解显示方法
CN105551083A (zh) * 2015-12-08 2016-05-04 中国石油天然气集团公司 一种刻画水平井井眼轨迹与地层关系的方法

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112145156A (zh) * 2020-07-16 2020-12-29 中国石油大学(华东) 一种井眼轨迹自适应测斜计算方法
WO2022011700A1 (zh) * 2020-07-16 2022-01-20 中国石油大学(华东) 一种井眼轨迹自适应测斜计算方法
US11319796B2 (en) 2020-07-16 2022-05-03 China University Of Petroleum (East China) Method for self-adaptive survey calculation of wellbore trajectory
CN112346129A (zh) * 2020-10-29 2021-02-09 中国石油天然气集团有限公司 测井曲线提取及合成地震记录制作方法及装置
CN117523035A (zh) * 2024-01-05 2024-02-06 中国石油集团东方地球物理勘探有限责任公司 水平井剖面中曲线变密填充方法、装置、设备及存储介质
CN117523035B (zh) * 2024-01-05 2024-03-26 中国石油集团东方地球物理勘探有限责任公司 水平井剖面中曲线变密填充方法、装置、设备及存储介质

Also Published As

Publication number Publication date
CN108961352B (zh) 2021-07-20

Similar Documents

Publication Publication Date Title
CN102622526B (zh) 数字矿山中的巷道开掘搜索预测方法
CN102116871B (zh) 一种提高地震反演数据预测储层能力的方法
AU2012388240B2 (en) Methods and systems of incorporating pseudo-surface pick locations in seismic velocity models
CN109003330A (zh) 一种基于基岩边界约束的三维地层建模方法
US10598817B2 (en) Local layer geometry engine with work zone generated from buffer defined relative to a wellbore trajectory
CN108961352A (zh) 一种随钻测井曲线的绘制方法
CN108898670B (zh) 一种基于实体和剖面的三维地质建模方法
Huang et al. Modeling well interference and optimal well spacing in unconventional reservoirs using the fast marching method
US8593462B2 (en) Method and system of correlating a measured log to a predicted log
CN104251135B (zh) 大斜度井空间归位方法
Ozkaya FRACOR-software toolbox for deterministic mapping of fracture corridors in oil fields on AutoCAD platform
Ozkaya Geometry of layer-bound fractures based on fracture density and aperture-depth plots from resistivity image logs of deviated wells
Pietraszsek-Mattner et al. Predicting the impact of natural fractures on reservoir performance
CN110469321B (zh) 一种确定地层破裂压力梯度的录井方法
US11727191B2 (en) Process for highlighting text with varied orientation
Garvin Bridger Lake field, Summit County, Utah
CN111562614B (zh) 基于穿断构造解释的油气运聚模拟方法
Stright et al. An Integrated Approach to Evaluation of Horizontal Well Prospects in the Niobrara Shale
CN108510586B (zh) 一种构建地质导向地层模型的方法
Al-Maskeen et al. Horizontal Well Correlation Using Real Time Data and Log Prediction in Geosteering Complex Reservoirs of Saudi Arabia
Ozgen et al. Fracture Characterization for Integrated Studies: A new approach and its applications
Kowalchuk et al. A New Approach to Geosteering in Emerging Unconventional Plays
CN114814941A (zh) 推移中间层位断棱位置的方法
Shou et al. Automatic Formation Method of Prospecting Line Profile Map in Complex Coal Geology Based on AutoCAD
Suleiman et al. Maximizing Reservoir Exposure in the Thin Limestone of Ahmadi Formation, Bahrain

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