CN110244385B - 一种大斜度井测井曲线的校正方法 - Google Patents

一种大斜度井测井曲线的校正方法 Download PDF

Info

Publication number
CN110244385B
CN110244385B CN201810186434.3A CN201810186434A CN110244385B CN 110244385 B CN110244385 B CN 110244385B CN 201810186434 A CN201810186434 A CN 201810186434A CN 110244385 B CN110244385 B CN 110244385B
Authority
CN
China
Prior art keywords
value
well
density
depth
highly
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
CN201810186434.3A
Other languages
English (en)
Other versions
CN110244385A (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 Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
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 Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201810186434.3A priority Critical patent/CN110244385B/zh
Publication of CN110244385A publication Critical patent/CN110244385A/zh
Application granted granted Critical
Publication of CN110244385B publication Critical patent/CN110244385B/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
    • G01V13/00Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00
    • 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)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Engineering & Computer Science (AREA)
  • Manufacturing & Machinery (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明公开了大斜度井测井曲线的校正方法,包括以下步骤:获取区域标准层内直井标准层的补偿中子值、密度值和深度值;建立直井标准层的补偿中子值和密度值随深度变化的关系式;获取大斜度井标准层补偿中子的理论值以及密度的理论值;获取大斜度井标准层补偿中子以及密度的实际测量值;获取大斜度井斜井段补偿中子以及密度的校正量;对大斜度井斜井段的补偿中子曲线和密度曲线进行校正。本发明提高了大斜度井中补偿中子曲线和密度曲线环境校正的精度,进而可以进一步指导油田的勘探开发工作。

Description

一种大斜度井测井曲线的校正方法
技术领域
本发明涉及油气勘探开发技术领域,尤其是涉及一种大斜度井测井曲线的校正方法。
背景技术
随着油田勘探开发技术的发展,大斜度井在复杂油气藏的应用日益广泛。实际应用中,由于受重力或井轨迹的影响,大斜度井斜井段中的中子、密度测井仪器易发生推靠井壁不紧密的现象,造成补偿中子和密度测井曲线的录取质量不高,和直井的测井响应特征存在较大差异,进而影响了大斜度井中储层孔隙度计算和测井储层评价的精度。而如何快速、准确地进行补偿中子、密度曲线校正,一直是大斜度井常规测井曲线环境校正的难题。
标准层法是目前进行补偿中子、密度等常规测井曲线校正的常用方法。该方法多用于直井常规测井曲线校正,前提是认为在工区范围内标准层的补偿中子、密度测井响应值是相同的,不会随着不同井标准层深度的差异而变化,而忽略了地层压实作用对标准层地层性质的影响,即不同深度的标准层应该具有不同的测井响应值,特别是在岩石成岩较弱的地层,如果不考虑地层深度的变化对测井响应值的影响,将极大影响测井曲线校正的精度,给储层孔隙度计算和测井储层评价带来很大困难。
发明内容
针对现有技术中所存在的上述技术问题,本发明提出了一种大斜度井测井曲线的校正方法,包括以下步骤:
S1、选取区域标准层,并获取所述区域标准层内直井标准层的补偿中子值、密度值和深度值;
S2、建立直井标准层的补偿中子值随深度变化的关系式,并建立直井标准层的密度值随深度变化的关系式;
S3、将大斜度井标准层内测量深度校正到垂直深度,获取大斜度井标准层补偿中子的理论值,并获取大斜度井标准层密度的理论值;
S4、获取大斜度井标准层补偿中子的实际测量值以及密度的实际测量值;
S5、根据大斜度井标准层补偿中子的理论值和实际测量值获取大斜度井斜井段补偿中子的校正量,并根据大斜度井标准层密度的理论值和实际测量值获取大斜度井斜井段密度的校正量;
S6、对大斜度井斜井段的补偿中子曲线和密度曲线进行校正。
优选地,步骤S1中,选取直井标准层的至少5个层段进行读值,补偿中子值和密度值选取每个层段的补偿中子值和密度值的平均值,深度值选取每个层段的中部深度值。
优选地,步骤S2中,直井标准层的补偿中子值随深度变化的关系式为:
NPHI=A*DEPTH+B
式中,NPHI为补偿中子值,DEPTH为深度值,A和B为常数。
优选地,步骤S2中,直井标准层的密度值随深度变化的关系式为:
RHOB=C*DEPTH+D
式中,RHOB为密度值,DEPTH为深度值,C和D为常数。
优选地,步骤S3包括:
S31、将大斜度井标准层内测量深度校正到垂直深度;
S32、根据获取的所述垂直深度和所述直井标准层的补偿中子值随深度变化的关系式,获取大斜度井标准层补偿中子的理论值;
S33、根据获取的所述垂直深度和所述直井标准层的密度值随深度变化的关系式,获取大斜度井标准层密度的理论值。
优选地,步骤S4中大斜度井标准层补偿中子的实际测量值的采样间隔、密度的实际测量值的采样间隔与步骤S3中测量深度值的采样间隔相同。
优选地,步骤S5中,大斜度井斜井段补偿中子的校正量为:
ΔNPHI==NPHI1-NPHI2
式中,ΔNPHI为大斜度井斜井段补偿中子的校正量,NPHI1为大斜度井斜井段补偿中子的理论值,NPHI2为大斜度井斜井段补偿中子的实际测量值。
优选地,步骤S5中,大斜度井斜井段密度的校正量为:
ΔRHOB=RHOB1-RHOB2
式中,ΔRHOB为大斜度井斜井段密度的校正量,RHOB1为大斜度井斜井段密度的理论值,RHOB2为大斜度井斜井段密度的实际测量值。
优选地,步骤S6中,将大斜度井斜井段补偿中子的实际测量值加上补偿中子的校正量ΔNPHI,以完成对大斜度井斜井段的补偿中子值的校正。
优选地,步骤S6中,将大斜度井斜井段密度的实际测量值加上密度的校正量ΔRHOB,以完成对大斜度井斜井段的密度值的校正。
与现有技术相比,本发明的优点在于,本发明提出的大斜度井测井曲线的校正方法,在标准层法的基础上,考虑了地层压实作用对不同深度标准层补偿中子、密度测井响应的影响,建立了补偿中子值、密度值随标准层深度变化的关系式,避免了以往标准层法忽略了压实作用对标准层性质变化影响的做法,能够准确、快速实现补偿中子和密度曲线的环境校正,提高了大斜度井中补偿中子和密度曲线环境校正的精度,进而提升了储层孔隙度计算和测井储层评价的可靠程度,从而指导油田的勘探开发工作。
附图说明
下面将结合附图来对本发明的优选实施例进行详细地描述。在图中:
图1是根据本发明实施例的一种大斜度井测井曲线的校正方法的流程图。
图2是根据本发明的一个实施例建立的大斜度井标准层的补偿中子值与深度变化关系图。
图3是根据本发明的一个实施例建立的大斜度井标准层的密度值与深度变化关系图。
图4是根据本发明的一个实施例对大斜度井补偿中子曲线和密度曲线的环境校正效果图。
在附图中,相同的部件使用相同的附图标记。附图并未按照实际的比例绘制。
具体实施方式
下面将结合附图对本发明做进一步说明。
图1显示了根据本发明实施例的提出了一种大斜度井测井曲线的校正方法,包括以下步骤:
S1、选取区域标准层,并获取所述区域标准层内直井标准层的补偿中子值、密度值和深度值;
S2、建立直井标准层的补偿中子值随深度变化的关系式,并建立直井标准层的密度值随深度变化的关系式;
S3、将大斜度井标准层内测量深度校正到垂直深度,获取大斜度井标准层补偿中子的理论值,并获取大斜度井标准层密度的理论值;
S4、获取大斜度井标准层补偿中子的实际测量值以及密度的实际测量值;
S5、根据大斜度井标准层补偿中子的理论值和实际测量值获取大斜度井斜井段补偿中子的校正量,并根据大斜度井标准层密度的理论值和实际测量值获取大斜度井斜井段密度的校正量;
S6、对大斜度井斜井段的补偿中子曲线和密度曲线进行校正。
本发明公开的大斜度井测井曲线的校正方法,在标准层法的基础上,考虑了地层压实作用对不同深度标准层补偿中子、密度测井响应的影响,建立了补偿中子值、密度值随标准层深度变化的关系式,避免了以往标准层法忽略了压实作用对标准层性质变化影响的做法,能够准确、快速实现补偿中子和密度曲线的环境校正,提高了大斜度井中补偿中子和密度曲线环境校正的精度,进而提升了储层孔隙度计算和测井储层评价的可靠程度,从而指导油田的勘探开发工作。
在一个实施例中,步骤S1中,需选取岩性纯、厚度大且分布稳定的页岩层作为区域标准层。为保证每口直井标准层读值的准确性,每口井的读值层段不少于5个。按照深度顺序,读取每个层段内补偿中子曲线和密度曲线的平均值,深度值取每个层段的中部深度值。
在一个实施例中,步骤S2中,将补偿中子值与深度值做线性回归,建立直井标准层的补偿中子值随深度变化的关系式为:
NPHI=A*DEPTH+B
式中,NPHI为补偿中子值,DEPTH为深度值,A和B为常数。
在一个实施例中,步骤S2中,将密度值与深度值做线性回归,建立直井标准层的密度值随深度变化的关系式为:
RHOB=C*DEPTH+D
式中,RHOB为密度值,DEPTH为深度值,C和D为常数。
在一个实施例中,步骤S3包括:
S31、将大斜度井标准层内测量深度校正到垂直深度;
S32、根据获取的所述垂直深度和所述直井标准层的补偿中子值随深度变化的关系式,获取大斜度井标准层补偿中子的理论值;
S33、根据获取的所述垂直深度和所述直井标准层的密度值随深度变化的关系式,获取大斜度井标准层密度的理论值。
在一个实施例中,步骤S4中大斜度井标准层补偿中子的实际测量值的采样间隔、密度的实际测量值的采样间隔与步骤S3中测量深度值的采样间隔相同。优选地,测量深度值的采样间隔为0.125m。
在一个实施例中,步骤S5中,大斜度井斜井段补偿中子的校正量为:
ΔNPHI==NPHI1-NPHI2
式中,ΔNPHI为大斜度井斜井段补偿中子的校正量,NPHI1为大斜度井斜井段补偿中子的理论值,NPHI2为大斜度井斜井段补偿中子的实际测量值。
在一个实施例中,步骤S5中,大斜度井斜井段密度的校正量为:
ΔRHOB=RHOB1-RHOB2
式中,ΔRHOB为大斜度井斜井段密度的校正量,RHOB1为大斜度井斜井段密度的理论值,RHOB2为大斜度井斜井段密度的实际测量值。
在一个实施例中,步骤S6中,将大斜度井斜井段补偿中子的实际测量值加上补偿中子的校正量ΔNPHI,以完成对大斜度井斜井段的补偿中子值的校正。
在一个实施例中,步骤S6中,将大斜度井斜井段密度的实际测量值加上密度的校正量ΔRHOB,以完成对大斜度井斜井段的密度值的校正。
以中国北部湾涠西断块区大斜度井涠4井为例,对本发明做进一步描述。选取延性纯、厚度大且分布稳定的流沙港组流二段泥页岩层作为区域标准层,读取研究区域内涠1、涠2、涠3、涠5、涠6等5口直井流二段的补偿中子值、密度值和深度值,每口井的读值层段均不少于5个。按照深度顺序,读取每个层段内补偿中子曲线和密度曲线的平均值,深度值取每个层段的中部深度值。
然后,如图2所示,将上述得到的补偿中子值与深度值做线性回归,建立直井标准层补偿中子值随深度变化的关系式:
NPHI=-0.0109*DEPTH+59.475
式中,NPHI为补偿中子值,%;DEPTH为深度值,m。
如图3所示,将上述得到的密度值与深度值做线性回归,建立直井标准层密度值随深度变化的关系式:
RHOB=0.0001*DEPTH+2.1372
式中:RHOB为密度值,g/cm3;DEPTH为深度值,m。
将涠4井标准层内任一深度点的测量深度校正到垂直深度。在一个优选的实施例中,测量深度值的采样间隔为0.125m。将得到的每个垂直深度代入上述直井标准层补偿中子值随深度变化的关系式,得到直井标准层每个垂直深度点的补偿中子的理论值,再根据每个垂直深度点的补偿中子的理论值,得到涠4井标准层补偿中子理论值的平均值NPHI1为29.5%。将得到的每个垂直深度代入上述直井标准层密度值随深度变化的关系式,得到直井标准层每个垂直深度点的密度值的理论值,再根据每个垂直深度点的密度值的理论值,得到涠4井标准层密度理论值的平均值RHOB1为2.49g/cm3
求取涠4井标准层补偿中子实际测量值的平均值NPHI2和密度实际测量值的平均值RHOB2,计算得到NPHI2为24.1%,RHOB2为2.60g/cm3
根据获取的涠4井标准层补偿中子理论值NPHI1和补偿中子实际测量值 NPHI2,得到涠4井标准层补偿中子的校正量:
ΔNPHI=NPHI1-NPHI2=29.5%-24.1%=5.4%
并且可以根据获取的涠4井标准层密度的理论值RHOB1和密度实际测量值 RHOB2,得到涠4井标准层密度的校正量:
ΔRHOB=RHOB1-RHOB2=2.49g/cm3-2.60g/cm3=-0.11g/cm3
如图4所示,对大斜度井涠4井斜井段的补偿中子和密度曲线进行环境校正。将上述步骤中得到的补偿中子的校正量ΔNPHI=5.4%加到涠4井斜井段补偿中子的实际测量值,即完成补偿中子曲线的环境校正。将上述步骤中得到的密度校正量ΔRHOB=-0.11g/cm3加到涠4井斜井段密度实际测量值,即完成密度曲线的环境校正。图4中,左侧实线为实际测量的补偿中子曲线,左侧虚线为矫正后的补偿中子曲线。右侧实线为实际测量的密度曲线,右侧虚线为矫正后的密度曲线。
以上所述仅为本发明的优选实施方式,但本发明保护范围并不局限于此,任何本领域的技术人员在本发明公开的技术范围内,可容易地进行改变或变化,而这种改变或变化都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求书的保护范围为准。

Claims (10)

1.一种大斜度井测井曲线的校正方法,其特征在于,包括以下步骤:
S1、选取区域标准层,并获取所述区域标准层内直井标准层的补偿中子值、密度值和深度值;
S2、建立直井标准层的补偿中子值随深度变化的关系式,并建立直井标准层的密度值随深度变化的关系式;
S3、将大斜度井标准层内测量深度校正到垂直深度,获取大斜度井标准层补偿中子的理论值,并获取大斜度井标准层密度的理论值;
S4、获取大斜度井标准层补偿中子的实际测量值以及密度的实际测量值;
S5、根据大斜度井标准层补偿中子的理论值和实际测量值获取大斜度井斜井段补偿中子的校正量,并根据大斜度井标准层密度的理论值和实际测量值获取大斜度井斜井段密度的校正量;
S6、对大斜度井斜井段的补偿中子曲线和密度曲线进行校正。
2.根据权利要求1所述的校正方法,其特征在于,步骤S1中,选取直井标准层的至少5个层段进行读值,补偿中子值和密度值选取每个层段的补偿中子值和密度值的平均值,深度值选取每个层段的中部深度值。
3.根据权利要求1所述的校正方法,其特征在于,步骤S2中,直井标准层的补偿中子值随深度变化的关系式为:
NPHI=A*DEPTH+B
式中,NPHI为补偿中子值,DEPTH为深度值,A和B为常数。
4.根据权利要求3所述的校正方法,其特征在于,步骤S2中,直井标准层的密度值随深度变化的关系式为:
RHOB=C*DEPTH+D
式中,RHOB为密度值,DEPTH为深度值,C和D为常数。
5.根据权利要求4所述的校正方法,其特征在于,步骤S3包括:
S31、将大斜度井标准层内测量深度校正到垂直深度;
S32、根据获取的所述垂直深度和所述直井标准层的补偿中子值随深度变化的关系式,获取大斜度井标准层补偿中子的理论值;
S33、根据获取的所述垂直深度和所述直井标准层的密度值随深度变化的关系式,获取大斜度井标准层密度的理论值。
6.根据权利要求1所述的校正方法,其特征在于,步骤S4中大斜度井标准层补偿中子的实际测量值的采样间隔、密度的实际测量值的采样间隔与步骤S3中测量深度值的采样间隔相同。
7.根据权利要求1所述的校正方法,其特征在于,步骤S5中,大斜度井斜井段补偿中子的校正量为:
ΔNPHI==NPHI1-NPHI2
式中,ΔNPHI为大斜度井斜井段补偿中子的校正量,NPHI1为大斜度井斜井段补偿中子的理论值,NPHI2为大斜度井斜井段补偿中子的实际测量值。
8.根据权利要求1所述的校正方法,其特征在于,步骤S5中,大斜度井斜井段密度的校正量为:
ΔRHOB=RHOB1-RHOB2
式中,ΔRHOB为大斜度井斜井段密度的校正量,RHOB1为大斜度井斜井段密度的理论值,RHOB2为大斜度井斜井段密度的实际测量值。
9.根据权利要求1所述的校正方法,其特征在于,步骤S6中,将大斜度井斜井段补偿中子的实际测量值加上补偿中子的校正量ΔNPHI,以完成对大斜度井斜井段的补偿中子值的校正。
10.根据权利要求1所述的校正方法,其特征在于,步骤S6中,将大斜度井斜井段密度的实际测量值加上密度的校正量ΔRHOB,以完成对大斜度井斜井段的密度值的校正。
CN201810186434.3A 2018-03-07 2018-03-07 一种大斜度井测井曲线的校正方法 Active CN110244385B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810186434.3A CN110244385B (zh) 2018-03-07 2018-03-07 一种大斜度井测井曲线的校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810186434.3A CN110244385B (zh) 2018-03-07 2018-03-07 一种大斜度井测井曲线的校正方法

Publications (2)

Publication Number Publication Date
CN110244385A CN110244385A (zh) 2019-09-17
CN110244385B true CN110244385B (zh) 2020-10-23

Family

ID=67882478

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810186434.3A Active CN110244385B (zh) 2018-03-07 2018-03-07 一种大斜度井测井曲线的校正方法

Country Status (1)

Country Link
CN (1) CN110244385B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112412435A (zh) * 2020-10-23 2021-02-26 中国石油集团工程咨询有限责任公司 一种用于螺旋井眼的密度曲线失真校正方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7191850B2 (en) * 2004-10-28 2007-03-20 Williams Danny T Formation dip geo-steering method
US7313479B2 (en) * 2005-01-31 2007-12-25 Baker Hughes Incorporated Method for real-time well-site interpretation of array resistivity log data in vertical and deviated wells
CN101981439A (zh) * 2008-01-04 2011-02-23 特克斯勒电子实验室公司 核子测量仪及核子测量仪的配置和校准方法
EP2241906B1 (en) * 2009-04-17 2015-04-01 Services Pétroliers Schlumberger Method of determining density of underground formations using neutron-gamma ray measurements
CN102662671B (zh) * 2012-04-06 2014-03-19 中国石油大学(华东) 大斜度、水平井测井资料二维分解显示方法
CN103984022B (zh) * 2014-04-16 2015-12-02 刘致水 大斜度井纵、横波测井速度各向异性同时校正方法
WO2015198137A2 (en) * 2014-06-25 2015-12-30 Cgg Services Sa Improving drilling operations using disparate well data types
CN104975853B (zh) * 2015-06-29 2018-01-05 中国石油天然气股份有限公司 一种获取地层孔隙度的方法及装置
CN106054254B (zh) * 2016-08-01 2018-09-04 中国石油天然气集团公司 一种测井曲线校正方法及装置

Also Published As

Publication number Publication date
CN110244385A (zh) 2019-09-17

Similar Documents

Publication Publication Date Title
CN105317375B (zh) 水平井引导入靶方法及装置
CN107907910B (zh) 一种不同岩性油藏横波测井确定方法
CN105114063B (zh) 一种斜井双侧向测井曲线校正方法
US8527204B2 (en) Volume of investigation based density image processing
CN105223617A (zh) 利用地震属性控制初始波阻抗的建模方法及系统
CN102428392A (zh) 使用中子伽玛射线测量结果确定地下地层的密度的方法
Mobach et al. In-situ reservoir compaction monitoring in the Groningen field
CN104516991B (zh) 一种伽马传感器全温度范围补偿方法
CN110244385B (zh) 一种大斜度井测井曲线的校正方法
CN105911592A (zh) 一种实钻约束的三维地震资料长波长静校正方法
CN110821403A (zh) 随钻井轨迹控制方法、装置及存储介质
CN112415614B (zh) 一种测井探测尺度三维精细地层模型构建方法
CN112115592A (zh) 一种基于岩石物理实验的大斜度/水平井电阻率校正方法
CN110488367B (zh) 一种基于阵列侧向测井资料的电阻率反演初值选取方法
CN109343099B (zh) 一种用gr和sp合并计算含放射性物质地层泥质含量的方法
CN111379551B (zh) 自然电位测井方法
CN107481298B (zh) 包含斜井的油藏剖面编绘方法
CN104879121A (zh) 总线结构随钻测量工具磁方位角补偿技术实现方法
CN105178950B (zh) 一种确定细粒沉积岩中总有机碳的方法
CN108734781B (zh) 一种地层模型构建方法
CN112131531B (zh) 基于最小二乘法的井温校正方法及装置
CN114109349B (zh) 一种致密砂岩储层孔隙度指数/饱和度指数的确定方法
CN113482596B (zh) 一种基于lwd随钻仪器实时三维伽马成像数据处理方法
CN109408837B (zh) 一种校正地质导向模型的方法
CN116816336A (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
GR01 Patent grant
GR01 Patent grant