CN111456726A - 声波时差校正方法及装置 - Google Patents
声波时差校正方法及装置 Download PDFInfo
- Publication number
- CN111456726A CN111456726A CN202010380368.0A CN202010380368A CN111456726A CN 111456726 A CN111456726 A CN 111456726A CN 202010380368 A CN202010380368 A CN 202010380368A CN 111456726 A CN111456726 A CN 111456726A
- Authority
- CN
- China
- Prior art keywords
- point
- time difference
- diameter
- acoustic
- diameter expansion
- 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
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing 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
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
- E21B47/12—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
- E21B47/14—Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling using acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
本发明提供了一种声波时差校正方法及装置,其中该方法包括:根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点;根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值;根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子;根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段的声波时差校正值。该方法得到了更加贴近实际的声波时差校正值,提高了声波时差的精度,降低了测井所得的声波时差与实际的误差。
Description
技术领域
本发明涉及石油勘探技术领域,尤其涉及一种声波时差校正方法及装置。
背景技术
声波时差是井震标定、岩石物理建模、横波预测、正演和反演技术的必备参数。声波时差测井曲线的正确与否对于地震合成记录,标定层位、地震反演有重要的影响。但测井资料在纵向上具有较高的分辨率,但是质量易受测井环境因素(测井仪器,井径,泥浆、地层温度和压力等)的影响,偏离原状地层真实值,造成数据畸变失真,从而影响储层预测的准确性。裸眼井声波测井测量是由发射器和接收器间沿井壁传播的滑行纵波的时间差。因此,井眼垮塌、井壁凹凸不平、仪器倾斜等都会给声波时差测井曲线带来误差。
目前常用的双发双收声波补偿测井,能够对仪器倾斜、洞穴或者井眼倾斜造成的测量误差起到明显的补偿作用。但是当井径扩径或不规则严重时,声波时差的异常增大不可避免。目前常用的校正方法分为两大类:
第一大类为需要声波时差和井径以外其他曲线参与的校正方法,如声波时差上限约束法;Gardner公式和Faust公式转换法,多参数拟合法,岩石物理模型重构法以及近年来发展迅猛的人工智能等方法。以上常用方法建立和验证非扩径处声波时差与其他曲线的关系,以此关系为桥梁,利用扩径处的其他曲线,重构扩径处声波时差。扩径校正的正确性依赖于参与计算的其他测井曲线在扩径段的正确性,但是无论密度、电阻率或者中子等曲线,扩径时的测量值均受到扩径处填充的泥浆的影响,使用失真的测井曲线来反推声波时差值,精度低且误差高。
第二类为仅有声波时差和井径参与的校正方法,如路径校正法和漫反射声路法、频率法等,但是由于此类方法常常只侧重于扩径时泥浆的影响,忽视了井径的不规则程度对测量带来的影响,导致得到的声波时差精度低且误差高。
发明内容
本发明实施例提供一种声波时差校正方法,用以提高声波时差的精度,降低测井所得的声波时差的误差,该方法包括:
获取解释层段中每个深度点处的声波时差测量值、泥质含量、测井直径、测井所用钻头直径以及仪器测量数据;
根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点;
根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值;
根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子;
根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段的声波时差校正值。
本发明实施例还提供一种声波时差校正装置,用以提高声波时差的精度,降低测井所得的声波时差的误差,该装置包括:
数据获取模块,用于获取解释层段中每个深度点处的声波时差测量值、泥质含量、测井直径、测井所用钻头直径以及仪器测量数据;
扩径确定模块,用于根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点;
第一校正模块,用于根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值;
影响因子确定模块,用于根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子;
第二校正模块,用于根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段的声波时差校正值。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述声波时差校正方法。
本发明实施例也提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述声波时差校正方法的计算机程序。
本发明实施例中,通过获取解释层段中每个深度点处的声波时差测量值、泥质含量、测井直径、测井所用钻头直径以及仪器测量数据;根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点;根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值;根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子;根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段的声波时差校正值。实现了无需井径以外其他曲线参与,仅需要声波时差和井径参与的校正;通过增加扩径点处的井壁不规则影响因子和泥浆路径影响因子,在声波时差校正时,既考虑了井径的不规则程度对测量带来的影响,也考虑了井眼不规则严重时泥浆和纵波传播对测量值的影响,从而得到更加贴近实际的声波时差校正值,提高了声波时差的精度,降低了测井所得的声波时差与实际的误差。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中声波时差校正方法的示意图。
图2为本发明实施例中步骤104的具体实施方法示意图。
图3为本发明实施例中步骤105的具体实施方法示意图。
图4为本发明一具体应用实施中测井的部分层段声波时差校正前后与合成记录对比图。
图5为本发明实施例中声波时差校正装置的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供了一种声波时差校正方法,用以提高声波时差的精度,降低测井所得的声波时差的误差,如图1所示,该方法包括:
步骤101:获取解释层段中每个深度点处的声波时差测量值、泥质含量、测井直径、测井所用钻头直径以及仪器测量数据;
步骤102:根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点;
步骤103:根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值;
步骤104:根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子;
步骤105:根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段的声波时差校正值。
由图1所示流程可以得知,本发明实施例中,通过获取解释层段中每个深度点处的声波时差测量值、泥质含量、测井直径、测井所用钻头直径以及仪器测量数据;根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点;根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值;根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子;根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段的声波时差校正值。实现了无需井径以外其他曲线参与,仅需要声波时差和井径参与的校正;通过增加扩径点处的井壁不规则影响因子和泥浆路径影响因子,在声波时差校正时,既考虑了井径的不规则程度对测量带来的影响,也考虑了井眼不规则严重时泥浆和纵波传播对测量值的影响,从而得到更加贴近实际的声波时差校正值,提高了声波时差的精度,降低了测井所得的声波时差与实际的误差。
具体实施时,首先获取解释层段中每个深度点处的声波时差测量值、泥质含量、测井直径、测井所用钻头直径以及仪器测量数据。其中,解释层段在声波时差测量时被分为若干个深度点。仪器测量数据例如可以包括:发射器间距、接收器间距和仪器的垂直分辨率。
根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点,具体过程包括:若一深度点的测井直径和测井所用钻头直径之间的差值的绝对值,大于或等于预设误差,则确定该深度点为扩径点;若一深度点的测井直径和测井所用钻头直径之间的差值的绝对值,小于预设误差,则确定该深度点为标准井径点。预设误差根据实际需求设置,例如可以为0.001cm,本领域技术人员可以理解,上述预设误差值仅为举例,不用于限定本发明的保护范围。
确定出所有深度点中的扩径点后,根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值。具体实施时,根据泥质含量最大的标准井径点处的声波时差测量值、孔隙度最大的标准井径点处的声波时差测量值以及每个扩径点处的泥质含量,确定每个扩径点处的声波时差上限值。具体按照如下公式,确定每个扩径点处的声波时差上限值:
Δtmax=Vsh×Δtsh+(1-Vsh)×Δtp (1)
其中,Δtmax表示扩径点处的声波时差上限值;Vsh表示扩径点处的泥质含量;
Δtsh表示泥质含量最大的标准井径点处的声波时差测量值;
Δtp表示孔隙度最大的标准井径点处的声波时差测量值。
确定每个扩径点处的声波时差上限值后,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值,具体过程,包括:
若一扩径点处的声波时差测量值大于声波时差上限值,则该扩径点处的声波时差的校正测量值为该扩径点处的声波时差上限值;
若一扩径点处的声波时差测量值小于或等于声波时差上限值,则该扩径点处的声波时差的校正测量值为该扩径点处的声波时差测量值。
得到每个扩径点处的声波时差的校正测量值后,根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子。具体过程如图2所示,包括:
步骤201:根据每个扩径点处的发射器间距和垂直分辨率,确定每个扩径点处的接收器间滑行纵波的传播距离;
步骤202:根据每个扩径点处的接收器间滑行纵波的传播距离和接收器间距,确定每个扩径点处的井壁不规则影响因子;
步骤203:根据每个扩径点处的接收器间距,确定每个扩径点处的接收器在泥浆路径中接受声波的差异系数;
步骤204:根据差异系数和声波在泥浆中传播的时差,确定每个扩径点处的泥浆路径影响因子。
其中,步骤202具体实施时,按照如下公式,确定每个扩径点处的井壁不规则影响因子:
ratio_enlarge=L/LRR (2)
其中,ratio_enlarge表示扩径点处的井壁不规则影响因子;L表示扩径点处的接收器间滑行纵波的传播距离;LRR表示扩径点处的接收器间距。
步骤204具体实施时,按照如下公式,确定每个扩径点处的泥浆路径影响因子:
Delta_DTmud=Δtmud×ΔLmud (3)
其中,Delta_DTmud表示扩径点处的泥浆路径影响因子;Δtmud表示扩径点处的声波在泥浆中传播的时差;ΔLmud表示扩径点处的接收器在泥浆路径中接受声波的差异系数,可以根据理论的接收器间距LRR内的扩径程度和不规则程度共同确定。
确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子后,根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段的声波时差校正值,具体过程如图3所示,包括:
步骤301:根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定每个扩径点处的声波时差校正值;
步骤302:整合每个扩径点处的声波时差校正值和每个标准井径点处的声波时差测量值,得到解释层段地层的声波时差校正值。
步骤301具体实施时,按照如下公式,确定每个扩径点处的声波时差校正值:
Δt_cor=(Δt-Delta_DTmud)/ratio_enlarge (4)
其中,Δt_cor表示扩径点处的声波时差校正值;Delta_DTmud表示扩径点处的泥浆路径影响因子;ratio_enlarge表示扩径点处的泥浆路径影响因子;Δt表示扩径点处的声波时差的校正测量值。
下面给出一具体实例说明本发明实施例如何进行声波时差校正。本例应用于一具体测井。
首先统计待校正地区测井的声波时差、泥质含量和测井直径、测井所用钻头直径等测量信息和仪器发射器、接收器间距、垂直分辨率等仪器信息。
判断井眼是否扩径,井眼直径Cal,钻头直径为Bits,给定误差限ε,当某个深度点的|Cal-Bits|<ε,则标记为标准井径点;否则,标记为扩径点。
计算每个扩径点处的地层声波时差上限Δtmax。当扩径点处的声波时差测量值小于Δtmax时,扩径点处的声波时差的校正测量值确定为声波时差测量值;当声波时差测量值大于Δtmax,扩径点处的声波时差的校正测量值等于Δtmax。
补偿声波测量为接收器间滑行纵波的平均时差,井壁的不规则度越高,测井直径的变化越剧烈,纵波沿地层的实际传播距离越大。因此,在不规则严重的地层,滑行纵波实际的传播距离L要大于理论的接收器间距LRR。通过实际滑行距离L和理论间距LRR利用公式(2)衡量井壁的不规则程度ratio_enlarge。
扩径或井眼形状不规则严重时,补偿声波测井的补偿作用有限,声波从地层穿过泥浆返回接收器的距离差距较大,通过公式(3)计算扩径点处的泥浆路径影响因子Delta_DTmud。
声波时差为声波通过单位距离的传播时间,通过公式(4)得到扩径点处的声波时差校正值Δt_cor。若该深度点为标准井径点,声波时差校正值即为声波时差测量值。
对解释层段中的每一个深度点,重复上述步骤,得到每个深度点的声波时差校正值,根据每个深度点的声波时差校正值,完成解释层段内的声波时差校正曲线。
图4所示为此具体实例中测井的部分层段声波时差校正前后与合成记录对比图。图中第1道和第2道分别是原始的测井直径(单位cm)、声波时差测量值的曲线(单位us/m),可以观察到大部分情况下,井中扩径十分严重、测井直径变化剧烈时,扩径段声波时差测量值出现明显增大。第3道为井眼坍塌时泥浆对声波时差测量值的影响(单位us/m),总体来看,井径的扩径程度与校正量(黑色实线)成正比;第4道为井径不规则程度(黑色虚线,无量纲),值越高表示井壁凹凸不平越严重;第5道为不规则程度对于测量值的影响(黑色实线,单位us/m),二者成正比关系。第6道为校正后的声波时差校正值(黑色实线)与校正前的声波时差测量值(黑色虚线),第7道为不同测井直径下声波时差校正量(单位us/m)。第8道为测井曲线对应的深度/时间。可以看到,扩径严重的位置,校正前后声波时差出现明显变化,并且相同测井直径下,井壁不规则化严重地层,声波时差校正量大于扩径而平缓的层段的校正量。第9、10和11道分别为声波时差校正前、校正后的合成记录和井旁道的对比,校正后的合成记录与井旁道的吻合程度明显提高,证明了利用本申请声波时差校正方法后得到的校正值的合理性。
通过分析声波时差与测井直径变化的关系,联合井壁不规则程度评估扩径处的滑行纵波沿地层和泥浆的传播路径的影响,消除井径变化对声波时差测井的影响,得到校正后的真实声波时差。同时考虑了扩径或井眼不规则严重时泥浆和纵波传播对测量值的影响,且不需要除测井直径以外其他曲线的参与,无须在校正前保证扩径处其他测井参数的正确性,应用更加灵活广泛。提高了声波时差测量的正确性,在井震标定、岩石物理建模、正演和反演方面具有推广价值。
上述具体应用的实施仅为举例,其余实施方式不再一一赘述。
基于同一发明构思,本发明实施例还提供一种声波时差校正装置,由于声波时差校正装置所解决问题的原理与声波时差校正方法相似,因此声波时差校正装置的实施可以参见声波时差校正方法的实施,重复之处不再赘述,具体结构如图5所示:
数据获取模块501,用于获取解释层段中每个深度点处的声波时差测量值、泥质含量、测井直径、测井所用钻头直径以及仪器测量数据;
扩径确定模块502,用于根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点;
第一校正模块503,用于根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值;
影响因子确定模块504,用于根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子;
第二校正模块505,用于根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段的声波时差校正值。
具体实施例中,扩径确定模块502具体用于:
若一深度点的测井直径和测井所用钻头直径之间的差值的绝对值,大于或等于预设误差,则确定所述深度点为扩径点;
若一深度点的测井直径和测井所用钻头直径之间的差值的绝对值,小于预设误差,则确定所述深度点为标准井径点。
具体实施例中,第一校正模块503具体用于:
根据泥质含量最大的标准井径点处的声波时差测量值、孔隙度最大的标准井径点处的声波时差测量值以及每个扩径点处的泥质含量,确定每个扩径点处的声波时差上限值。
具体实施例中,第一校正模块503具体用于:
若一扩径点处的声波时差测量值大于声波时差上限值,则该扩径点处的声波时差的校正测量值为该扩径点处的声波时差上限值;
若一扩径点处的声波时差测量值小于或等于声波时差上限值,则该扩径点处的声波时差的校正测量值为该扩径点处的声波时差测量值。
具体实施时,仪器测量数据包括:发射器间距、接收器间距和仪器的垂直分辨率。影响因子确定模块504具体用于:
根据每个扩径点处的发射器间距和垂直分辨率,确定每个扩径点处的接收器间滑行纵波的传播距离;
根据每个扩径点处的接收器间滑行纵波的传播距离和接收器间距,确定每个扩径点处的井壁不规则影响因子;
根据每个扩径点处的接收器间距,确定每个扩径点处的接收器在泥浆路径中接受声波的差异系数;
根据上述差异系数和声波在泥浆中传播的时差,确定每个扩径点处的泥浆路径影响因子。
具体实施例中,第二校正模块505具体用于:
根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定每个扩径点处的声波时差校正值;
整合每个扩径点处的声波时差校正值和每个标准井径点处的声波时差测量值,得到解释层段地层的声波时差校正值。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述声波时差校正方法。
本发明实施例还提供一种计算机可读存储介质,存储有执行上述声波时差校正方法的计算机程序。
综上所述,本发明实施例提供的声波时差校正方法及装置具有如下优点:
通过获取解释层段中每个深度点处的声波时差测量值、泥质含量、测井直径、测井所用钻头直径以及仪器测量数据;根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点;根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值;根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子;根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段的声波时差校正值。实现了无需井径以外其他曲线参与,仅需要声波时差和井径参与的校正;通过增加扩径点处的井壁不规则影响因子和泥浆路径影响因子,在声波时差校正时,既考虑了井径的不规则程度对测量带来的影响,也考虑了井眼不规则严重时泥浆和纵波传播对测量值的影响,从而得到更加贴近实际的声波时差校正值,提高了声波时差的精度,降低了测井所得的声波时差与实际的误差。且校正时不需要其他地质参量测井曲线的参与,无须在校正前保证扩径处其他测井参数的正确性,从而应用更加灵活广泛。
本领域内的技术人员应明白,本发明的实施例可提供为方法、装置、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、装置和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明实施例可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (17)
1.一种声波时差校正方法,其特征在于,包括:
获取解释层段中每个深度点处的声波时差测量值、泥质含量、测井直径、测井所用钻头直径以及仪器测量数据;
根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点;
根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值;
根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子;
根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段的声波时差校正值。
2.如权利要求1所述的方法,其特征在于,根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点,包括:
若一深度点的测井直径和测井所用钻头直径之间的差值的绝对值,大于或等于预设误差,则确定所述深度点为扩径点;
若一深度点的测井直径和测井所用钻头直径之间的差值的绝对值,小于预设误差,则确定所述深度点为标准井径点。
3.如权利要求2所述的方法,其特征在于,根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值,包括:
根据泥质含量最大的标准井径点处的声波时差测量值、孔隙度最大的标准井径点处的声波时差测量值以及每个扩径点处的泥质含量,确定每个扩径点处的声波时差上限值。
4.如权利要求3所述的方法,其特征在于,按照如下公式,根据泥质含量最大的标准井径点处的声波时差测量值、孔隙度最大的标准井径点处的声波时差测量值以及每个扩径点处的泥质含量,确定每个扩径点处的声波时差上限值:
Δtmax=Vsh×Δtsh+(1-Vsh)×Δtp
其中,Δtmax表示扩径点处的声波时差上限值;Vsh表示扩径点处的泥质含量;
Δtsh表示泥质含量最大的标准井径点处的声波时差测量值;
Δtp表示孔隙度最大的标准井径点处的声波时差测量值。
5.如权利要求2所述的方法,其特征在于,根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段地层的声波时差校正值,包括:
根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定每个扩径点处的声波时差校正值;
整合每个扩径点处的声波时差校正值和每个标准井径点处的声波时差测量值,得到解释层段地层的声波时差校正值。
6.如权利要求5所述的方法,其特征在于,按照如下公式,根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定每个扩径点处的声波时差校正值:
Δt_cor=(Δt-Delta_DTmud)/ratio_enlarge
其中,Δt_cor表示扩径点处的声波时差校正值;Delta_DTmud表示扩径点处的泥浆路径影响因子;ratio_enlarge表示扩径点处的泥浆路径影响因子;Δt表示扩径点处的声波时差的校正测量值。
7.如权利要求1所述的方法,其特征在于,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值,包括:
若一扩径点处的声波时差测量值大于声波时差上限值,则所述扩径点处的声波时差的校正测量值为所述扩径点处的声波时差上限值;
若一扩径点处的声波时差测量值小于或等于声波时差上限值,则所述扩径点处的声波时差的校正测量值为所述扩径点处的声波时差测量值。
8.如权利要求1所述的方法,其特征在于,所述仪器测量数据包括:
发射器间距、接收器间距和仪器的垂直分辨率。
9.如权利要求8所述的方法,其特征在于,根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子,包括:
根据每个扩径点处的发射器间距和垂直分辨率,确定每个扩径点处的接收器间滑行纵波的传播距离;
根据每个扩径点处的接收器间滑行纵波的传播距离和接收器间距,确定每个扩径点处的井壁不规则影响因子;
根据每个扩径点处的接收器间距,确定每个扩径点处的接收器在泥浆路径中接受声波的差异系数;
根据所述差异系数和声波在泥浆中传播的时差,确定每个扩径点处的泥浆路径影响因子。
10.如权利要求9所述的方法,其特征在于,按照如下公式,根据每个扩径点处的接收器间滑行纵波的传播距离和接收器间距,确定每个扩径点处的井壁不规则影响因子:
ratio_enlarge=L/LRR
其中,ratio_enlarge表示扩径点处的井壁不规则影响因子;L表示扩径点处的接收器间滑行纵波的传播距离;LRR表示扩径点处的接收器间距。
11.如权利要求9所述的方法,其特征在于,按照如下公式,根据所述差异系数和声波在泥浆中传播的时差,确定每个扩径点处的泥浆路径影响因子:
Delta_DTmud=Δtmud×ΔLmud
其中,Delta_DTmud表示扩径点处的泥浆路径影响因子;Δtmud表示扩径点处的声波在泥浆中传播的时差;ΔLmud表示扩径点处的接收器在泥浆路径中接受声波的差异系数。
12.一种声波时差校正装置,其特征在于,包括:
数据获取模块,用于获取解释层段中每个深度点处的声波时差测量值、泥质含量、测井直径、测井所用钻头直径以及仪器测量数据;
扩径确定模块,用于根据每个深度点处的测井直径和测井所用钻头直径,确定每个深度点是否为扩径点;
第一校正模块,用于根据每个扩径点处的泥质含量确定每个扩径点处的声波时差上限值,根据每个扩径点处的声波时差上限值和声波时差测量值,得到每个扩径点处的声波时差的校正测量值;
影响因子确定模块,用于根据每个扩径点处的仪器测量数据,确定每个扩径点处的井壁不规则影响因子和泥浆路径影响因子;
第二校正模块,用于根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定解释层段的声波时差校正值。
13.如权利要求12所述的装置,其特征在于,所述扩径确定模块具体用于:
若一深度点的测井直径和测井所用钻头直径之间的差值的绝对值,大于或等于预设误差,则确定所述深度点为扩径点;
若一深度点的测井直径和测井所用钻头直径之间的差值的绝对值,小于预设误差,则确定所述深度点为标准井径点。
14.如权利要求13所述的装置,其特征在于,所述第一校正模块具体用于:
根据泥质含量最大的标准井径点处的声波时差测量值、孔隙度最大的标准井径点处的声波时差测量值以及每个扩径点处的泥质含量,确定每个扩径点处的声波时差上限值。
15.如权利要求13所述的装置,其特征在于,所述第二校正模块具体用于:
根据每个扩径点处的声波时差的校正测量值、井壁不规则影响因子和泥浆路径影响因子,确定每个扩径点处的声波时差校正值;
整合每个扩径点处的声波时差校正值和每个标准井径点处的声波时差测量值,得到解释层段地层的声波时差校正值。
16.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至11任一所述方法。
17.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至11任一所述方法的计算机程序。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010380368.0A CN111456726B (zh) | 2020-05-08 | 2020-05-08 | 声波时差校正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010380368.0A CN111456726B (zh) | 2020-05-08 | 2020-05-08 | 声波时差校正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111456726A true CN111456726A (zh) | 2020-07-28 |
CN111456726B CN111456726B (zh) | 2023-06-30 |
Family
ID=71677783
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010380368.0A Active CN111456726B (zh) | 2020-05-08 | 2020-05-08 | 声波时差校正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111456726B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112228037A (zh) * | 2020-09-04 | 2021-01-15 | 长江大学 | 一种基于声波传播理论的孔隙度解释方法和装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1053036A1 (ru) * | 1981-07-14 | 1983-11-07 | Всесоюзный научно-исследовательский институт нефтепромысловой геофизики | Способ поверки аппаратуры акустического каротажа |
CN102454399A (zh) * | 2010-10-26 | 2012-05-16 | 中国石油化工股份有限公司 | 测井声波时差信号校正方法 |
CN104453873A (zh) * | 2014-10-23 | 2015-03-25 | 中国石油天然气股份有限公司 | 页岩油气经济有效层段的评价方法 |
CN107092031A (zh) * | 2017-06-01 | 2017-08-25 | 中国石油天然气集团公司 | 一种声波测井数据校正方法及装置 |
-
2020
- 2020-05-08 CN CN202010380368.0A patent/CN111456726B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1053036A1 (ru) * | 1981-07-14 | 1983-11-07 | Всесоюзный научно-исследовательский институт нефтепромысловой геофизики | Способ поверки аппаратуры акустического каротажа |
CN102454399A (zh) * | 2010-10-26 | 2012-05-16 | 中国石油化工股份有限公司 | 测井声波时差信号校正方法 |
CN104453873A (zh) * | 2014-10-23 | 2015-03-25 | 中国石油天然气股份有限公司 | 页岩油气经济有效层段的评价方法 |
CN107092031A (zh) * | 2017-06-01 | 2017-08-25 | 中国石油天然气集团公司 | 一种声波测井数据校正方法及装置 |
Non-Patent Citations (4)
Title |
---|
刘国强;谭廷栋;: "声波和密度测井资料的井径校正" * |
刘浩杰;王延光;王慧;: "基于井径变化的声波测井时差校正" * |
吉伟平;王志章;姜爽;曾经;: "井眼垮塌条件下声波和密度曲线的校正方法" * |
王莹,陈兆山,麻凤海: "声波测井曲线环境校正方法探讨" * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112228037A (zh) * | 2020-09-04 | 2021-01-15 | 长江大学 | 一种基于声波传播理论的孔隙度解释方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN111456726B (zh) | 2023-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7840394B2 (en) | Method for generating a 3D earth model | |
CN106405651B (zh) | 一种基于测井匹配的全波形反演初始速度模型构建方法 | |
CN111520126A (zh) | 密度测井扩径校正方法及装置 | |
CN110703330B (zh) | 基于岩石物理的碳酸盐岩测井校正和横波速度预测方法 | |
CN111650646B (zh) | 一种用于地震处理中的声波曲线校正方法 | |
CN106548456A (zh) | 校正超声图像轮廓中的偏心的方法 | |
CN111208564A (zh) | 一种深度域层位标定方法及装置 | |
CN111456726A (zh) | 声波时差校正方法及装置 | |
CN111308558B (zh) | 页岩气水平井纵波时差校正方法 | |
CN114428355A (zh) | 储层物性参数直接反演方法、装置、电子设备及介质 | |
CN113189672B (zh) | 一种基于多属性反演的隧道超前地质预报方法 | |
CN109061737B (zh) | 一种基于合成地震记录的储层预测方法及装置 | |
CN110219642B (zh) | 基于声波传播路径的声波时差校正方法 | |
CN113267809A (zh) | I类页岩储层预测方法及装置 | |
CN112906465A (zh) | 一种基于地层因素的煤系地层声波曲线重构方法及系统 | |
CN112859172A (zh) | 纵横波微测井资料处理方法及装置 | |
CN105178950B (zh) | 一种确定细粒沉积岩中总有机碳的方法 | |
CN113187457B (zh) | 一种确定储层压裂效果的方法、装置和存储介质 | |
CN116500686A (zh) | 复杂井径扩径情况下的声波曲线校正方法 | |
CN114893166B (zh) | 一种地层压力系数计算方法 | |
CN112228046B (zh) | 一种随钻超声井径数据校正方法 | |
CN111198401B (zh) | 基于多地球物理参数相关性的测井曲线校正方法与装置 | |
CN113466943B (zh) | 变偏移距vsp地表一致性振幅补偿方法和装置 | |
CN114063193A (zh) | 基于lld和gr联合刻度的致密砂岩密度测井曲线校正方法 | |
CN113009579B (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 |