CN106054254A - 一种测井曲线校正方法及装置 - Google Patents

一种测井曲线校正方法及装置 Download PDF

Info

Publication number
CN106054254A
CN106054254A CN201610621073.1A CN201610621073A CN106054254A CN 106054254 A CN106054254 A CN 106054254A CN 201610621073 A CN201610621073 A CN 201610621073A CN 106054254 A CN106054254 A CN 106054254A
Authority
CN
China
Prior art keywords
data
curve data
curve
correction
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
CN201610621073.1A
Other languages
English (en)
Other versions
CN106054254B (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 National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201610621073.1A priority Critical patent/CN106054254B/zh
Publication of CN106054254A publication Critical patent/CN106054254A/zh
Application granted granted Critical
Publication of CN106054254B publication Critical patent/CN106054254B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation
    • G01V2210/52Move-out correction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using well-logging

Landscapes

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

Abstract

本申请提供一种测井曲线校正方法及装置,其中,所述方法包括:获取预设测井的预设测井曲线数据;对所述预设测井曲线数据中待校正的当前曲线数据进行深度校正;建立交会图并从所述交会图中确定待校正的目标数据点,并对确定的目标数据点进行扩径校正;将扩径校正后的数据点替换校正之前的目标数据点,得到经过扩径校正后的当前曲线数据;根据预设标准井的标准测井曲线数据以及经过扩径校正后的预设测井曲线数据,建立交会图,并根据建立的交会图进行井间一致性校正,得到校正后的预设测井曲线数据。本申请实施方式提供的一种测井曲线校正方法及装置,能够提高测井校正的精度。

Description

一种测井曲线校正方法及装置
技术领域
本申请涉及地质勘探技术领域,特别涉及一种测井曲线校正方法及装置。
背景技术
野外测井曲线的测量过程中,往往会受到当时测井环境、仪器、地层等多方面因素的影响,导致所测的曲线数据不能真实反映地下地质情况,因此需要对使用的测井曲线数据进行校正。
现有技术中往往存在多种对测井曲线数据进行校正的方法,其中,一种方法是通过使用声波时差曲线数据与伽马曲线数据的交会来区分扩径段,利用伽马曲线数据受井径影响较小的特点,将其与受井径影响较大的曲线数据做交会,来识别需要被校正的井段,但是这种方法容易将指示油气的层段当成扩径层段来校正,从而造成错误的测井校正结果。另一种方法是通过数学方程拟合曲面来表示地质数据在工区内的变化趋势,首先找到充分逼近标准层的最佳趋势面方程,然后再计算趋势值与真实值之差,就可以得到校正量,但这种方法的计算量偏大,比较复杂,也往往导致最终的校正结果不准确。
应该注意,上面对技术背景的介绍只是为了方便对本申请的技术方案进行清楚、完整的说明,并方便本领域技术人员的理解而阐述的。不能仅仅因为这些方案在本申请的背景技术部分进行了阐述而认为上述技术方案为本领域技术人员所公知。
发明内容
本申请实施方式的目的在于提供一种测井曲线校正方法及装置,能够提高测井校正的精度。
为实现上述目的,本申请实施方式提供一种测井曲线校正方法,所述方法包括:获取预设测井的预设测井曲线数据,所述预设测井曲线数据包括声波时差曲线数据、密度曲线数据、伽马曲线数据、钻头直径曲线数据、井径曲线数据和电阻率曲线数据;从所述预设测井曲线数据中确定出控制曲线数据,并利用所述控制曲线数据对所述预设测井曲线数据中待校正的当前曲线数据进行深度校正;根据所述钻头直径曲线数据和所述井径曲线数据,确定井径差值曲线数据;基于经过深度校正的当前曲线数据、所述井径曲线数据和所述井径差值曲线数据建立交会图,所述交会图中分布经过深度校正的当前曲线数据中的数据点;从所述交会图中确定待校正的目标数据点,并对确定的目标数据点进行扩径校正;将扩径校正后的数据点替换所述待校正的目标数据点,得到经过扩径校正后的当前曲线数据;根据预设标准井的标准测井曲线数据以及经过扩径校正后的预设测井曲线数据建立交会图,并根据建立的交会图进行井间一致性校正,得到校正后的预设测井曲线数据。
为实现上述目的,本申请实施方式还提供一种测井曲线校正装置,所述装置包括:预设测井曲线数据获取单元,用于获取预设测井的预设测井曲线数据,所述预设测井曲线数据包括声波时差曲线数据、密度曲线数据、伽马曲线数据、钻头直径曲线数据、井径曲线数据和电阻率曲线数据;深度校正单元,用于从所述预设测井曲线数据中确定出控制曲线数据,并利用所述控制曲线数据对所述预设测井曲线数据中待校正的当前曲线数据进行深度校正;井径差值曲线数据确定单元,用于根据所述钻头直径曲线数据和所述井径曲线数据,确定井径差值曲线数据;交会图建立单元,用于基于经过深度校正的当前曲线数据、所述井径曲线数据和所述井径差值曲线数据建立交会图,所述交会图中分布经过深度校正的当前曲线数据中的数据点;扩径校正单元,用于从所述交会图中确定待校正的目标数据点,并对确定的目标数据点进行扩径校正;数据替换单元,用于将扩径校正后的数据点替换所述待校正的目标数据点,得到经过扩径校正后的当前曲线数据;一致性校正单元,用于根据预设标准井的标准测井曲线数据以及经过扩径校正后的预设测井曲线数据建立交会图,并根据建立的交会图进行井间一致性校正,得到校正后的预设测井曲线数据。
本申请实施方式利用控制曲线数据对预设测井曲线数据中待校正的当前曲线数据进行深度校正,从而使得待校正的当前曲线数据与所述控制曲线数据相似的曲线形态位于同一深度位置。进而采用交会图的方式来识别需要进行扩径校正的数据,从而可以利用多元线性拟合的方法对识别出的数据进行扩径校正,使得经过扩径校正后的当前曲线数据更加符合实际地层特征。最后根据预设标准井的标准测井曲线数据以及经过扩径校正后的预设测井曲线数据建立交会图,并根据建立的交会图进行井间一致性校正,使得校正后的预设测井曲线数据能够更加精确地反映地质的实际情况。本申请实施方式结合深度校正、扩径校正以及井间一致性校正,能够提高测井校正的精度。
参照后文的说明和附图,详细公开了本申请的特定实施方式,指明了本申请的原理可以被采用的方式。应该理解,本申请的实施方式在范围上并不因而受到限制。在所附权利要求的精神和条款的范围内,本申请的实施方式包括许多改变、修改和等同。
针对一种实施方式描述和/或示出的特征可以以相同或类似的方式在一个或更多个其它实施方式中使用,与其它实施方式中的特征相组合,或替代其它实施方式中的特征。
应该强调,术语“包括/包含”在本文使用时指特征、整件、步骤或组件的存在,但并不排除一个或更多个其它特征、整件、步骤或组件的存在或附加。
附图说明
所包括的附图用来提供对本申请实施方式的进一步的理解,其构成了说明书的一部分,用于例示本申请的实施方式,并与文字描述一起来阐释本申请的原理。显而易见地,下面描述中的附图仅仅是本申请的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本申请实施方式提供的一种测井曲线校正方法的流程图;
图2为本申请实施方式中深度校正前后的曲线数据对比示意图;
图3为本申请实施方式中声波时差曲线数据在扩径校正之前的交会图;
图4为本申请实施方式中密度曲线数据在扩径校正之前的交会图;
图5为本申请实施方式中扩径校正后的曲线数据对比示意图;
图6为本申请实施方式提供的一种测井曲线校正装置的功能模块图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施方式中的附图,对本申请实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式仅仅是本申请一部分实施方式,而不是全部的实施方式。基于本申请中的实施方式,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施方式,都应当属于本申请保护的范围。
图1为本申请实施方式提供的一种测井曲线校正方法的流程图。虽然下文描述流程包括以特定顺序出现的多个操作,但是应该清楚了解,这些过程可以包括更多或更少的操作,这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。如图1所示,所述方法包括:
步骤S1:获取预设测井的预设测井曲线数据,所述预设测井曲线数据包括声波时差曲线数据、密度曲线数据、伽马曲线数据、钻头直径曲线数据、井径曲线数据和电阻率曲线数据。
在本实施方式中,所述预设测井可以是位于某工区内的多个测井中的其中一个。在本实施方式中仅对所述预设测井的预设测井曲线数据进行校正的过程进行阐述,对于其它测井的测井曲线数据的校正过程可以此类推。
在本实施方式中,所述预设测井可以位于储层为低孔、低渗的白云岩储层,该储层发育裂缝和溶蚀孔洞,白云岩中含少量泥质,使得井径略有扩大,从而影响了所述预设测井曲线数据的形态和幅值,需要对所述预设测井曲线数据进行校正。
在本实施方式中,所述预设测井曲线数据可以包括声波时差曲线数据、密度曲线数据、伽马曲线数据、钻头直径曲线数据、井径曲线数据和电阻率曲线数据。在本实施方式中,可以将校正地震反演所必需的声波时差曲线数据和密度曲线数据作为待校正的曲线数据。
步骤S2:从所述预设测井曲线数据中确定出控制曲线数据,并利用所述控制曲线数据对所述预设测井曲线数据中待校正的当前曲线数据进行深度校正。
在本实施方式中,所述控制曲线数据可以是受环境影响较小的测井曲线数据。具体地,可以将所述预设测井曲线数据中测量半径范围最大的曲线数据确定为所述控制曲线数据。在实际应用场景中,可以将伽马曲线数据作为所述控制曲线数据,这样,便可以利用所述伽马曲线数据对所述预设测井曲线数据中待校正的当前曲线数据进行深度校正。具体地,所述待校正的当前曲线数据例如可以是声波时差曲线数据或密度曲线数据。
在本实施方式中,在进行深度校正之前,声波时差曲线数据和密度曲线数据中在预设深度处的曲线峰值并不对应一致。而在深度校正之后,声波时差曲线数据和密度曲线数据中在所述预设深度处的曲线峰值可以与伽马曲线数据中在所述预设深度处的曲线峰值对齐。请参阅图2,图2左侧的测井曲线为伽马曲线,中间的测井曲线为声波差值曲线,最右侧的测井曲线为密度曲线。其中,实线可以表示校正之前的测井曲线,虚线可以表示校正之后的测井曲线。那么从图2中可以看出,在进行深度校正之后,声波时差曲线和密度曲线中在5175米处对齐了曲线峰值。具体的深度校正的过程可以在现有技术中成熟的软件中进行,这里便不再赘述。
在本实施方式中,所述控制曲线数据也可以是其它受环境影响较小的曲线数据,例如自然电位曲线数据或者电阻率曲线数据等。
步骤S3:根据所述钻头直径曲线数据和所述井径曲线数据,确定井径差值曲线数据;
步骤S4:基于经过深度校正的当前曲线数据、所述井径曲线数据和所述井径差值曲线数据建立交会图,所述交会图中分布经过深度校正的当前曲线数据中的数据点。
在本实施方式中,在对当前曲线数据进行深度校正之后,为了使得当前曲线数据能够更加符合实际地层特征,可以对经过深度校正后的当前曲线数据继续进行扩径校正。具体地,在本实施方式中可以采用交会图的方式来识别需要进行扩径校正的数据,从而可以利用多元线性拟合的方法对识别出的数据进行扩径校正。
在本实施方式中,首先可以根据所述钻头直径曲线数据和所述井径曲线数据,确定井径差值曲线数据。具体地,在本实施方式中,可以按照下述公式确定所述井径差值曲线数据:
CALO=CAL-BIT
其中,CALO表示所述井径差值曲线数据,CAL表示所述井径曲线数据,BIT表示所述钻头直径曲线数据。
在确定出所述井径差值曲线数据之后,可以基于经过深度校正的当前曲线数据、所述井径曲线数据和所述井径差值曲线数据建立交会图。具体地,在本实施方式中,以经过深度校正的当前曲线数据为X轴,以所述井径曲线数据为Y轴,以所述井径差值曲线数据为色标,以所述经过深度校正的当前曲线数据的校正半径为色标范围,可以建立交会图。
请参阅图3,对于声波时差曲线数据而言,可以将声波时差曲线数据AC作为X轴,井径曲线数据CAL为Y轴,井径差值曲线数据CALO作为色标,声波时差曲线数据的校正半径0-4.8英寸作为色标范围,建立交会图。从图3中可以看出,所述交会图中分布着经过深度校正的声波时差曲线数据中的数据点。其中,不同形状的数据点对应着不同测井的经过深度校正的声波时差曲线数据。
请参阅图4,对于密度曲线数据而言,可以将密度曲线数据DEN作为X轴,井径曲线数据CAL为Y轴,井径差值曲线数据CALO作为色标,密度曲线数据的校正半径0-1英寸作为色标范围,建立交会图。从图4中可以看出,所述交会图中分布着经过深度校正的密度曲线数据中的数据点。其中,不同形状的数据点对应着不同测井的经过深度校正的密度曲线数据。
这样,所述交会图中分布的数据点与当前曲线数据值、井径差值和校正半径值相对应。在图3和图4中,黑色的数据点可以表示井径差值大于校正半径值的数据点,而灰色的数据点可以表示井径差值小于校正半径值的数据点。
步骤S5:从所述交会图中确定待校正的目标数据点,并对确定的目标数据点进行扩径校正。
在本实施方式中,对于图3和图4而言,已知白云石的速度为7010m/s,声波时差值约为43.5us/ft,密度值约为2.87g/cm3。但是同一岩性不同埋深的地层有着不同的速度和密度,加之孔隙度填充泥质、流体或其他矿物的影响,不同深度地层的速度和密度变化会比较大。由图3和图4中可以看出,声波时差值的大部分数据都在42-62us/ft范围之内,密度值的大部分数据都在2-3g/cm3范围之内,这些大部分数据所处的范围可以为预设范围。
在本实施方式中,所述交会图中分布的数据点可以分为四类:第一类为井径差值小于校正半径值,且当前曲线数据值在预设范围之内的数据点,这些数据点无需校正;第二类为井径差值小于校正半径值,但是当前曲线数据值在预设范围以外的数据点,这些数据点可能指示油气,因此也不要校正,若对于进行归位校正,会引起错误的测井解释结果;第三类为井径差值大于校正半径值,但是当前曲线数据值在预设范围以内的数据点,可以认为这些数据点没受扩径影响或已经经过初步校正,因此无需校正;第四类为井径差值大于校正半径值,并且当前曲线数据值在预设范围之外的数据点,这些数据点才是受扩径影响严重的数据点,要对其进行校正。
因此,在本实施方式中,可以将所述交会图中井径差值大于校正半径值,并且当前曲线数据值在预设范围之外的数据点确定为所述待校正的目标数据点。
在确定了所述目标数据点之后,便可以对确定的目标数据点进行扩径校正,以使得当前曲线数据能够更加符合实际地层特征。在本实施方式中,可以利用多元线性拟合的方法来拟合目标数据点所在深度段的曲线数据值。具体地,可以由受井径影响较小的伽马曲线和电阻率曲线,通过建立多元线性关系式来拟合声波时差曲线和密度曲线。在本实施方式中,当所述经过深度校正的当前曲线数据为声波时差曲线数据时,可以按照下述公式对确定的目标数据点进行扩径校正:
AC=41.447+0.417*GR-6.434*10-5*RT-2.715*10-6*RXO
其中,AC表示扩径校正后的声波时差数据点,GR为与所述声波时差数据点相对应的伽马曲线数据,RT为与所述声波时差数据点相对应的第一电阻率曲线数据,RXO为与所述声波时差数据点相对应的第二电阻率曲线数据。
当所述经过深度校正的当前曲线数据为密度曲线数据时,可以按照下述公式对确定的目标数据点进行扩径校正:
DEN=4.053-0.019*AC-0.032*GR+1.282*10-6*RT-6.082*10-6*RXO
其中,DEN表示扩径校正后的密度数据点,AC表示扩径校正后的声波时差数据点,GR为与所述声波时差数据点相对应的伽马曲线数据,RT为与所述声波时差数据点相对应的第一电阻率曲线数据,RXO为与所述声波时差数据点相对应的第二电阻率曲线数据。
需要说明的是,上述两个公式中的各个系数可以根据实际情况进行改变,本实施方式的两个公式中的系数可以应用于具体的某个实际应用场景中。本领域技术人员应当知晓,对于改变系数后的其它公式,也同样属于本申请的保护范围之内。
步骤S6:将扩径校正后的数据点替换所述待校正的目标数据点,得到经过扩径校正后的当前曲线数据。
在本实施方式中,在拟合得到校正后的数据点之后,便可以将扩径校正后的数据点替换校正之前的目标数据点,从而得到经过扩径校正后的当前曲线数据。请参阅图5,图5左侧为声波差值曲线,右侧为密度曲线。从图5中可以看出,在5210-5215m深度范围内,由于泥质含量过多引起了扩径现象,使得速度降低,声波时差升高,密度降低,经过校正后,曲线更加合理,且由于密度曲线更易受扩径的影响,其校正范围比声波时差曲线大。
步骤S7:根据预设标准井的标准测井曲线数据以及经过扩径校正后的预设测井曲线数据建立交会图,并根据建立的交会图进行井间一致性校正,得到校正后的预设测井曲线数据。
在本实施方式中,通过步骤S2至S6,可以对所述预设测井曲线数据中的各个曲线数据进行深度校正和扩径校正。为了使得最终的校正结果更加符合实际地质情况,在本实施方式中可以对经过扩径校正后的预设测井曲线数据进行井间一致性校正。具体地,在本实施方式中,同样可以通过交会图的方式来进行井间一致性校正。
在本实施方式中,可以将工区内测井曲线数据最齐全的测井作为标准井,所述标准井中的各个曲线数据也可以经过本申请中的深度校正和扩径校正的处理步骤,所述标准井中的测井曲线数据可以作为井间一致性校正的依据。具体地,在本实施方式中,可以从标准井的测井曲线数据中确定标准声波时差曲线数据和标准密度曲线数据,其中,可以将标准声波时差曲线数据转换为标准速度曲线数据(标准速度曲线数据为标准声波时差曲线数据的倒数),然后可以建立标准速度曲线数据和标准密度曲线数据的交会图。在所述交会图中可以包括标准速度曲线数据和标准密度曲线数据的概率分布直方图。然后可以确定预设分布概率值对应的标准速度值和标准密度值。在本实施方式中,所述预设分布概率值例如可以为80%,那么在80%的分布概率值下,标准速度值可以为6719m/s,标准密度值可以为2.83g/cm3,所述标准速度值和标准密度值可以作为第一标准值和第二标准值。
在本实施方式中,对于经过扩径校正后的预设测井曲线数据可以进行同样的处理。首先,可以确定所述预设测井对应的经过扩径校正后的第一曲线数据和第二曲线数据,并根据所述第一曲线数据和所述第二曲线数据建立二维交会图,所述二维交会图中包括第一曲线数据和第二曲线数据对应的概率分布直方图。所述第一曲线数据可以是通过对经过扩径校正的声波差值曲线数据进行求倒数而得到的速度曲线数据,所述第二曲线数据可以是经过扩径校正的密度曲线数据。
接着,可以从所述二维交会图中确定预设分布概率值对应的第一曲线数据值和第二曲线数据值。所述预设分布概率值可以为80%,那么在80%的分布概率值下,待校正的速度值可以为6620m/s,待校正的密度值可以为2.8g/cm3
这样,可以将所述第一曲线数据值与所述第一标准值进行对比并且将所述第二曲线数据值与所述第二标准值进行对比,分别得到所述第一曲线数据对应的第一校正量和所述第二曲线数据对应的第二校正量。具体地,速度校正量可以为第一校正量,该校正量可以为6620-6719=99m/s;密度校正量可以为第二校正量,该校正量可以为2.83-2.8=0.03g/cm3
最后,便可以利用所述第一校正量和所述第二校正量对所述第一曲线数据和所述第二曲线数据进行井间一致性校正,从而能够使得校正后的预设测井曲线数据能够更加符合实际地质情况。
请参阅图6,本申请还提供一种测井曲线校正装置,所述装置包括:
预设测井曲线数据获取单元100,用于获取预设测井的预设测井曲线数据,所述预设测井曲线数据包括声波时差曲线数据、密度曲线数据、伽马曲线数据、钻头直径曲线数据、井径曲线数据和电阻率曲线数据;
深度校正单元200,用于从所述预设测井曲线数据中确定出控制曲线数据,并利用所述控制曲线数据对所述预设测井曲线数据中待校正的当前曲线数据进行深度校正;
井径差值曲线数据确定单元300,用于根据所述钻头直径曲线数据和所述井径曲线数据,确定井径差值曲线数据;
交会图建立单元400,用于基于经过深度校正的当前曲线数据、所述井径曲线数据和所述井径差值曲线数据建立交会图,所述交会图中分布经过深度校正的当前曲线数据中的数据点;
扩径校正单元500,用于从所述交会图中确定待校正的目标数据点,并对确定的目标数据点进行扩径校正;
数据替换单元600,用于将扩径校正后的数据点替换所述待校正的目标数据点,得到经过扩径校正后的当前曲线数据;
一致性校正单元700,用于根据预设标准井的标准测井曲线数据以及经过扩径校正后的预设测井曲线数据建立交会图,并根据建立的交会图进行井间一致性校正,得到校正后的预设测井曲线数据。
在本申请的一个优选实施方式中,所述一致性校正单元700具体包括:
二维交会图建立模块,用于确定所述预设测井对应的经过扩径校正后的第一曲线数据和第二曲线数据,并根据所述第一曲线数据和所述第二曲线数据建立二维交会图,所述二维交会图中包括第一曲线数据和第二曲线数据对应的概率分布直方图;
数据值确定模块,用于从所述二维交会图中确定预设分布概率值对应的第一曲线数据值和第二曲线数据值;
标准值确定模块,用于从所述预设标准井的标准测井曲线数据中确定在所述预设分布概率值下,与所述第一曲线数据相对应的第一标准值以及与所述第二曲线数据相对应的第二标准值;
校正量获取模块,用于将所述第一曲线数据值与所述第一标准值进行对比并且将所述第二曲线数据值与所述第二标准值进行对比,分别得到所述第一曲线数据对应的第一校正量和所述第二曲线数据对应的第二校正量;
校正模块,用于利用所述第一校正量和所述第二校正量对所述第一曲线数据和所述第二曲线数据进行井间一致性校正。
在本申请的一个优选实施方式中,当所述经过深度校正的当前曲线数据为声波时差曲线数据时,所述扩径校正单元500按照下述公式对确定的目标数据点进行扩径校正:
AC=41.447+0.417*GR-6.434*10-5*RT-2.715*10-6*RXO
其中,AC表示扩径校正后的声波时差数据点,GR为与所述声波时差数据点相对应的伽马曲线数据,RT为与所述声波时差数据点相对应的第一电阻率曲线数据,RXO为与所述声波时差数据点相对应的第二电阻率曲线数据。
在本申请的一个优选实施方式中,当所述经过深度校正的当前曲线数据为密度曲线数据时,所述扩径校正单元500按照下述公式对确定的目标数据点进行扩径校正:
DEN=4.053-0.019*AC-0.032*GR+1.282*10-6*RT-6.082*10-6*RXO
其中,DEN表示扩径校正后的密度数据点。
需要说明的是,上述各个功能模块的具体实现方式以及计算公式均与步骤S1至S7中的描述一致,这里便不再赘述。
本申请实施方式利用控制曲线数据对预设测井曲线数据中待校正的当前曲线数据进行深度校正,从而使得待校正的当前曲线数据与所述控制曲线数据相似的曲线形态位于同一深度位置。进而采用交会图的方式对经过深度校正的当前曲线数据进行扩径校正,使得经过扩径校正后的当前曲线数据更加符合实际地层特征。最后根据预设标准井的标准测井曲线数据以及经过扩径校正后的预设测井曲线数据建立交会图,并根据建立的交会图进行井间一致性校正,使得校正后的预设测井曲线数据能够更加精确地反映地质的实际情况。本申请实施方式结合深度校正、扩径校正以及井间一致性校正,能够提高测井校正的精度。
在本说明书中,诸如第一和第二这样的形容词仅可以用于将一个元素或动作与另一元素或动作进行区分,而不必要求或暗示任何实际的这种关系或顺序。在环境允许的情况下,参照元素或部件或步骤(等)不应解释为局限于仅元素、部件、或步骤中的一个,而可以是元素、部件、或步骤中的一个或多个等。
上面对本申请的各种实施方式的描述以描述的目的提供给本领域技术人员。其不旨在是穷举的、或者不旨在将本发明限制于单个公开的实施方式。如上所述,本申请的各种替代和变化对于上述技术所属领域技术人员而言将是显而易见的。因此,虽然已经具体讨论了一些另选的实施方式,但是其它实施方式将是显而易见的,或者本领域技术人员相对容易得出。本申请旨在包括在此已经讨论过的本发明的所有替代、修改、和变化,以及落在上述申请的精神和范围内的其它实施方式。
本说明书中的各个实施方式均采用递进的方式描述,各个实施方式之间相同相似的部分互相参见即可,每个实施方式重点说明的都是与其他实施方式的不同之处。尤其,对于装置实施方式而言,由于其基本相似于方法实施方式,所以描述的比较简单,相关之处参见方法实施方式的部分说明即可。
虽然通过实施方式描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的精神。

Claims (10)

1.一种测井曲线校正方法,其特征在于,所述方法包括:
获取预设测井的预设测井曲线数据,所述预设测井曲线数据包括声波时差曲线数据、密度曲线数据、伽马曲线数据、钻头直径曲线数据、井径曲线数据和电阻率曲线数据;
从所述预设测井曲线数据中确定出控制曲线数据,并利用所述控制曲线数据对所述预设测井曲线数据中待校正的当前曲线数据进行深度校正;
根据所述钻头直径曲线数据和所述井径曲线数据,确定井径差值曲线数据;
基于经过深度校正的当前曲线数据、所述井径曲线数据和所述井径差值曲线数据建立交会图,所述交会图中分布经过深度校正的当前曲线数据中的数据点;
从所述交会图中确定待校正的目标数据点,并对确定的目标数据点进行扩径校正;
将扩径校正后的数据点替换所述待校正的目标数据点,得到经过扩径校正后的当前曲线数据;
根据预设标准井的标准测井曲线数据以及经过扩径校正后的预设测井曲线数据建立交会图,并根据建立的交会图进行井间一致性校正,得到校正后的预设测井曲线数据。
2.根据权利要求1所述的方法,其特征在于,从所述预设测井曲线数据中确定控制曲线数据具体包括:
将所述预设测井曲线数据中测量半径范围最大的曲线数据确定为所述控制曲线数据。
3.根据权利要求1所述的方法,其特征在于,按照下述公式确定所述井径差值曲线数据:
CALO=CAL-BIT
其中,CALO表示所述井径差值曲线数据,CAL表示所述井径曲线数据,BIT表示所述钻头直径曲线数据。
4.根据权利要求1所述的方法,其特征在于,基于经过深度校正的当前曲线数据、所述井径曲线数据和所述井径差值曲线数据建立交会图具体包括:
以经过深度校正的当前曲线数据为X轴,以所述井径曲线数据为Y轴,以所述井径差值曲线数据为色标,以所述经过深度校正的当前曲线数据的校正半径为色标范围建立交会图;
相应地,所述交会图中分布的数据点与当前曲线数据值、井径差值和校正半径值相对应。
5.根据权利要求4所述的方法,其特征在于,从所述交会图中确定待校正的目标数据点具体包括:
将所述交会图中井径差值大于校正半径值,并且当前曲线数据值在预设范围之外的数据点确定为所述待校正的目标数据点。
6.根据权利要求1所述的方法,其特征在于,当所述经过深度校正的当前曲线数据为声波时差曲线数据时,按照下述公式对确定的目标数据点进行扩径校正:
AC=41.447+0.417*GR-6.434*10-5*RT-2.715*10-6*RXO
其中,AC表示扩径校正后的声波时差数据点,GR为与所述声波时差数据点相对应的伽马曲线数据,RT为与所述声波时差数据点相对应的第一电阻率曲线数据,RXO为与所述声波时差数据点相对应的第二电阻率曲线数据。
7.根据权利要求6所述的方法,其特征在于,当所述经过深度校正的当前曲线数据为密度曲线数据时,按照下述公式对确定的目标数据点进行扩径校正:
DEN=4.053-0.019*AC-0.032*GR+1.282*10-6*RT-6.082*10-6*RXO
其中,DEN表示扩径校正后的密度数据点。
8.根据权利要求1所述的方法,其特征在于,根据预设标准井的标准测井曲线数据,对经过扩径校正后的预设测井曲线数据进行井间一致性校正具体包括:
确定所述预设测井对应的经过扩径校正后的第一曲线数据和第二曲线数据,并根据所述第一曲线数据和所述第二曲线数据建立二维交会图,所述二维交会图中包括第一曲线数据和第二曲线数据对应的概率分布直方图;
从所述二维交会图中确定预设分布概率值对应的第一曲线数据值和第二曲线数据值;
从所述预设标准井的标准测井曲线数据中确定在所述预设分布概率值下,与所述第一曲线数据相对应的第一标准值以及与所述第二曲线数据相对应的第二标准值;
将所述第一曲线数据值与所述第一标准值进行对比并且将所述第二曲线数据值与所述第二标准值进行对比,分别得到所述第一曲线数据对应的第一校正量和所述第二曲线数据对应的第二校正量;
利用所述第一校正量和所述第二校正量对所述第一曲线数据和所述第二曲线数据进行井间一致性校正。
9.一种测井曲线校正装置,其特征在于,所述装置包括:
预设测井曲线数据获取单元,用于获取预设测井的预设测井曲线数据,所述预设测井曲线数据包括声波时差曲线数据、密度曲线数据、伽马曲线数据、钻头直径曲线数据、井径曲线数据和电阻率曲线数据;
深度校正单元,用于从所述预设测井曲线数据中确定控制曲线数据,并利用所述控制曲线数据对所述预设测井曲线数据中待校正的当前曲线数据进行深度校正;
井径差值曲线数据确定单元,用于根据所述钻头直径曲线数据和所述井径曲线数据,确定井径差值曲线数据;
交会图建立单元,用于基于经过深度校正的当前曲线数据、所述井径曲线数据和所述井径差值曲线数据建立交会图,所述交会图中分布经过深度校正的当前曲线数据中的数据点;
扩径校正单元,用于从所述交会图中确定待校正的目标数据点,并对确定的目标数据点进行扩径校正;
数据替换单元,用于将扩径校正后的数据点替换所述待校正的目标数据点,得到经过扩径校正后的当前曲线数据;
一致性校正单元,用于根据预设标准井的标准测井曲线数据以及经过扩径校正后的预设测井曲线数据建立交会图,并根据建立的交会图进行井间一致性校正,得到校正后的预设测井曲线数据。
10.根据权利要求9所述的装置,其特征在于,所述一致性校正单元具体包括:
二维交会图建立模块,用于确定所述预设测井对应的经过扩径校正后的第一曲线数据和第二曲线数据,并根据所述第一曲线数据和所述第二曲线数据建立二维交会图,所述二维交会图中包括第一曲线数据和第二曲线数据对应的概率分布直方图;
数据值确定模块,用于从所述二维交会图中确定预设分布概率值对应的第一曲线数据值和第二曲线数据值;
标准值确定模块,用于从所述预设标准井的标准测井曲线数据中确定在所述预设分布概率值下,与所述第一曲线数据相对应的第一标准值以及与所述第二曲线数据相对应的第二标准值;
校正量获取模块,用于将所述第一曲线数据值与所述第一标准值进行对比并且将所述第二曲线数据值与所述第二标准值进行对比,分别得到所述第一曲线数据对应的第一校正量和所述第二曲线数据对应的第二校正量;
校正模块,用于利用所述第一校正量和所述第二校正量对所述第一曲线数据和所述第二曲线数据进行井间一致性校正。
CN201610621073.1A 2016-08-01 2016-08-01 一种测井曲线校正方法及装置 Active CN106054254B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610621073.1A CN106054254B (zh) 2016-08-01 2016-08-01 一种测井曲线校正方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610621073.1A CN106054254B (zh) 2016-08-01 2016-08-01 一种测井曲线校正方法及装置

Publications (2)

Publication Number Publication Date
CN106054254A true CN106054254A (zh) 2016-10-26
CN106054254B CN106054254B (zh) 2018-09-04

Family

ID=57197019

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610621073.1A Active CN106054254B (zh) 2016-08-01 2016-08-01 一种测井曲线校正方法及装置

Country Status (1)

Country Link
CN (1) CN106054254B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107092031A (zh) * 2017-06-01 2017-08-25 中国石油天然气集团公司 一种声波测井数据校正方法及装置
CN107991705A (zh) * 2017-10-13 2018-05-04 中国石油天然气股份有限公司 基于二维统计特征的测井曲线校正方法和装置
CN109386284A (zh) * 2018-11-23 2019-02-26 中国石油化工股份有限公司 基于多曲线联合校正的储层物性分析方法
CN110005396A (zh) * 2018-12-11 2019-07-12 中国石油天然气集团有限公司 一种测井波列数据深度弹性校正方法
CN110244385A (zh) * 2018-03-07 2019-09-17 中国石油化工股份有限公司 一种大斜度井测井曲线的校正方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5666057A (en) * 1996-02-29 1997-09-09 Western Atlas International, Inc. Method of skin effect correction and data quality verification for a multi-frequency induction well logging instrument
CN1243959A (zh) * 1999-06-24 2000-02-09 佟广新 同位素吸水剖面测井方法
CN103775074A (zh) * 2014-01-24 2014-05-07 西安石油大学 一种煤层气储层密度测井扩径影响的校正方法
CN103837893A (zh) * 2014-03-10 2014-06-04 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于vsp资料的声波测井曲线校正方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5666057A (en) * 1996-02-29 1997-09-09 Western Atlas International, Inc. Method of skin effect correction and data quality verification for a multi-frequency induction well logging instrument
CN1243959A (zh) * 1999-06-24 2000-02-09 佟广新 同位素吸水剖面测井方法
CN103775074A (zh) * 2014-01-24 2014-05-07 西安石油大学 一种煤层气储层密度测井扩径影响的校正方法
CN103837893A (zh) * 2014-03-10 2014-06-04 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于vsp资料的声波测井曲线校正方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
姜勇: "多参数拟合测井曲线环境校正方法研究及应用", 《海洋石油》 *
杨斌 等: "地震标定中的测井资料精细处理方法及应用", 《物探化探计算技术》 *
王方雄 等: "烃源岩测井评价新进展", 《测井技术》 *
董维武: "煤层气储层测井评价方法研究", 《中国优秀硕士学位论文全文数据库·基础科学辑》 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107092031A (zh) * 2017-06-01 2017-08-25 中国石油天然气集团公司 一种声波测井数据校正方法及装置
CN107092031B (zh) * 2017-06-01 2019-03-15 中国石油天然气集团公司 一种声波测井数据校正方法及装置
CN107991705A (zh) * 2017-10-13 2018-05-04 中国石油天然气股份有限公司 基于二维统计特征的测井曲线校正方法和装置
CN107991705B (zh) * 2017-10-13 2019-09-06 中国石油天然气股份有限公司 基于二维统计特征的测井曲线校正方法和装置
CN110244385A (zh) * 2018-03-07 2019-09-17 中国石油化工股份有限公司 一种大斜度井测井曲线的校正方法
CN109386284A (zh) * 2018-11-23 2019-02-26 中国石油化工股份有限公司 基于多曲线联合校正的储层物性分析方法
CN109386284B (zh) * 2018-11-23 2021-09-28 中国石油化工股份有限公司 基于多曲线联合校正的储层物性分析方法
CN110005396A (zh) * 2018-12-11 2019-07-12 中国石油天然气集团有限公司 一种测井波列数据深度弹性校正方法
CN110005396B (zh) * 2018-12-11 2023-02-10 中国石油天然气集团有限公司 一种测井波列数据深度弹性校正方法

Also Published As

Publication number Publication date
CN106054254B (zh) 2018-09-04

Similar Documents

Publication Publication Date Title
CN106054254A (zh) 一种测井曲线校正方法及装置
CN104047597B (zh) 含油气泥页岩地层测井曲线标准化方法
CN105488248B (zh) 一种深部矿层和岩层的判定方法
CN103670388A (zh) 一种泥页岩有机碳含量的评价方法
CN106154322B (zh) 测井曲线校正方法和装置
CN105317431A (zh) 一种用于解释和评价水平井测井参数的方法
CN103573250A (zh) 一种计算水平井井眼到地层上下界面距离的方法
CN104632202B (zh) 确定干粘土三孔隙度测井参数值的方法及装置
CN107807410A (zh) 储层地球物理敏感参数优选方法及系统
CN105317429A (zh) 一种确定页岩地层干酪根含量的方法
CN105626056A (zh) 一种确定目标区三维地层孔隙压力的方法和装置
CN106154342B (zh) 一种确定洞穴充填物电阻率的方法
CN106501146B (zh) 一种致密油藏物性上限确定方法及装置
CN113189672B (zh) 一种基于多属性反演的隧道超前地质预报方法
CN112415614B (zh) 一种测井探测尺度三维精细地层模型构建方法
CN107748393B (zh) 一种基于数值模拟的地层倾角对电阻率影响的校正方法
CN104895561A (zh) 基于并行计算技术的随钻电阻率测井联合反演方法
CN112115592A (zh) 一种基于岩石物理实验的大斜度/水平井电阻率校正方法
CN105464650A (zh) 一种随钻测井解释方法
CN111257932A (zh) 一种更新地层地震速度的方法
CN116108772A (zh) 地下工程岩体质量修正方法、装置、设备及可读存储介质
Walker et al. Stochastic inversion for facies: A case study on the Schiehallion field
CN114060015B (zh) 一种致密砂岩含气性的评价方法及装置
CN105178950B (zh) 一种确定细粒沉积岩中总有机碳的方法
CN108734781B (zh) 一种地层模型构建方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant