CN112180436A - 一种新的声波测井曲线合成方法 - Google Patents

一种新的声波测井曲线合成方法 Download PDF

Info

Publication number
CN112180436A
CN112180436A CN202010927794.1A CN202010927794A CN112180436A CN 112180436 A CN112180436 A CN 112180436A CN 202010927794 A CN202010927794 A CN 202010927794A CN 112180436 A CN112180436 A CN 112180436A
Authority
CN
China
Prior art keywords
curve
acoustic
logging curve
sandstone
acoustic logging
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
CN202010927794.1A
Other languages
English (en)
Other versions
CN112180436B (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
Daqing Oilfield Co Ltd
Original Assignee
Petrochina Co Ltd
Daqing Oilfield 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, Daqing Oilfield Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010927794.1A priority Critical patent/CN112180436B/zh
Publication of CN112180436A publication Critical patent/CN112180436A/zh
Application granted granted Critical
Publication of CN112180436B publication Critical patent/CN112180436B/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
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/364Seismic filtering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种新的声波测井曲线合成方法,用于解决常规方法无法完全区分井间储层岩石类型的问题。具体步骤包括:利用研究区测井曲线和岩芯资料,从测井曲线中选择一条区分砂岩和泥岩较好的敏感测井曲线;基于敏感测井曲线和声波测井曲线的拟合公式,将敏感曲线转化成为中间曲线;对声波测井曲线进行低通滤波,并与中间曲线进行合并计算,生成一条初始合成声波测井曲线;根据砂岩和泥岩在初始合成声波测井曲线的值域重叠范围,分别将砂岩和泥岩对应的值域范围进行减小和增大,生成最终的合成声波测井曲线。所述合成的声波测井曲线能够完全区分砂岩和泥岩,利用最终合成的声波测井曲线通过反距离插值算法得到合成声波三维模型,可以提高井间砂岩预测精度。

Description

一种新的声波测井曲线合成方法
技术领域
本发明涉及油田储层描述领域,特别涉及一种新的声波测井曲线合成方法。
背景技术
目前国内绝大多数油田已经进入开发中晚期,油藏含水率高、产量递减快,准确预测优质储层发育位置是保证油田稳产的关键,而测井曲线是准确预测储层发育位置的必备基础数据。测井曲线是钻井之后,通过专业的仪器测量地下地层的岩石类型而得到的曲线。目前测井曲线有不同分类,包括声波时差、密度、自然电位、电阻率、伽马等。其中声波测井曲线非常重要,可以用来描述储层类型、孔隙结构及作为连接地震资料的桥梁,而原始声波测井曲线由于分辨能力有限,用来预测储层时具有误差,需要技术人员利用其它多种类型的测井曲线合成一条新的声波测井曲线。
目前常用的声波测井曲线合成方法存在的问题是:不同类型岩石在合成的声波测井曲线值域范围内存在局部重叠,合成曲线无法完全区分储层岩石类型,导致无法准确预测有利储层分布位置,影响了油田油气描述及开发。
发明内容
本发明在于克服背景技术中存在的现有常规声波合成曲线方法无法完全区分储层岩石类型导致无法准确预测有利储层分布位置的问题,而提供一种新的声波测井曲线合成方法。该新的声波测井曲线合成方法,能够合成具有较高分辨能力的声波测井曲线,合成的声波测井曲线能够完全区分储层岩石类型;并且,利用新合成的声波测井曲线通过数学插值得到合成声波三维模型,可以提高井间有利储层预测精度。
本发明解决其问题可通过如下技术方案来达到:一种新的声波测井曲线合成方法,包括以下步骤:
(1)利用研究区测井曲线和岩芯资料,得到岩石类型解释结果,所述岩石类型包括砂岩和泥岩;
(2)根据步骤(1)得到的岩石类型解释结果,从测井曲线中选择出区分岩石类型最准确的一条测井曲线作为敏感曲线;
(3)将步骤(2)选择的敏感曲线和测井曲线作相关性曲线拟合,得到拟合公式,并利用拟合公式将敏感曲线转化成与声波测井曲线类型、量纲一致的曲线,称之为中间曲线;
(4)对原始声波测井曲线进行低通滤波,将滤波后的声波曲线与中间曲线进行合并计算,生成一条初始合成声波测井曲线;
(5)针对初始合成声波测井曲线进行校正,生成最终的合成声波测井曲线。
利用步骤(5)生成的最终的合成声波测井曲线通过反距离插值算法得到合成声波三维模型,实现有利储层砂岩分布的预测。
所述步骤(5)初始合成声波测井曲线进行校正的过程为:
1)砂岩和泥岩在初始合成声波曲线值域范围内存在局部重叠,计算重叠区的值域范围h=max(砂岩)-min(泥岩);
2)根据重叠区值域范围,分别将砂岩和泥岩对应的声波曲线值进行减小和增大,生成最终的合成声波测井曲线。
所述步骤4)低通滤波的公式为:
Y(n)=bX(n) + (1-b)Y(n-1),
式中b为滤波系数,取值0.1;X(n)为对应点的声波值;Y(n-1)为上次声波滤波输出值;Y(n)本次声波滤波输出值。
本发明与上述背景技术相比较可具有如下有益效果:本发明提供了一种新的声波测井曲线合成方法 ,该方法得到一条新的声波测井曲线,能够区分砂岩和泥岩的特征,并且利用该曲线通过反距离插值得到三维模型,可以提高井间砂岩的预测精度,剖面结果与后验井上的砂岩剖面和泥岩剖面对应关系好,符合率达到90%以上,比基于原始声波插值计算的模型提高15%以上。为油田高效开发提供了精准的储层描述结果。
附图说明:
附图1是本发明新的声波测井曲线合成方法流程图;
附图2是本发明实施例研究区测井曲线和储层岩石类型解释结果剖面;
附图3是本发明实施例敏感曲线(GR)和原始声波曲线(AC)关系图及拟合公式;
附图4是本发明实施例原始声波曲线(低通滤波)、中间曲线和初始合成声波曲线;
附图5是本发明实施例初始合成声波曲线(左)与最终合成声波曲线(右)直方图;
附图6是本发明实施例初始合成声波曲线(左)与最终合成声波曲线(右)单井剖面图;
附图7是本发明实施例合成声波测井曲线插值剖面图。
具体实施方式
以下就本发明步骤及所能产生的效果,配合附图、实施例详细说明如下:
实施例1
本发明专利针对大庆长垣油田北二西区块萨尔图油层,开展一种新的声波测井曲线合成方法研究,如图1所示,具体步骤包括:
(1)在北二西区块萨尔图油层,通过测井和钻井取芯的手段,得到测井曲线和钻井岩芯资料,测井曲线包括AC(声波曲线)、DEN、LLD、LLS、GR(伽马曲线);
(2)利用上述得到的测井曲线和钻井岩芯资料得到北二西区块所用井的砂岩剖面分布和泥岩剖面分布(如图2所示,图2是其中某一口井的研究区测井曲线和储层岩石类型解释结果剖面,最右列为砂岩剖面分布和泥岩剖面分布);
(3)从上述测井曲线中选择出一条与砂岩剖面和泥岩剖面对应关系最好的、最敏感的测井曲线,从图2中可以看到GR测井曲线与砂岩剖面和泥岩剖面对应关系最好、最敏感,GR测井曲线值分布在50~169之间,其中GR测井曲线高值区(134~169)对应的是泥岩剖面,GR曲线低值区(50~162)对应的是砂岩剖面;
(4)将GR测井曲线和声波测井曲线做散点交汇图(横坐标为GR测井曲线、纵坐标为声波测井曲线),可见GR曲线值分布在50~169之间,声波测井曲线分布在201~330之间,散点的声波曲线值随GR曲线值的增大而增大,大多数的散点呈条带状分布,可以在条带中心画出一个趋势线,利用趋势线上的两个散点坐标(80,230)和(150,300),计算得到线性拟合关系方程y=x+150,利用该关系将GR测井曲线转化成为与声波曲线(AC)数据类型、量纲一致的曲线,称为中间曲线,如图3所示;
(5)对声波测井曲线(AC)进行低通滤波,得到低通滤波声波曲线,低通滤波的算法公式为:Y(n)=bX(n) + (1-b)Y(n-1),
式中b为滤波系数,取值0.1;
X(n)为对应点的声波曲线(AC)值;
Y(n-1)为上次声波滤波输出值;
Y(n)本次声波滤波输出值,即低通滤波声波曲线;
(6)如图4所示,将低通滤波声波曲线与中间曲线进行合并计算,生成一条初始合成声波测井曲线,算法公式为:
G(n)=Y(n)-min{Y(1),Y(2)...Y(n)}+ Z(n),
Y(n)为对应点的低通滤波声波曲线值,
Z(n)为对应点的中间曲线值,
G(n)初始合成声波测井曲线;
(7)将初始合成声波测井曲线与砂岩剖面、泥岩剖面做直方分布图(见图5左图),可见,砂岩和泥岩在初始合成声波曲线值域范围内存在局部重叠,计算重叠区的值域范围h=max(砂岩)-min(泥岩),即h=312-284=28;
(8)根据重叠区值域范围,分别将砂岩和泥岩对应的初始合成声波曲线段进行减小和增大,对减小和增大的幅度进行了实际测试,得到了砂岩段初始合成声波曲线值减掉h/2=28/2=14,泥岩段初始合成声波曲线值加上h/2=28/2=14,生成最终的合成声波测井曲线,图5右图为最终合成声波曲线的直方图,砂岩段最终合成声波曲线的值域范围为<298,砂岩段最终合成声波曲线的值域范围为>298,可见砂岩和泥岩在最终合成声波曲线值域范围内不存在重叠了。图6为初始合成声波曲线与最终合成声波曲线对比剖面图,可见最终合成声波曲线与砂岩剖面和泥岩剖面对应关系更好,最终合成声波测井曲线能够完全将砂岩和泥岩区分开,进一步提高了声波曲线与砂岩剖面和泥岩剖面对应关系;
(9)利用上述方法将北二西区块萨尔图油层的所有井进行声波测井曲线合成,并以合成的声波曲线为基础,通过反距离加权算法进行空间插值,得到合成声波曲线预测模型。图7为三维模型中的一条二维剖面,可见剖面结果与后验井上的砂岩剖面和泥岩剖面对应关系好(参与井为参与反距离插值计算的井,后验井为没有参与计算的井,用来检验预测精度),符合率达到90%以上,比基于原始声波插值计算的模型提高15%以上。说明本发明方法提高了长垣油田北二西区块井间砂岩与泥岩的预测精度,能够为油田高效开发提供精准的储层描述结果。

