CN113356834B - 井眼声波时差的校正方法、校正装置和电子设备 - Google Patents

井眼声波时差的校正方法、校正装置和电子设备 Download PDF

Info

Publication number
CN113356834B
CN113356834B CN202110797553.4A CN202110797553A CN113356834B CN 113356834 B CN113356834 B CN 113356834B CN 202110797553 A CN202110797553 A CN 202110797553A CN 113356834 B CN113356834 B CN 113356834B
Authority
CN
China
Prior art keywords
borehole
acoustic
time difference
invasion
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.)
Active
Application number
CN202110797553.4A
Other languages
English (en)
Other versions
CN113356834A (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
Original Assignee
Petrochina 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 filed Critical Petrochina Co Ltd
Priority to CN202110797553.4A priority Critical patent/CN113356834B/zh
Publication of CN113356834A publication Critical patent/CN113356834A/zh
Application granted granted Critical
Publication of CN113356834B publication Critical patent/CN113356834B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本申请实施例公开一种井眼声波时差的校正方法、校正装置和电子设备。井眼声波时差的校正方法包括:获取井眼内杂质的侵入半径和井眼声波时差的理论值;基于井眼内杂质的侵入半径和井眼声波时差的理论值,建立侵入校正模型;获取井眼声波时差的斜度校正值;将井眼声波时差的斜度校正值输入侵入校正模型,得到校正后的井眼声波时差。本申请实施例能够在校正地层界面的倾斜角度对井眼声波时产生的影响的基础上,进一步校正井眼内杂质的侵入对井眼声波时差产生的影响,确保了校正后的井眼声波时差的准确性,使得井眼声波时差能够真实地反应井眼周围的地质情况,并且还能够提高井眼声波时差的校正效率。

Description

井眼声波时差的校正方法、校正装置和电子设备
技术领域
本申请实施例涉及井眼声波时差校正的技术领域,尤其涉及一种井眼声波时差的校正方法、一种井眼声波时差的校正装置和一种电子设备。
背景技术
页岩气已经从一种潜力资源逐渐转变为一种正在大规模开发利用的清洁能源,成为常规油气资源的有效接替。不同地区页岩物性的纵向分布、各向异性变化及构造等各不相同,如何更精准地获取不同地区页岩的各项异性等参数变得更为重要。
相关技术中,通过声波测井的方式来获取声波时差,将声波时差作为后续井震标定、页岩物理建模和压裂优化设计等过程的基础数据。但是,声波时差会受到井眼周围环境的影响,例如井眼内泥浆的侵入、地层的倾斜角度等因素,都会导致声波时差出现偏离,从而影响了通过声波时差得到的井眼周围页岩的地质情况的准确性。
发明内容
为了解决上述技术问题中至少之一,本申请实施例提供及一种井眼声波时差的校正方法、一种井眼声波时差的校正装置和一种电子设备。
第一方面,本申请实施例提供了一种井眼声波时差的校正方法,包括:获取井眼内杂质的侵入半径和井眼声波时差的理论值;基于井眼内杂质的侵入半径和井眼声波时差的理论值,建立侵入校正模型;获取井眼声波时差的斜度校正值;将井眼声波时差的斜度校正值输入侵入校正模型,得到校正后的井眼声波时差。
在一种可行的实施方式中,基于井眼内杂质的侵入半径和井眼声波时差的理论值,建立侵入校正模型的步骤包括:基于井眼内杂质的侵入半径,获取侵入声波时差;基于侵入声波时差和井眼声波时差的理论值,建立侵入校正模型。
在一种可行的实施方式中,基于井眼内杂质的侵入半径,获取侵入声波时差的步骤包括:根据公式
Figure BDA0003163286940000021
获取侵入声波时差;其中,y为侵入声波时差,x为井眼内杂质的侵入半径,y0为常数,A为常数,w为常数,xc为常数。
在一种可行的实施方式中,基于侵入声波时差和井眼声波时差的理论值,建立侵入校正模型的步骤包括:计算侵入声波时差和井眼声波时差的理论值的比值;根据比值和井眼内杂质的侵入半径,建立侵入校正模型。
在一种可行的实施方式中,获取井眼声波时差的斜度校正值的步骤包括:根据井眼的半径、声源的位置、声波接收器的位置、声波在杂质中的传播速度、声波在地层中的传播速度和地层的倾斜角度中至少之一,获取井眼声波时差的斜度校正值。
在一种可行的实施方式中,根据井眼的半径、声源的位置、声波接收器的位置、声波在杂质中的传播速度、声波在地层中的传播速度和地层的倾斜角度中至少之一,获取井眼声波时差的斜度校正值的步骤包括:根据公式
Figure BDA0003163286940000022
获取井眼声波时差的斜度校正值;其中,t为井眼声波时差的斜度校正值,r为井眼的半径,L为声源与声波接收器之间的直线距离,b为声波接收器与井眼和地层交点之间的直线距离,vf为声波在杂质中的传播速度,v为声波在地层中的传播速度,C为与地层的倾斜角度相关的变量,D为与地层的倾斜角度相关的变量,E为与地层的倾斜角度相关的变量。
在一种可行的实施方式中,声波接收器包括第一声波接收器和第二声波接收器,获取井眼声波时差的斜度校正值的步骤包括:通过第一声波接收器的位置,获取第一斜度校正值;通过第二声波接收器的位置,获取第二斜度校正值;基于第一斜度校正值和第二斜度校正值,获取井眼声波时差的斜度校正值。
在一种可行的实施方式中,基于第一斜度校正值和第二斜度校正值,获取井眼声波时差的斜度校正值的步骤包括:根据公式
Figure BDA0003163286940000031
获取井眼声波时差的斜度校正值;其中,T为井眼声波时差的斜度校正值,L1为第一声波接收器与声源的直线距离,L2为第二声波接收器与声源的直线距离,t1为第一斜度校正值,t2为第二斜度校正值。
第二方面,本申请实施例提供了一种井眼声波时差的校正装置,包括:获取装置,获取装置用于获取井眼内杂质的侵入半径和井眼声波时差的理论值;模型建立装置,模型建立装置基于井眼内杂质的侵入半径和井眼声波时差的理论值,建立侵入校正模型;斜度校正装置,斜度校正装置用于获取井眼声波时差的斜度校正值;输入装置,输入装置用于将井眼声波时差的斜度校正值输入侵入校正模型,得到校正后的井眼声波时差。
第三方面,本申请实施例提供了一种电子设备,包括处理器,存储器及存储在存储器上并可在处理器上运行的程序或指令,程序或指令被处理器执行时实现上述第一方面的井眼声波时差的校正方法的步骤。
本申请实施例有益效果如下:
通过侵入校正模型对井眼声波时差的斜度校正值进行校正,也即是在校正地层的倾斜角度对井眼声波时差产生的影响的基础上,进一步校正井眼内杂质的侵入对井眼声波时差产生的影响,确保了校正后的井眼声波时差的准确性,使得井眼声波时差能够真实地反应井眼周围的地质情况,为后续的地质建模、压裂优化等提供可靠数据。
此外,通过建立侵入校正模型的方式,来校正井眼声波时差的斜度校正值,进一步提高了对于井眼声波时差的校正效率。
本申请实施例的附加方面和优点将在下面的描述部分中变得明显,或通过本申请的实践了解到。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本申请的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1为本申请提供的一种实施例的井眼声波时差的校正方法的步骤流程图之一;
图2为本申请提供的一种实施例的井眼声波时差的校正方法的步骤流程图之二;
图3为本申请提供的一种实施例的侵入声波时差变化曲线示意图;
图4为本申请提供的一种实施例的井眼声波时差的校正方法的步骤流程图之三;
图5为本申请提供的一种实施例的侵入校正模型曲线示意图;
图6为本申请提供的一种实施例的井眼声波时差的测量示意图;
图7为本申请提供的一种实施例的井眼声波时差的校正方法的步骤流程图之四;
图8为本申请提供的一种实施例的井眼声波时差的校正装置结构示意框图;
图9为本申请提供的一种实施例的校正后的井眼声波时差的波形示意图。
其中,图6和图8中附图标记与部件名称之间的对应关系为:
T:声源,R1:第一声波接收器,R2:第二声波接收器,100:井眼声波时差的校正装置,110:获取装置,120:模型建立装置,130:斜度校正装置,140:输入装置。
具体实施方式
为了能够更清楚地理解本申请的上述目的、特征和优点,下面结合附图和具体实施方式对本申请进行进一步的详细描述。需要说明的是,在不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本申请,但是,本申请还可以采用其他不同于在此描述的其他方式来实施,因此,本申请的保护范围并不受下面公开的具体实施例的限制。
第一方面,如图1所示,本申请实施例提供了井眼声波时差的校正方法,包括:
步骤S101,获取井眼内杂质的侵入半径和井眼声波时差的理论值;
步骤S102,基于井眼内杂质的侵入半径和井眼声波时差的理论值,建立侵入校正模型;
步骤S103,获取井眼声波时差的斜度校正值;
步骤S104,将井眼声波时差的斜度校正值输入侵入校正模型,得到校正后的井眼声波时差。
可以理解地,井眼的声波时差为声源位置和声波接收器之间沿井壁传播的滑行纵波的时间差。获取井眼内杂质的侵入半径和井眼声波时差的理论值,可以理解地,井眼内的杂质可以为泥浆。井眼声波时差的理论值为不受环境因素的影响下,井眼的声波时差。
在一些示例中,井眼内杂质的侵入半径和井眼声波时差的理论值可以通过查阅资料、实地勘验、计算或者模拟仿真等方式获取。
基于井眼内杂质的侵入半径和井眼声波时差的理论值,建立侵入校正模型。可以理解地,侵入校正模型用于校正井眼内的杂质对于井眼声波时差的影响。
在一些示例中,根据侵入井眼内的杂质种类的不同,可以建立不同的侵入校正模型,提高侵入校正模型对于井眼声波时差校正的准确性和可靠性。
获取井眼声波时差的斜度校正值,可以理解地,井眼声波时差的斜度校正值为校正地层的倾斜角度影响之后的井眼声波时差。在一些示例中,井眼声波时差的斜度校正值可以通过查阅资料、计算、实地勘验或者模拟仿真等方式获取。
在一些示例中,井眼声波时差的斜度校正值随着地层的倾斜角度的变化而变化。
将井眼声波时差的斜度校正值输入侵入校正模型,可以理解地,通过侵入校正模型对井眼声波时差的斜度校正值进行校正,也即是在校正地层倾斜角度对井眼声波时差的影响的基础上,进一步校正井眼内杂质的侵入对井眼声波时差产生的影响,确保了校正后的井眼声波时差的准确性,使得井眼声波时差能够真实地反应井眼周围的地质情况,为后续的地质建模、压裂优化等提供可靠数据。此外,通过建立侵入校正模型的方式,来校正井眼声波时差的斜度校正值,进一步提高了对于井眼声波时差的校正效率。
在一些示例中,如图2所示,基于井眼内杂质的侵入半径和井眼声波时差的理论值,建立侵入校正模型的步骤包括:
步骤S201,基于井眼内杂质的侵入半径,获取侵入声波时差;
步骤S202,基于侵入声波时差和井眼声波时差的理论值,建立侵入校正模型。
可以理解地,侵入声波时差为在不同的侵入半径的影响下,井眼的声波时差。在一些示例中,可以通过查阅资料、实地勘验、计算或者模拟仿真等方式,来获取不同侵入半径下,井眼的侵入声波时差。
基于侵入声波时差和井眼声波时差的理论值,建立侵入校正模型,进一步提高了侵入校正模型的准确性,从而提高侵入校正模型对于井眼声波时差的斜度校正值的校正效果,确保了校正后的井眼声波时差的准确性,使得校正后的井眼声波时差能够准确地反应井眼的地质参数,提高井眼声波时差的校正方法的可靠性。
在一些示例中,基于井眼内杂质的侵入半径,获取侵入声波时差的步骤包括:
根据公式
Figure BDA0003163286940000071
获取侵入声波时差;
其中,y为侵入声波时差,x为井眼内杂质的侵入半径,y0为常数,A为常数,w为常数,xc为常数。
可以理解地,如图3所示,井眼的侵入声波时差与井眼内杂质的侵入半径存在对应的变化关系。在不同的侵入半径下,对井眼的侵入声波时差进行拟合,能够得到侵入声波时差随侵入半径的变化曲线,进而得到侵入声波时差的计算公式。
具体地,公式
Figure BDA0003163286940000072
中,y为侵入声波时差,x为井眼内杂质的侵入半径,y0、A、w和xc均为常数。
通过拟合公式,即可计算出侵入声波时差y,确保了得到了侵入声波时差的准确性,从而提高建立的侵入校正模型的准确性,提高井眼声波时差的校正方法的可靠性。
可以理解地,可以根据井眼的不同环境,设置y0、A、w和xc的不同取值。在一些示例中,y0的取值为260.827,xc的取值为0.01,w的取值为0.94,A的取值为1.815。
在一些示例中,如图4所示,基于侵入声波时差和井眼声波时差的理论值,建立侵入校正模型的步骤包括:
步骤S301,计算侵入声波时差和井眼声波时差的理论值的比值;
步骤S302,根据比值和井眼内杂质的侵入半径,建立侵入校正模型。
可以理解地,不同地质参数,对应不同的井眼声波时差的理论值,故而在不同的地质参数下,侵入声波时差与理论值的比值也不同。如图5所示,不同曲线代表不同地质参数下,比值与侵入半径之间的变化关系,提高了建立的侵入模型的适用性。
计算侵入声波时差和井眼声波时差的理论值的比值,根据侵入声波时差和井眼声波时差的比值与侵入半径,建立侵入校正模型,提高了侵入模型的准确性,从而提高校正后的声波时差的准确性,使得声波时差能够更为准确地反应地质情况,提高井眼声波时差的校正方法的可靠性。同时,还能够根据不同的地质参数,建立不同的侵入模型,提高了侵入模型的适用性,方法简单,易于实现,提高井眼声波的校正方法的效率。
在一些示例中,获取井眼声波时差的斜度校正值的步骤包括:根据井眼的半径、声源的位置、声波接收器的位置、声波在杂质中的传播速度、声波在地层中的传播速度和地层的倾斜角度中至少之一,获取井眼声波时差的斜度校正值。
可以理解地,在测量井眼的声波时差时,井眼的半径、声源的位置、声波接收器的位置、声波的传播速度和地层的倾斜角度都会对井眼声波时差造成影响。故而,根据井眼的半径、声源的位置、声波接收器的位置、声波在杂质中的传播速度、声波在地层中的传播速度和地层的倾斜角度中至少之一,来校正井眼声波时差,得到井眼声波时差的斜度校正值,减小了地层的倾斜角度对井眼声波时差造成的影响,提高了获取到的井眼声波时差的斜度校正值的准确性,从而提高井眼声波时差的校正方法的准确性。
在一些示例中,井眼周围存在多种不同地质,声波在不同地层中的传播速度并不相同。通过声波在不同地层中的传播速度来校正井眼声波时差,能够进一步提高井眼声波时差的斜度校正值的准确性,减小地层的倾斜角度对井眼声波时差造成的影响。
在一些示例中,根据井眼的半径、声源的位置、声波接收器的位置、声波在杂质中的传播速度、声波在地层中的传播速度和地层的倾斜角度中至少之一,获取井眼声波时差的斜度校正值的步骤包括:
根据公式
Figure BDA0003163286940000091
获取井眼声波时差的斜度校正值;
其中,t为井眼声波时差的斜度校正值,r为井眼的半径,L为声源与声波接收器之间的直线距离,b为声波接收器与井眼和地层交点之间的直线距离,vf为声波在杂质中的传播速度,v为声波在地层中的传播速度,C为与地层的倾斜角度相关的变量,D为与地层的倾斜角度相关的变量,E为与地层的倾斜角度相关的变量。
在计算井眼声波时差的斜度校正值时,r为井眼的半径,L为声源与声波接收器之间的直线距离,b为声波接收器与井眼和地层交点之间的直线距离。可以理解地,井眼与地层的交点可以为井眼中心轴线与地层之间的交点。vf为声波在杂质中的传播速度,在一些示例中,杂质可以为泥浆。v为声波在地层中的传播速度,在一些示例中,由于井眼周围存在多种不同的地质环境,声波在地层中的传播速度v的可以为多个不同数值。C、D和E均为与地层的倾斜角度相关的变量。
具体地,如图6所示,夹角θ为地层界面的倾斜角度,声源T发射出声波之后,声波在不同地层之间折射传播,最终被第一声波接收器R1接收。声波在折射传播过程中,形成夹角θ1、θ2、θ3和θ4。井眼声波时差的斜度校正值的计算公式为
Figure BDA0003163286940000101
其中,t为井眼声波时差的斜度校正值,r为井眼的半径,vf为声波在杂质中的传播速度,L1为声源T与第一声波接收器R1之间的直线距离,b1为第一声波接收器R1与井眼和地层交点之间的直线距离。v1和v2为声波在不同地层中的传播速度。
通过公式
Figure BDA0003163286940000102
计算获取井眼声波时差的斜度校正值,进一步提高了井眼声波时差的斜度校正值的准确性,从而提高井眼声波时差的校正方法的准确性和可靠性。
在一些示例中,声波接收器的数量可以为多个,多个声波接收器之间相互补偿校正,进一步提高井眼声波时差的斜度校正值的准确性。
在一些示例中,如图7所示,声波接收器包括第一声波接收器和第二声波接收器,获取井眼声波时差的斜度校正值的步骤包括:
步骤S401,通过第一声波接收器的位置,获取第一斜度校正值;
步骤S402,通过第二声波接收器的位置,获取第二斜度校正值;
步骤S403,基于第一斜度校正值和第二斜度校正值,获取井眼声波时差的斜度校正值。
可以理解地,如图6所示,声波接收器包括第一声波接收器R1和第二声波接收器R2。声源T发射出声波之后,声波在不同地层之间折射传播,最终被第一声波接收器R1和第二声波接收器R2接收。声波在折射传播过程中,形成夹角θ1、θ2、θ3、θ4、θ1’、θ2’、θ3’和θ4’。
具体地,可以根据公式
Figure BDA0003163286940000103
计算第一斜度校正值,其中t1为第一斜度校正值,r为井眼的半径,vf为声波在杂质中的传播速度,L1为声源T与第一声波接收器R1之间的直线距离,b1为声波接收器与井眼和地层交点之间的直线距离。v1和v2为声波在不同地层中的传播速度。
根据公式
Figure BDA0003163286940000111
计算第二斜度校正值,其中t2为第二斜度校正值,r为井眼的半径,vf为声波在杂质中的传播速度,L2为声源T与第二声波接收器R2之间的直线距离,b2为第二声波接收器R2与井眼和地层交点之间的直线距离。v1和v2为声波在不同地层中的传播速度。
通过第一声波接收器R1的位置获取第一斜度校正值,通过第二声波接收器R2的位置获取第二斜度校正值,第二斜度校正值与第一斜度校正值之间能够进行校正补偿,进而提高获取的井眼声波时差的斜度校正值的准确性,进一步提高井眼声波时差的校正方法的可靠性。
在一些示例中,基于第一斜度校正值和第二斜度校正值,获取井眼声波时差的斜度校正值的步骤包括:
根据公式
Figure BDA0003163286940000112
获取井眼声波时差的斜度校正值;
其中,T为井眼声波时差的斜度校正值,L1为第一声波接收器与声源的直线距离,L2为第二声波接收器与声源的直线距离,t1为第一斜度校正值,t2为第二斜度校正值。
通过公式
Figure BDA0003163286940000113
计算井眼声波时差的斜度校正值,使得第一斜度校正值和第二斜度校正值之间能够进行校正补偿,从而提高获取到的井眼声波时差的斜度校正值准确性,减小了地层倾斜角度对于声波时差的影响,提高井眼声波时差的校正方法的可靠性。
第二方面,如图8所示,基于上述方法的同一发明构思,本申请实施例提供一种井眼声波时差的校正装置100,包括获取装置110、模型建立装置120、斜度校正装置130和输入装置140。获取装置110用于获取井眼内杂质的侵入半径和井眼声波时差的理论值。模型建立装置120基于井眼内杂质的侵入半径和井眼声波时差的理论值,建立侵入校正模型。斜度校正装置130用于获取井眼声波时差的斜度校正值。输入装置140用于将井眼声波时差的斜度校正值输入侵入校正模型,得到校正后的井眼声波时差。本申请实施例中的井眼声波时差的校正装置能够实现上述第一方面中井眼声波时差的校正方法的步骤,因此具有上述第一方面的全部有益效果,在此不再赘述。
第三方面,基于上述方法的同一发明构思,本申请实施例提供一种电子设备,包括处理器,存储器及存储在存储器上并可在处理器上运行的程序或指令,程序或指令被处理器执行时实现上述第一方面的井眼声波时差的校正方法的步骤,因此具有上述第一方面的全部有益效果,在此不再赘述。
在本申请一个具体实施例中,如图6所示,对井眼声波时差进行校正时,首先校正地层的倾斜角度对于井眼声波时差的影响。
声源T发出的声波被不同的地层介质折射之后,被第一声波接收器R1和第二声波接收器R2接收。声波在折射过程中,形成夹角θ1、θ2、θ3、θ4、θ1’、θ2’、θ3’和θ4’。
根据第一声波接收器R1的设置位置,可以计算出第一斜度校正值
Figure BDA0003163286940000121
其中t1为第一斜度校正值,r为井眼的半径,vf为声波在杂质中的传播速度,L1为声源T与第一声波接收器R1之间的直线距离,b1为声波接收器与井眼和地层交点之间的直线距离。v1和v2为声波在不同地层中的传播速度。
根据第二声波接收器R2的设置位置,可以计算出第二斜度校正值
Figure BDA0003163286940000122
其中t2为第二斜度校正值,r为井眼的半径,vf为声波在杂质中的传播速度,L2为声源T与第二声波接收器R2之间的直线距离,b2为第二声波接收器R2与井眼和地层交点之间的直线距离。v1和v2为声波在不同地层中的传播速度。
根据公式
Figure BDA0003163286940000123
可以计算出井眼声波时差的斜度校正值,也即是对地层的倾斜角度进行了校正。
其中,T为井眼声波时差的斜度校正值,L1为第一声波接收器与声源的直线距离,L2为第二声波接收器与声源的直线距离,t1为第一斜度校正值,t2为第二斜度校正值。
根据公式
Figure BDA0003163286940000131
计算侵入声波时差,其中y为侵入声波时差,x为井眼内杂质的侵入半径,y0、A、w和xc均为常数。在一些示例中,y0的取值为260.827,xc的取值为0.01,w的取值为0.94,A的取值为1.815。
计算侵入声波时差和井眼声波时差的理论值的比值,根据比值和侵入半径,建立侵入校正模型。如图5所示,不同地质参数下,井眼声波时差的理论值不同,故而侵入声波时差与井眼声波时差理论值的比值也就不同。
将井眼声波时差的斜度校正值输入侵入校正模型,进一步校正侵入半径对于井眼声波时差的影响。具体地,可以将井眼声波时差的斜度校正值作为侵入声波时差带入到侵入校正模型中,查找对应侵入半径下,井眼声波时差的理论值,实现对与井眼声波时差的校正,得到校正后的井眼声波时差,减小了杂质侵入和地层倾斜角度对于井眼声波时差的影响。如图9所示,通过校正,能够较小井眼声波时差的偏离,提高井眼声波时差的准确性,为后续的地质建模、压裂优化等提供可靠数据。
在本发明中,术语“第一”、“第二”、“第三”仅用于描述的目的,而不能理解为指示或暗示相对重要性;术语“多个”则指两个或两个以上,除非另有明确的限定。术语“安装”、“相连”、“连接”、“固定”等术语均应做广义理解,例如,“连接”可以是固定连接,也可以是可拆卸连接,或一体地连接;“相连”可以是直接相连,也可以通过中间媒介间接相连。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
本发明的描述中,需要理解的是,术语“上”、“下”、“左”、“右”、“前”、“后”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或单元必须具有特定的方向、以特定的方位构造和操作,因此,不能理解为对本发明的限制。
在本申请的描述中,术语“一个实施例”、“一些实施例”、“具体实施例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或特点包含于本发明的至少一个实施例或示例中。在本申请中,对上述术语的示意性表述不一定指的是相同的实施例或实例。而且,描述的具体特征、结构、材料或特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种井眼声波时差的校正方法,其特征在于,包括:
获取井眼内杂质的侵入半径和井眼声波时差的理论值;
基于所述井眼内杂质的侵入半径和所述井眼声波时差的理论值,建立侵入校正模型;
获取井眼声波时差的斜度校正值;
将所述井眼声波时差的斜度校正值输入所述侵入校正模型,得到校正后的井眼声波时差;
所述基于所述井眼内杂质的侵入半径和所述井眼声波时差的理论值,建立所述侵入校正模型的步骤包括:
基于所述井眼内杂质的侵入半径,获取侵入声波时差;
基于所述侵入声波时差和所述井眼声波时差的理论值,建立所述侵入校正模型;
所述基于所述井眼内杂质的侵入半径,获取所述侵入声波时差的步骤包括:
根据公式获取所述侵入声波时差;
其中,y为所述侵入声波时差,x为所述井眼内杂质的侵入半径,y0为常数,A为常数,w为常数,xc为常数;
所述获取所述井眼声波时差的斜度校正值的步骤包括:
根据井眼的半径、声源的位置、声波接收器的位置、声波在杂质中的传播速度、声波在地层中的传播速度和地层的倾斜角度中至少之一,获取所述井眼声波时差的斜度校正值;
所述根据所述井眼的半径、所述声源的位置、所述声波接收器的位置、所述声波在杂质中的传播速度、所述声波在地层中的传播速度和所述地层的倾斜角度中至少之一,获取所述井眼声波时差的斜度校正值的步骤包括:
根据公式获取所述井眼声波时差的斜度校正值;
其中,t为所述井眼声波时差的斜度校正值,r为所述井眼的半径,L为所述声源与所述声波接收器之间的直线距离,b为所述声波接收器与井眼和地层交点之间的直线距离,vf为所述声波在杂质中的传播速度,v为所述声波在地层中的传播速度,C为与所述地层的倾斜角度相关的变量,D为与所述地层的倾斜角度相关的变量,E为与所述地层的倾斜角度相关的变量。
2.根据权利要求1所述的井眼声波时差的校正方法,其特征在于,所述基于所述侵入声波时差和所述井眼声波时差的理论值,建立所述侵入校正模型的步骤包括:
计算所述侵入声波时差和所述井眼声波时差的理论值的比值;
根据所述比值和所述井眼内杂质的侵入半径,建立所述侵入校正模型。
3.根据权利要求1所述的井眼声波时差的校正方法,其特征在于,所述声波接收器包括第一声波接收器和第二声波接收器,所述获取所述井眼声波时差的斜度校正值的步骤包括:
通过所述第一声波接收器的位置,获取第一斜度校正值;
通过所述第二声波接收器的位置,获取第二斜度校正值;
基于所述第一斜度校正值和所述第二斜度校正值,获取所述井眼声波时差的斜度校正值。
4.根据权利要求3所述的井眼声波时差的校正方法,其特征在于,所述基于所述第一斜度校正值和所述第二斜度校正值,获取所述井眼声波时差的斜度校正值的步骤包括:
根据公式获取所述井眼声波时差的斜度校正值;
其中,T为所述井眼声波时差的斜度校正值,L1为所述第一声波接收器与所述声源的直线距离,L2为所述第二声波接收器与所述声源的直线距离,t1为所述第一斜度校正值,t2为所述第二斜度校正值。
5.一种井眼声波时差的校正装置,其特征在于,用于执行如权利要求1至4中任一项所述的井眼声波时差的校正方法,包括:
获取装置,所述获取装置用于获取所述井眼内杂质的侵入半径和所述井眼声波时差的理论值;
模型建立装置,所述模型建立装置基于所述井眼内杂质的侵入半径和所述井眼声波时差的理论值,建立所述侵入校正模型;
斜度校正装置,所述斜度校正装置用于获取所述井眼声波时差的斜度校正值;
输入装置,所述输入装置用于将所述井眼声波时差的斜度校正值输入所述侵入校正模型,得到所述校正后的井眼声波时差;
所述基于所述井眼内杂质的侵入半径和所述井眼声波时差的理论值,建立所述侵入校正模型的步骤包括:
基于所述井眼内杂质的侵入半径,获取侵入声波时差;
基于所述侵入声波时差和所述井眼声波时差的理论值,建立所述侵入校正模型;
所述基于所述井眼内杂质的侵入半径,获取所述侵入声波时差的步骤包括:
根据公式获取所述侵入声波时差;
其中,y为所述侵入声波时差,x为所述井眼内杂质的侵入半径,y0为常数,A为常数,w为常数,xc为常数;
所述获取所述井眼声波时差的斜度校正值的步骤包括:
根据井眼的半径、声源的位置、声波接收器的位置、声波在杂质中的传播速度、声波在地层中的传播速度和地层的倾斜角度中至少之一,获取所述井眼声波时差的斜度校正值;
所述根据所述井眼的半径、所述声源的位置、所述声波接收器的位置、所述声波在杂质中的传播速度、所述声波在地层中的传播速度和所述地层的倾斜角度中至少之一,获取所述井眼声波时差的斜度校正值的步骤包括:
根据公式获取所述井眼声波时差的斜度校正值;
其中,t为所述井眼声波时差的斜度校正值,r为所述井眼的半径,L为所述声源与所述声波接收器之间的直线距离,b为所述声波接收器与井眼和地层交点之间的直线距离,vf为所述声波在杂质中的传播速度,v为所述声波在地层中的传播速度,C为与所述地层的倾斜角度相关的变量,D为与所述地层的倾斜角度相关的变量,E为与所述地层的倾斜角度相关的变量。
6.一种电子设备,其特征在于,包括处理器,存储器及存储在所述存储器上并可在所述处理器上运行的程序或指令,所述程序或指令被所述处理器执行时实现如权利要求1至4中任一项所述的井眼声波时差的校正方法的步骤。
CN202110797553.4A 2021-07-14 2021-07-14 井眼声波时差的校正方法、校正装置和电子设备 Active CN113356834B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110797553.4A CN113356834B (zh) 2021-07-14 2021-07-14 井眼声波时差的校正方法、校正装置和电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110797553.4A CN113356834B (zh) 2021-07-14 2021-07-14 井眼声波时差的校正方法、校正装置和电子设备

Publications (2)

Publication Number Publication Date
CN113356834A CN113356834A (zh) 2021-09-07
CN113356834B true CN113356834B (zh) 2023-04-18

Family

ID=77539506

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110797553.4A Active CN113356834B (zh) 2021-07-14 2021-07-14 井眼声波时差的校正方法、校正装置和电子设备

Country Status (1)

Country Link
CN (1) CN113356834B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101787884A (zh) * 2010-01-28 2010-07-28 中国石油集团川庆钻探工程有限公司 声波孔隙度和中子孔隙度差值储层流体类型判别方法
CN105089663A (zh) * 2015-07-31 2015-11-25 中国石油天然气集团公司 一种大斜度/水平井地层电阻率各向异性校正方法
CN106951660A (zh) * 2017-04-05 2017-07-14 中国石油天然气股份有限公司 一种海相碎屑岩水平井储层测井解释方法及装置
CN107678064A (zh) * 2016-08-02 2018-02-09 中石化石油工程技术服务有限公司 一种声波时差实时提取方法
CN110320562A (zh) * 2018-03-30 2019-10-11 中国石油化工股份有限公司 页岩气储层水平井中声波时差的校正方法
CN110826249A (zh) * 2019-11-18 2020-02-21 核工业北京地质研究院 一种铀矿综合测井曲线井眼环境校正方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10385678B2 (en) * 2014-03-21 2019-08-20 Conocophillips Company Method for analysing pore pressure in shale formations
US10012749B2 (en) * 2014-10-20 2018-07-03 Schlumberger Technology Corporation Fast model based inversion of acoustic impedance of annulus behind casing

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101787884A (zh) * 2010-01-28 2010-07-28 中国石油集团川庆钻探工程有限公司 声波孔隙度和中子孔隙度差值储层流体类型判别方法
CN105089663A (zh) * 2015-07-31 2015-11-25 中国石油天然气集团公司 一种大斜度/水平井地层电阻率各向异性校正方法
CN107678064A (zh) * 2016-08-02 2018-02-09 中石化石油工程技术服务有限公司 一种声波时差实时提取方法
CN106951660A (zh) * 2017-04-05 2017-07-14 中国石油天然气股份有限公司 一种海相碎屑岩水平井储层测井解释方法及装置
CN110320562A (zh) * 2018-03-30 2019-10-11 中国石油化工股份有限公司 页岩气储层水平井中声波时差的校正方法
CN110826249A (zh) * 2019-11-18 2020-02-21 核工业北京地质研究院 一种铀矿综合测井曲线井眼环境校正方法

Also Published As

Publication number Publication date
CN113356834A (zh) 2021-09-07

Similar Documents

Publication Publication Date Title
EP3433643B1 (en) Method and device for estimating sonic slowness in a subterranean formation
EP2093591B1 (en) Method for Three Dimensional Seismic Travel Time Tomography in Transversely Isotropic Media
US9529105B2 (en) Acoustic logging systems and methods employing multi-mode inversion for anisotropy and shear slowness
US8670288B2 (en) Velocity model for well time-depth conversion
US20070162249A1 (en) Traveltime calculation in three dimensional transversely isotropic (3D TTI) media by the fast marching method
WO2017210231A1 (en) Estimation of horizontal stresses and nonlinear constants in anisotropic formations such as interbedded carbonate layers in organic-shale reservoirs
RU2611096C2 (ru) Расчет параметра анизотропии
CN111257946B (zh) 一种地球物理指导钻井方法及更新地层地震速度的方法
CN103018778B (zh) 校正声波测井曲线进行储层预测的方法及设备
CN101285381B (zh) 一种泄漏模式波反演软地层横波速度的方法
US11237288B2 (en) Verifying measurements of elastic anisotropy parameters in an anisotropic wellbore environment
CN112034517B (zh) 一种地应力数据预测方法、装置、储存介质及电子设备
Tsang et al. Modeling shear rigidity of stratified bedrock in site response analysis
CA2731780C (en) Estimating formation stresses using radial profiles of three shear moduli
CN113356834B (zh) 井眼声波时差的校正方法、校正装置和电子设备
CN111257932B (zh) 一种更新地层地震速度的方法
US10989824B2 (en) Depth-dependent mud density determination and processing for horizontal shear slowness in vertical transverse isotropy environment using full-waveform sonic data
CN114740532B (zh) 一种水泥声阻抗、泥浆声阻抗的计算方法
US20140350857A1 (en) Method Of Mapping A Subterranean Formation Based Upon Wellbore Position And Seismic Data And Related System
CN112305595B (zh) 基于折射波分析地质体结构的方法及存储介质
CN113031067A (zh) 一种基于Rytov-WKBJ近似的叠前地震反演方法
CN109061737A (zh) 一种基于合成地震记录的储层预测方法及装置
US20240044239A1 (en) A method of and apparatus for optimizing a drilling strategy
US20230273332A1 (en) Generating low frequency models for seismic waveform inversion in formation regions with limited control wells
CN114718555A (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