CN111856571A - 测井频段速度向地震频段速度转换方法、装置及设备 - Google Patents
测井频段速度向地震频段速度转换方法、装置及设备 Download PDFInfo
- Publication number
- CN111856571A CN111856571A CN202010543114.6A CN202010543114A CN111856571A CN 111856571 A CN111856571 A CN 111856571A CN 202010543114 A CN202010543114 A CN 202010543114A CN 111856571 A CN111856571 A CN 111856571A
- Authority
- CN
- China
- Prior art keywords
- longitudinal wave
- velocity
- frequency
- core
- frequency band
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 239000012530 fluid Substances 0.000 claims abstract description 82
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 59
- 239000006185 dispersion Substances 0.000 claims abstract description 53
- 239000011435 rock Substances 0.000 claims abstract description 53
- 238000005070 sampling Methods 0.000 claims abstract description 39
- 238000005259 measurement Methods 0.000 claims description 51
- 229920006395 saturated elastomer Polymers 0.000 claims description 23
- 238000004590 computer program Methods 0.000 claims description 13
- 238000012937 correction Methods 0.000 claims description 13
- 238000006243 chemical reaction Methods 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 12
- 230000008569 process Effects 0.000 description 9
- 238000003860 storage Methods 0.000 description 9
- 230000006870 function Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 3
- 238000000691 measurement method Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010952 in-situ formation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XOOUIPVCVHRTMJ-UHFFFAOYSA-L zinc stearate Chemical compound [Zn+2].CCCCCCCCCCCCCCCCCC([O-])=O.CCCCCCCCCCCCCCCCCC([O-])=O XOOUIPVCVHRTMJ-UHFFFAOYSA-L 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/32—Transforming one recording into another or one representation into another
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6163—Electromagnetic
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (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
技术领域
本说明书涉及油气勘探技术领域,尤其是涉及一种测井频段速度向地震频段速度转换方法、装置及设备。
背景技术
含油气地下储层介质速度的大小与其储层岩石孔隙度渗透率参数、含流体类型等有很强的相关性,因此准确获取地下储层介质的速度信息对油气田储层位置的判断、储量预测、含流体类型预测具有重要指导意义。
在油气田的勘探开发中,对地下介质速度的测量与估算方法有很多种,最为常用的有声波时差测井方法、叠前地震反演方法和实验室岩芯测量方法。但是,声波时差测井方法与实验室岩芯测量方法的尺度都较小,无法反应整个区域地下介质的速度信息;在大尺度范围上,叠前地震反演方法是比较合适的速度获取方法。因此通常使用测井数据约束,验证地震反演得到的速度信息,这样能够使反演速度尽可能地准确。而当声波时差测井曲线缺失时,可以通过其余的数据质量可靠的测井曲线进行模拟和预测,实验室岩芯测量的速度常被用来验证这些预测得到的测井曲线。
在实现本申请的过程中,本申请的发明人发现:由于饱和流体多孔介质的频散性,会使波在地下介质传播的过程中产生频率依赖性。在上述三种方法测量中,叠前地震反演方法的测量频率在10-100Hz左右;声波时差测井方法的测量频率在10KHz左右;实验室岩芯测量方法的测量频率在兆赫兹级别;考虑到频散的影响,这三种不同频率测量方法测得的速度会有较大的差别,不能直接相互比对。
然而,在常用的叠前反演方法中,井震标定是连接测井与地震的桥梁,是叠前反演的重要步骤,传统做法将测井曲线(纵波速度、横波速度和密度)与地震子波褶积得到合成地震记录,然后与实际地震记录相匹配。这忽视了频散现象的存在,直接将两个不同频段测量的数据进行匹配,会因为人为地调整两者间的时深关系而改变了真实的时深关系,最终对反演结果产生影响。不仅如此,在叠前反演中,常用测井频段速度作为反演过程的约束条件及验证手段,但是同样未考虑速度频散的影响,这样会使地震反演的结果更加靠近测井频段的测量结果,丢失了地震频段测量中原本包含的一些信息。
综上所述,现有技术所获得的储层介质速度未考虑速度频散的影响,不够准确。如果后续将这种不准确储层介质速度的应用于井震匹配和反演约束等工作,会产生不利影响。
发明内容
本说明书实施方案的目的在于提供一种测井频段速度向地震频段速度转换方法、装置及设备,以获得更为准确的储层介质速度。
为达到上述目的,一方面,本说明书实施方案提供了一种测井频段速度向地震频段速度转换方法,包括:
获取不同深度取样的岩芯在实验室条件下的低、高频纵波速度,并从测井曲线中提取所述岩芯在对应取芯深度的中频纵波速度;
根据所述岩芯的实际地层压力、实际含流体性校正所述低、高频纵波速度,获得校正后的低、高频纵波速度;
根据所述中频纵波速度以及所述校正后的低、高频纵波速度,生成所述岩芯的全频段速度频散曲线;
根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
本说明书实施方案的测井频段速度向地震频段速度转换方法中,所述根据所述岩芯的实际地层压力校正所述低、高频纵波速度,包括:
根据所述岩芯对应的取样深度及地层压力梯度确定所述岩芯的实际地层压力;
确认所述低、高频纵波速度是否为,在围压为所述实际地层压力下测得的纵波速度;
其中,σ为围压,e为自然常数,为压力系数,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kdry-ini、Gdry-ini分别为干燥岩芯在零围压下的体积模量、剪切模量,Kdry-HP、Gdry-HP分别为干燥岩芯在上极限围压下的体积模量、剪切模量,VP为纵波速度,VS为横波速度;
将所述实际地层压力输入所述对应关系式,获得校正后的低、高频纵波速度。
本说明书实施方案的测井频段速度向地震频段速度转换方法中,所述根据所述岩芯的实际含流体性校正所述低、高频纵波速度,包括:
根据测井数据解释成果与录井信息确定所述岩芯的实际含流体性;
确认所述低、高频纵波速度是否为,在含流体性为所述实际含流体性下测得的纵波速度;
其中,Ksat(σ)、Gsat(σ)分别为饱和流体岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kf为饱和流体的体积模量,φ为饱和流体岩芯的孔隙度,Kma为饱和流体岩芯的岩石骨架体积模量,σ为围压,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,VP为纵波速度,VS为横波速度。
本说明书实施方案的测井频段速度向地震频段速度转换方法中,所述根据所述中频纵波速度以及所述校正后的低、高频纵波速度,生成所述岩芯的全频段速度频散曲线,包括:
拟合所述中频纵波速度以及所述校正后的低、高频纵波速度,获得所述岩芯的全频段速度频散曲线。
本说明书实施方案的测井频段速度向地震频段速度转换方法中,所述根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段,包括:
确定所述全频段速度频散曲线中,测井频段与地震主频间的纵波速度差平均值;
对所有岩芯取样深度范围内的测井曲线减去所述纵波速度差平均值,从而将所述岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
另一方面,本说明书实施方案还提供了一种测井频段速度向地震频段速度转换装置,包括:
获取模块,用于获取不同深度取样的岩芯在实验室条件下的低、高频纵波速度,并从测井曲线中提取所述岩芯在对应取芯深度的中频纵波速度;
校正模块,用于根据所述岩芯的实际地层压力、实际含流体性校正所述低、高频纵波速度,获得校正后的低、高频纵波速度;
生成模块,用于根据所述中频纵波速度以及所述校正后的低、高频纵波速度,生成所述岩芯的全频段速度频散曲线;
转换模块,用于根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
本说明书实施方案的测井频段速度向地震频段速度转换装置中,所述根据所述岩芯的实际地层压力校正所述低、高频纵波速度,包括:
根据所述岩芯对应的取样深度及地层压力梯度确定所述岩芯的实际地层压力;
确认所述低、高频纵波速度是否为,在围压为所述实际地层压力下测得的纵波速度;
其中,σ为围压,e为自然常数,为压力系数,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kdry-ini、Gdry-ini分别为干燥岩芯在零围压下的体积模量、剪切模量,Kdry-HP、Gdry-HP分别为干燥岩芯在上极限围压下的体积模量、剪切模量,VP为纵波速度,VS为横波速度;
将所述实际地层压力输入所述对应关系式,获得校正后的低、高频纵波速度。
本说明书实施方案的测井频段速度向地震频段速度转换装置中,所述根据所述岩芯的实际含流体性校正所述低、高频纵波速度,包括:
根据测井数据解释成果与录井信息确定所述岩芯的实际含流体性;
确认所述低、高频纵波速度是否为,在含流体性为所述实际含流体性下测得的纵波速度;
其中,Ksat(σ)、Gsat(σ)分别为饱和流体岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kf为饱和流体的体积模量,φ为饱和流体岩芯的孔隙度,Kma为饱和流体岩芯的岩石骨架体积模量,σ为围压,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,VP为纵波速度,VS为横波速度。
本说明书实施方案的测井频段速度向地震频段速度转换装置中,所述根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段,包括:
确定所述全频段速度频散曲线中,测井频段与地震主频间的纵波速度差平均值;
对所有岩芯取样深度范围内的测井曲线减去所述纵波速度差平均值,从而将所述岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
另一方面,本说明书实施方案还提供了一种电子设备,包括存储器、处理器、以及存储在所述存储器上的计算机程序,所述计算机程序被所述处理器运行时执行如下步骤:
获取不同深度取样的岩芯在实验室条件下的低、高频纵波速度,并从测井曲线中提取所述岩芯在对应取芯深度的中频纵波速度;
根据所述岩芯的实际地层压力、实际含流体性校正所述低、高频纵波速度,获得校正后的低、高频纵波速度;
根据所述中频纵波速度以及所述校正后的低、高频纵波速度,生成所述岩芯的全频段速度频散曲线;
根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
由以上本说明书实施方案提供的技术方案可见,在本说明书实施方案中,由于可以根据岩芯的实际地层压力、实际含流体性,将不同深度取样的岩芯在实验室条件下的低、高频纵波速度进行校正,从而使得校正后的低、高频纵波速度可以与从测井曲线中提取岩芯在对应取芯深度的中频纵波速度,均是在相同的地层压力和含流体性下得到的。在此基础上,就可以根据中频纵波速度以及校正后的低、高频纵波速度生成全频段速度频散曲线,进而可以根据全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段。从而实现在考虑到储层介质速度的频散现象基础上的速度校正,因而有利于获得更为准确的储层介质速度,进而为后续精确地完成井震匹配等工作提供了基础。
附图说明
为了更清楚地说明本说明书实施方案或现有技术中的技术方案,下面将对实施方案或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本说明书中记载的一些实施方案,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本说明书提供的实施方案中测井频段速度向地震频段速度转换方法的流程图;
图2a-图2d为本说明书提供的实施方案中不同岩芯在饱和水、不同压力下的频散曲线示意图;
图3a为本说明书提供的实施方案中饱和水条件下低频(30Hz)纵波速度随压力变化示意图;
图3b为本说明书提供的实施方案中饱和水条件下高频纵波速度随压力变化示意图;
图4为本说明书提供的实施方案中取芯位置(取芯深度)对应的测井曲线示意图;
图5a-图5d为本说明书提供的实施方案中不同岩芯的全频段频散曲线示意图;
图6为本说明书提供的实施方案中转换前的井震标定示意图;
图7为本说明书提供的实施方案中转换后的井震标定示意图;
图8为本说明书提供的实施方案中测井频段速度向地震频段速度转换装置的结构框图;
图9为本说明书提供的实施方案中电子设备的结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案,下面将结合本说明书实施方案中的附图,对本说明书实施方案中的技术方案进行清楚、完整地描述,显然,所描述的实施方案仅仅是本说明书一部分实施方案,而不是全部的实施方案。基于本说明书中的实施方案,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方案,都应当属于本说明书保护的范围。
参考图1所示,本说明书一些实施方案的测井频段速度向地震频段速度转换方法可以包括以下步骤:
S101、获取不同深度取样的岩芯在实验室条件下的低、高频纵波速度,并从测井曲线中提取所述岩芯在对应取芯深度的中频纵波速度。
在本说明书一些实施方案中,可以将不同深度取样的岩芯进行实验室岩芯应力-应变低频测量,从而可以得到不同流体特性(流体类型、流体饱和度)、不同压力条件下的低频岩石纵波速度。本说明书中的低频范围一般为3~3000Hz。此外,还可以通过超声脉冲透射设备(例如高温高压超声测定系统等),测量不同深度取样的岩芯在这些条件下的高频纵波速度。所取岩芯都来自于储层,所取岩芯的取样深度不能相差太远。本说明书中的高频一般不低于106Hz,例如,在一示例性实施方案中,高频可以为30MHz。
如图2a~图2d所示,在一示例性实施方案中,通过上述方式可以获取58号、59号、62号和63号岩芯在不同条件下的低频纵波速度。从图2a~图2d中可以看出,随着频率的增大,纵波速度会变大,从频率为10HZ到100HZ速度差值在50~100m/s,从而体现了纵波速度的频散现象。此外,通过超声脉冲透射设备还可以获取58号、59号、62号和63号岩芯的高频纵波速度。
在本说明书一些实施方案中,根据岩芯的取芯位置(即取芯深度),可以从中频(例如10KHz)声波时差测井曲线中提取岩芯在对应取芯深度的中频纵波速度。如图4所示,在一示例性实施方案中,一些岩芯(图4中部纵向分布的黑色圆点)的取芯深度大概为3150米~3200米之间,从该取芯深度所对应的测井曲线中即可以得到对应的中频纵波速度。为了提高测井曲线的精度,在提取中频纵波速度前,还可以对测井曲线原始测井曲线进行常规预处理及矫正(例如单井深度校正和多井一致性校正等)工作。
S102、根据所述岩芯的实际地层压力、实际含流体性校正所述低、高频纵波速度,获得校正后的低、高频纵波速度。
由于从测井曲线中提取出的中频纵波速度是通过实际地下的测井数据得到,即中频纵波速度是在储层原位地层压力(实际地层压力)条件下测量获取的。而在实验室条件下,岩芯的低、高频纵波速度的测量是将地下岩芯钻取至地面进行测量,所以理想条件下实验室的岩石物理测量其压力条件需要完全与中频纵波速度的测量吻合。然而,在实验室实际测量时,考虑到设备、样品等客观实验条件,一般难以完全满足实际地层压力条件。因此,需要通过实验室测量的纵波速度来估算获取岩芯对应深度处实际地层压力下的低、高频纵波速度(例如图3a和图3b所示)。
在本说明书一些实施方案中,所述根据所述岩芯的实际地层压力校正所述低、高频纵波速度,可以如下步骤:
11)、根据所述岩芯对应的取样深度及地层压力梯度确定所述岩芯的实际地层压力。其中,地层压力梯度可以根据实际情况确定。例如,在一示例性实施方案中,地层压力梯度例如可以为0.01Mpa/m。
12)、确认所述低、高频纵波速度是否为,在围压为所述实际地层压力下测得的纵波速度。
其中,σ为围压,e为自然常数,为压力系数,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kdry-ini、Gdry-ini分别为干燥岩芯在零围压下的体积模量、剪切模量,Kdry-HP、Gdry-HP分别为干燥岩芯在上极限围压下的体积模量、剪切模量,VP为纵波速度,VS为横波速度。
当然,如果所述低、高频纵波速度是在围压为所述实际地层压力下测得的纵波速度,则可以不用对其进行校正。但是,一般情况下,不需要进行校正的情况几乎不存在。
14)、将所述实际地层压力输入所述对应关系式,获得校正后的低、高频纵波速度。
不同流体饱和情况对纵波速度频散有较大的影响,因此需要在低频、中频、高频测量纵波速度时,保证储层岩石的含流体性条件完全一致(即(即具有相同的流体类型及流体饱和度)。从测井曲线中提取出的中频纵波速度是通过实际地下的测井数据得到,因此其对应的含流体性反映了储层原位地层的实际含流体性。而取样岩芯的实验室测量中受限于各种实验条件,一般难以使岩芯的含流体性完全符合储层原位地层的实际含流体性。因此,需要对低、高频纵波速度进行含流体性校正,从而使其含流体性可以与储层原位地层的实际含流体性保持一致。
在本说明书一些实施方案中,所述根据所述岩芯的实际含流体性校正所述低、高频纵波速度可以包括如下步骤:
21)、根据测井数据解释成果与录井信息确定所述岩芯的实际含流体性。本说明书中提及的含流体性可以包括流体类型和流体饱和度等。
22)、确认所述低、高频纵波速度是否为,在含流体性为所述实际含流体性下测得的纵波速度。
其中,Ksat(σ)、Gsat(σ)分别为饱和流体岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kf为饱和流体的体积模量,φ为饱和流体岩芯的孔隙度,Kma为饱和流体岩芯的岩石骨架体积模量,可以通过测井矿物解释或岩芯的实验室矿物分析得到,σ为围压,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,VP为纵波速度,VS为横波速度。
类似地,如果所述低、高频纵波速度是在含流体性为所述实际含流体性下测得的纵波速度,则可以不用对其进行校正。
本领域技术人员可以理解,上述根据岩芯的实际地层压力、实际含流体性校正低、高频纵波速度,实际上可以是在先进行其中一项校正的基础上,再进行另一项校正,本说明书对于具体校正顺序不作限定,具体可以根据需要选择。例如,可以先根据岩芯的实际地层压力校正低、高频纵波速度,再在此基础上根据岩芯的实际含流体性进行进一步校正。
S103、根据所述中频纵波速度以及所述校正后的低、高频纵波速度,生成所述岩芯的全频段速度频散曲线。
在本说明书一些实施方案中,通过拟合所述中频纵波速度以及所述校正后的低、高频纵波速度,可以获得所述岩芯的全频段速度频散曲线。如图5a~图5d所示,在一示例性实施方案中,通过上述的方式可以得到58号、59号、62号和63号岩芯全频段速度频散曲线。其中,图5a~图5d中,圆点表示对应地层压力下的低频段测量数据,星号表示岩芯深度处测井纵波速度,三角点表示高频测量数据。
S104、根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
在本说明书一些实施方案中,所述根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段,可以包括以下步骤:
1)、确定所述全频段速度频散曲线中,测井频段与地震主频间的纵波速度差平均值。
2)、对所有岩芯取样深度范围内的测井曲线减去所述纵波速度差平均值,从而将所述岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
由此可见,在本说明书实施方案中,由于可以根据岩芯的实际地层压力、实际含流体性,将不同深度取样的岩芯在实验室条件下的低、高频纵波速度进行校正,从而使得校正后的低、高频纵波速度可以与从测井曲线中提取岩芯在对应取芯深度的中频纵波速度,均是在相同的地层压力和含流体性下得到的。在此基础上,就可以根据中频纵波速度以及校正后的低、高频纵波速度生成全频段速度频散曲线,进而可以根据全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段。从而实现在考虑到储层介质速度的频散现象基础上的速度校正,因而有利于获得更为准确的储层介质速度,进而为后续精确地完成井震匹配等工作提供了基础。
参考图6和图7所示,图6和图7中分别示出了转换前后的井震标定结果。对比图6和图7可以看出,图7中基于转换后的井震标定结果得到的时深关系更加合理,井震标定精确性更高,从而了验证考虑频散现象的速度转换方法的有效性。
与上述测井频段速度向地震频段速度转换方法对应,本说明书还提供了测井频段速度向地震频段速度转换装置。参考图8所示,在本说明书一些实施方案中,所述测井频段速度向地震频段速度转换装置可以包括:
获取模块81,可以用于获取不同深度取样的岩芯在实验室条件下的低、高频纵波速度,并从测井曲线中提取所述岩芯在对应取芯深度的中频纵波速度。
校正模块82,可以用于根据所述岩芯的实际地层压力、实际含流体性校正所述低、高频纵波速度,获得校正后的低、高频纵波速度。
生成模块83,可以用于根据所述中频纵波速度以及所述校正后的低、高频纵波速度,生成所述岩芯的全频段速度频散曲线。
转换模块84,可以用于根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
本说明书一些实施方案的测井频段速度向地震频段速度转换装置中,所述根据所述岩芯的实际地层压力校正所述低、高频纵波速度,包括:
根据所述岩芯对应的取样深度及地层压力梯度确定所述岩芯的实际地层压力;
确认所述低、高频纵波速度是否为,在围压为所述实际地层压力下测得的纵波速度;
其中,σ为围压,e为自然常数,为压力系数,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kdry-ini、Gdry-ini分别为干燥岩芯在零围压下的体积模量、剪切模量,Kdry-HP、Gdry-HP分别为干燥岩芯在上极限围压下的体积模量、剪切模量,VP为纵波速度,VS为横波速度;
将所述实际地层压力输入所述对应关系式,获得校正后的低、高频纵波速度。
本说明书一些实施方案的测井频段速度向地震频段速度转换装置中,所述根据所述岩芯的实际含流体性校正所述低、高频纵波速度,包括:
根据测井数据解释成果与录井信息确定所述岩芯的实际含流体性;
确认所述低、高频纵波速度是否为,在含流体性为所述实际含流体性下测得的纵波速度;
其中,Ksat(σ)、Gsat(σ)分别为饱和流体岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kf为饱和流体的体积模量,φ为饱和流体岩芯的孔隙度,Kma为饱和流体岩芯的岩石骨架体积模量,σ为围压,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,VP为纵波速度,VS为横波速度。
本说明书一些实施方案的测井频段速度向地震频段速度转换装置中,所述根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段,包括:
确定所述全频段速度频散曲线中,测井频段与地震主频间的纵波速度差平均值;
对所有岩芯取样深度范围内的测井曲线减去所述纵波速度差平均值,从而将所述岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本说明书时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
与上述测井频段速度向地震频段速度转换方法对应,本说明书还提供了电子设备。参考图9所示,在本说明书一些实施方案中,所述电子设备可以包括存储器、处理器、以及存储在所述存储器上的计算机程序,所述计算机程序被所述处理器运行时执行如下步骤:
获取不同深度取样的岩芯在实验室条件下的低、高频纵波速度,并从测井曲线中提取所述岩芯在对应取芯深度的中频纵波速度;
根据所述岩芯的实际地层压力、实际含流体性校正所述低、高频纵波速度,获得校正后的低、高频纵波速度;
根据所述中频纵波速度以及所述校正后的低、高频纵波速度,生成所述岩芯的全频段速度频散曲线;
根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
虽然上文描述的过程流程包括以特定顺序出现的多个操作,但是,应当清楚了解,这些过程可以包括更多或更少的操作,这些操作可以顺序执行或并行执行(例如使用并行处理器或多线程环境)。
本申请是参照根据本说明书实施方案的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
内存可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。内存是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
本领域技术人员应明白,本说明书的实施方案可提供为方法、系统或计算机程序产品。因此,本说明书实施方案可采用完全硬件实施方案、完全软件实施方案或结合软件和硬件方面的实施方案的形式。而且,本说明书实施方案可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本说明书实施方案可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本说明书实施方案,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施方案均采用递进的方式描述,各个实施方案之间相同相似的部分互相参见即可,每个实施方案重点说明的都是与其他实施方案的不同之处。尤其,对于系统实施方案而言,由于其基本相似于方法实施方案,所以描述的比较简单,相关之处参见方法实施方案的部分说明即可。在本说明书的描述中,参考术语“一个实施方案”、“一些实施方案”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施方案或示例描述的具体特征、结构、材料或者特点包含于本说明书实施方案的至少一个实施方案或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施方案或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施方案或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施方案或示例以及不同实施方案或示例的特征进行结合和组合。
以上所述仅为本申请的实施方案而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。
Claims (10)
1.一种测井频段速度向地震频段速度转换方法,其特征在于,包括:
获取不同深度取样的岩芯在实验室条件下的低、高频纵波速度,并从测井曲线中提取所述岩芯在对应取芯深度的中频纵波速度;
根据所述岩芯的实际地层压力、实际含流体性校正所述低、高频纵波速度,获得校正后的低、高频纵波速度;
根据所述中频纵波速度以及所述校正后的低、高频纵波速度,生成所述岩芯的全频段速度频散曲线;
根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
2.如权利要求1所述的测井频段速度向地震频段速度转换方法,其特征在于,所述根据所述岩芯的实际地层压力校正所述低、高频纵波速度,包括:
根据所述岩芯对应的取样深度及地层压力梯度确定所述岩芯的实际地层压力;
确认所述低、高频纵波速度是否为,在围压为所述实际地层压力下测得的纵波速度;
其中,σ为围压,e为自然常数,为压力系数,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kdry-ini、Gdry-ini分别为干燥岩芯在零围压下的体积模量、剪切模量,Kdry-HP、Gdry-HP分别为干燥岩芯在上极限围压下的体积模量、剪切模量,VP为纵波速度,VS为横波速度;
将所述实际地层压力输入所述对应关系式,获得校正后的低、高频纵波速度。
3.如权利要求1所述的测井频段速度向地震频段速度转换方法,其特征在于,所述根据所述岩芯的实际含流体性校正所述低、高频纵波速度,包括:
根据测井数据解释成果与录井信息确定所述岩芯的实际含流体性;
确认所述低、高频纵波速度是否为,在含流体性为所述实际含流体性下测得的纵波速度;
其中,Ksat(σ)、Gsat(σ)分别为饱和流体岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kf为饱和流体的体积模量,φ为饱和流体岩芯的孔隙度,Kma为饱和流体岩芯的岩石骨架体积模量,σ为围压,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,VP为纵波速度,VS为横波速度。
4.如权利要求1所述的测井频段速度向地震频段速度转换方法,其特征在于,所述根据所述中频纵波速度以及所述校正后的低、高频纵波速度,生成所述岩芯的全频段速度频散曲线,包括:
拟合所述中频纵波速度以及所述校正后的低、高频纵波速度,获得所述岩芯的全频段速度频散曲线。
5.如权利要求1所述的测井频段速度向地震频段速度转换方法,其特征在于,所述根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段,包括:
确定所述全频段速度频散曲线中,测井频段与地震主频间的纵波速度差平均值;
对所有岩芯取样深度范围内的测井曲线减去所述纵波速度差平均值,从而将所述岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
6.一种测井频段速度向地震频段速度转换装置,其特征在于,包括:
获取模块,用于获取不同深度取样的岩芯在实验室条件下的低、高频纵波速度,并从测井曲线中提取所述岩芯在对应取芯深度的中频纵波速度;
校正模块,用于根据所述岩芯的实际地层压力、实际含流体性校正所述低、高频纵波速度,获得校正后的低、高频纵波速度;
生成模块,用于根据所述中频纵波速度以及所述校正后的低、高频纵波速度,生成所述岩芯的全频段速度频散曲线;
转换模块,用于根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
7.如权利要求6所述的测井频段速度向地震频段速度转换装置,其特征在于,所述根据所述岩芯的实际地层压力校正所述低、高频纵波速度,包括:
根据所述岩芯对应的取样深度及地层压力梯度确定所述岩芯的实际地层压力;
确认所述低、高频纵波速度是否为,在围压为所述实际地层压力下测得的纵波速度;
其中,σ为围压,e为自然常数,为压力系数,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kdry-ini、Gdry-ini分别为干燥岩芯在零围压下的体积模量、剪切模量,Kdry-HP、Gdry-HP分别为干燥岩芯在上极限围压下的体积模量、剪切模量,VP为纵波速度,VS为横波速度;
将所述实际地层压力输入所述对应关系式,获得校正后的低、高频纵波速度。
8.如权利要求6所述的测井频段速度向地震频段速度转换装置,其特征在于,所述根据所述岩芯的实际含流体性校正所述低、高频纵波速度,包括:
根据测井数据解释成果与录井信息确定所述岩芯的实际含流体性;
确认所述低、高频纵波速度是否为,在含流体性为所述实际含流体性下测得的纵波速度;
其中,Ksat(σ)、Gsat(σ)分别为饱和流体岩芯在围压σ下的体积模量测量值、剪切模量测量值,Kf为饱和流体的体积模量,φ为饱和流体岩芯的孔隙度,Kma为饱和流体岩芯的岩石骨架体积模量,σ为围压,Kdry(σ)、Gdry(σ)分别为干燥岩芯在围压σ下的体积模量测量值、剪切模量测量值,VP为纵波速度,VS为横波速度。
9.如权利要求6所述的测井频段速度向地震频段速度转换装置,其特征在于,所述根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段,包括:
确定所述全频段速度频散曲线中,测井频段与地震主频间的纵波速度差平均值;
对所有岩芯取样深度范围内的测井曲线减去所述纵波速度差平均值,从而将所述岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
10.一种电子设备,包括存储器、处理器、以及存储在所述存储器上的计算机程序,其特征在于,所述计算机程序被所述处理器运行时执行如下步骤:
获取不同深度取样的岩芯在实验室条件下的低、高频纵波速度,并从测井曲线中提取所述岩芯在对应取芯深度的中频纵波速度;
根据所述岩芯的实际地层压力、实际含流体性校正所述低、高频纵波速度,获得校正后的低、高频纵波速度;
根据所述中频纵波速度以及所述校正后的低、高频纵波速度,生成所述岩芯的全频段速度频散曲线;
根据所述全频段速度频散曲线,将岩芯取样深度范围内的测井频段纵波速度转换至地震频段。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010543114.6A CN111856571B (zh) | 2020-06-15 | 2020-06-15 | 测井频段速度向地震频段速度转换方法、装置及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010543114.6A CN111856571B (zh) | 2020-06-15 | 2020-06-15 | 测井频段速度向地震频段速度转换方法、装置及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111856571A true CN111856571A (zh) | 2020-10-30 |
CN111856571B CN111856571B (zh) | 2022-04-08 |
Family
ID=72986623
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010543114.6A Active CN111856571B (zh) | 2020-06-15 | 2020-06-15 | 测井频段速度向地震频段速度转换方法、装置及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111856571B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060285437A1 (en) * | 2005-06-03 | 2006-12-21 | Schlumberger Technology Corporation | Radial profiling of formation mobility using horizontal and vertical shear slowness profiles |
CN102445709A (zh) * | 2010-10-14 | 2012-05-09 | 中国石油大学(北京) | 与孔隙结构相关的全频带速度预测模型 |
CN102508296A (zh) * | 2011-11-14 | 2012-06-20 | 中国石油天然气股份有限公司 | 一种非饱和双重孔隙介质地震波频散衰减分析方法及装置 |
CN102854531A (zh) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | 用于储层烃类检测的多尺度岩石物理图版方法及装置 |
CN104656142A (zh) * | 2013-11-19 | 2015-05-27 | 中国石油天然气集团公司 | 一种利用垂直地震剖面与测井联合的地震层位标定方法 |
CN107065013A (zh) * | 2016-12-27 | 2017-08-18 | 中国石油大学(华东) | 一种地震尺度下的层速度确定方法及装置 |
CN108562938A (zh) * | 2018-03-23 | 2018-09-21 | 中国石油天然气股份有限公司 | 一种消除频散效应的方法、装置及系统 |
CN110261898A (zh) * | 2019-04-26 | 2019-09-20 | 中国石油化工股份有限公司 | 基于地震岩石物理实验分析的测井与地震速度匹配方法 |
CN110333549A (zh) * | 2019-04-26 | 2019-10-15 | 中国石油化工股份有限公司 | 基于测井资料频散特征的井点裂隙孔隙度预测方法 |
-
2020
- 2020-06-15 CN CN202010543114.6A patent/CN111856571B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060285437A1 (en) * | 2005-06-03 | 2006-12-21 | Schlumberger Technology Corporation | Radial profiling of formation mobility using horizontal and vertical shear slowness profiles |
CN102445709A (zh) * | 2010-10-14 | 2012-05-09 | 中国石油大学(北京) | 与孔隙结构相关的全频带速度预测模型 |
CN102508296A (zh) * | 2011-11-14 | 2012-06-20 | 中国石油天然气股份有限公司 | 一种非饱和双重孔隙介质地震波频散衰减分析方法及装置 |
CN102854531A (zh) * | 2012-09-11 | 2013-01-02 | 中国石油天然气股份有限公司 | 用于储层烃类检测的多尺度岩石物理图版方法及装置 |
CN104656142A (zh) * | 2013-11-19 | 2015-05-27 | 中国石油天然气集团公司 | 一种利用垂直地震剖面与测井联合的地震层位标定方法 |
CN107065013A (zh) * | 2016-12-27 | 2017-08-18 | 中国石油大学(华东) | 一种地震尺度下的层速度确定方法及装置 |
CN108562938A (zh) * | 2018-03-23 | 2018-09-21 | 中国石油天然气股份有限公司 | 一种消除频散效应的方法、装置及系统 |
CN110261898A (zh) * | 2019-04-26 | 2019-09-20 | 中国石油化工股份有限公司 | 基于地震岩石物理实验分析的测井与地震速度匹配方法 |
CN110333549A (zh) * | 2019-04-26 | 2019-10-15 | 中国石油化工股份有限公司 | 基于测井资料频散特征的井点裂隙孔隙度预测方法 |
Non-Patent Citations (6)
Title |
---|
MA XIAO-YI 等: "Velocity dispersion and fluid substitution in sandstone under partially saturated conditions", 《APPLIED GEOPHYSICS》 * |
刘志远 等: "基于饱水页岩地震频段岩石物理测试的地震响应频散特征分析", 《物探化探计算技术》 * |
张程远 等: "基于岩体完整性评价的超声–地震波速度跨尺度转换方法研究", 《岩石力学与工程学报》 * |
赵建国 等: "跨频段岩石物理测量技术进展及在储层预测中的应用", 《中国地球科学联合学术年会2017》 * |
邓继新 等: "基于储层砂岩微观孔隙结构特征的弹性波频散响应分析", 《地球物理学报》 * |
陆永杰 等: "测井与地震数据匹配方法研究", 《CPS/SEG北京2018国际地球物理会议暨展览》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111856571B (zh) | 2022-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11098565B2 (en) | Method for estimating permeability of fractured rock formations from induced slow fluid pressure waves | |
US6954402B2 (en) | Method of wave diagnostics of the oil-and-gas-deposit | |
US10921477B2 (en) | Method, apparatus, and system for eliminating frequency dispersion effect | |
CN103064112B (zh) | 一种检测地层含油气性的吸收方法 | |
Huang et al. | Improving production history matching using time-lapse seismic data | |
Kianoush et al. | Compressional and shear interval velocity modeling to determine formation pressures in an oilfield of SW Iran | |
CN113050189B (zh) | 一种测井曲线的重构方法、装置、设备及存储介质 | |
CN106355571A (zh) | 一种白云岩储层质量的确定方法及装置 | |
CN108873065B (zh) | 砂岩优质储层预测方法及装置 | |
CN111208565B (zh) | 基于kt模型的孔缝参数反演方法、装置及存储介质 | |
CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
CN111812716B (zh) | 页岩气储层总有机碳含量叠前定量预测方法、装置及设备 | |
Liu et al. | Combined petrophysics and 3D seismic attributes to predict shale reservoirs favourable areas | |
CN103628866A (zh) | 获取地层的地层压力系数的方法 | |
CN111208566B (zh) | 基于sca模型的孔缝参数反演方法、装置及存储介质 | |
CN111856571B (zh) | 测井频段速度向地震频段速度转换方法、装置及设备 | |
CN111239809B (zh) | 一种基于属性融合原理的优势储层确定方法及系统 | |
CN113219531A (zh) | 致密砂岩气水分布的识别方法及装置 | |
CN112861061A (zh) | 地质分层方法及系统、处理设备 | |
CN111502647A (zh) | 一种钻井地质环境因素确定方法、装置及存储介质 | |
CN112305594A (zh) | 非均质储层的油气分布确定方法及系统 | |
CN112180443A (zh) | 页岩气二维地震甜点区优选方法及装置 | |
CN108931813B (zh) | 一种识别礁盖、礁核、礁基的方法 | |
CN107589457A (zh) | 测井资料尺度粗化的方法和装置 | |
CN113253346B (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 |