Claims (5)

1.一种新的声波测井曲线合成方法,包括以下步骤:
(1)利用研究区测井曲线和岩芯资料,得到岩石类型解释结果,所述岩石类型包括砂岩和泥岩;
(2)根据步骤(1)得到的岩石类型解释结果,从测井曲线中选择出区分岩石类型最准确的一条测井曲线作为敏感曲线;
(3)将步骤(2)选择的敏感曲线和测井曲线作相关性曲线拟合,得到拟合公式,并利用拟合公式将敏感曲线转化成与声波测井曲线类型、量纲一致的曲线,称之为中间曲线;
(4)对原始声波测井曲线进行低通滤波,将滤波后的声波曲线与中间曲线进行合并计算,生成一条初始合成声波测井曲线;
(5)针对初始合成声波测井曲线进行校正,生成最终的合成声波测井曲线。
2.根据权利要求1所述的一种新的声波测井曲线合成方法,其特征在于:利用步骤(5)生成的最终的合成声波测井曲线通过反距离插值算法得到合成声波三维模型,实现有利储层砂岩分布的预测。
3.根据权利要求1或2所述的一种新的声波测井曲线合成方法,其特征在于:初始合成声波测井曲线进行校正的过程为:
(1)砂岩和泥岩在初始合成声波曲线值域范围内存在局部重叠,计算重叠区的值域范围h=max(砂岩)-min(泥岩);
(2)根据重叠区值域范围,分别将砂岩和泥岩对应的声波曲线值进行减小和增大,生成最终的合成声波测井曲线。
4.根据权利要求3所述的一种新的声波测井曲线合成方法,其特征在于:所述步骤(2)砂岩段对应的声波曲线值减掉h/2,泥岩段对应的声波曲线值加上h/2,生成最终的合成声波测井曲线。
5.根据权利要求1所述的一种新的声波测井曲线合成方法,其特征在于:所述步骤4)低通滤波的公式为:
Y(n)=bX(n) + (1-b)Y(n-1);
式中b为滤波系数,取值0.1;X(n)为对应点的声波值;Y(n-1)为上次声波滤波输出值;Y(n)本次声波滤波输出值。
CN202010927794.1A 2020-09-07 2020-09-07 一种声波测井曲线合成方法 Active CN112180436B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010927794.1A CN112180436B (zh) 2020-09-07 2020-09-07 一种声波测井曲线合成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010927794.1A CN112180436B (zh) 2020-09-07 2020-09-07 一种声波测井曲线合成方法

Publications (2)

Publication Number Publication Date
CN112180436A true CN112180436A (zh) 2021-01-05
CN112180436B CN112180436B (zh) 2023-11-28

Family

ID=73925859

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010927794.1A Active CN112180436B (zh) 2020-09-07 2020-09-07 一种声波测井曲线合成方法

Country Status (1)

Country Link
CN (1) CN112180436B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110261898A (zh) * 2019-04-26 2019-09-20 中国石油化工股份有限公司 基于地震岩石物理实验分析的测井与地震速度匹配方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110261898A (zh) * 2019-04-26 2019-09-20 中国石油化工股份有限公司 基于地震岩石物理实验分析的测井与地震速度匹配方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
付志方;张君;邢卫新;许浩;: "拟声波构建技术在砂泥岩薄互层储层预测中的应用", 石油物探, no. 04, pages 1 - 3 *
卢伟;: "拟声波重构技术在巴彦呼舒凹陷储层预测中的应用", 断块油气田, no. 03 *
李素贞;商建瓴;: "拟声波曲线重构反演在柳屯储层预测中的应用", 复杂油气藏 *
贺懿;刘怀山;毛传龙;姜绍辉;王建花;: "多曲线声波重构技术在储层预测中的应用研究", 石油地球物理勘探, no. 05 *
赵继龙;熊冉;陈戈;常少英;曹鹏;: "伽马拟声波分频重构反演在储层预测中的应用", 地球物理学进展, no. 04 *

Also Published As

Publication number Publication date
CN112180436B (zh) 2023-11-28

Similar Documents

Publication Publication Date Title
CN109102180B (zh) 致密砂岩储层双甜点评价的综合参数评价方法
CN108802812A (zh) 一种井震融合的地层岩性反演方法
CN109541685B (zh) 一种河道砂体识别方法
CN110441823B (zh) 基于多源数据融合的地层对比不确定性可视化方法
CN112505784B (zh) 一种煤矿井下盾构机掘进路线选择方法
CN109143399B (zh) 一种识别碳酸盐岩层序界面的方法
CN105863628A (zh) 一种油田开发阶段储层精细预测方法
CN116168224A (zh) 基于成像砾石含量的机器学习岩相自动识别方法
CN109100791B (zh) 基于纵横向空间约束的速度反演方法
CN112946743B (zh) 区分储层类型的方法
CN115793094B (zh) 一种曲线叠合重构识别复杂页岩层系岩性的方法及应用
CN112180436B (zh) 一种声波测井曲线合成方法
CN109188522B (zh) 速度场构建方法及装置
CN112906465B (zh) 一种基于地层因素的煤系地层声波曲线重构方法及系统
CN112147687B (zh) 一种储层含气性预测方法及预测系统
CN111308550B (zh) 一种页岩vti储层的各向异性参数多波联合反演方法
CN114200524A (zh) 一种基于人工智能深度学习的测井密度曲线校正方法
CN110687591B (zh) 基于先验数据的波形匹配确定煤层及围岩物性参数的方法
CN112987126A (zh) 一种岩石物理相的地球物理预测方法
CN110795513B (zh) 河流相源储异位型致密油气甜点区分布的预测方法
CN104375171B (zh) 一种高分辨率地震反演方法
CN113720745A (zh) 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法
CN116953797B (zh) 基于典型高频旋回等时对比的湖相薄层砂岩储层预测方法
Lavenu et al. Modeling Fractures in a Heterogeneous Carbonate Reservoir Onshore UAE: A Case Study
CN110320572B (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