CN113884580A - 激波管实验斯通利波幅度提取方法及装置 - Google Patents

激波管实验斯通利波幅度提取方法及装置 Download PDF

Info

Publication number
CN113884580A
CN113884580A CN202010625075.4A CN202010625075A CN113884580A CN 113884580 A CN113884580 A CN 113884580A CN 202010625075 A CN202010625075 A CN 202010625075A CN 113884580 A CN113884580 A CN 113884580A
Authority
CN
China
Prior art keywords
waveform data
amplitude
waveform
time
stoneley wave
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
CN202010625075.4A
Other languages
English (en)
Other versions
CN113884580B (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 CN202010625075.4A priority Critical patent/CN113884580B/zh
Publication of CN113884580A publication Critical patent/CN113884580A/zh
Application granted granted Critical
Publication of CN113884580B publication Critical patent/CN113884580B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/4463—Signal correction, e.g. distance amplitude correction [DAC], distance gain size [DGS], noise filtering
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/4409—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
    • G01N29/449—Statistical methods not provided for in G01N29/4409, e.g. averaging, smoothing and interpolation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/023—Solids
    • G01N2291/0232—Glass, ceramics, concrete or stone
    • 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
    • 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
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Signal Processing (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Probability & Statistics with Applications (AREA)
  • Acoustics & Sound (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明提供了一种激波管实验斯通利波幅度提取方法及装置，所述方法包含：获取激波管实验过程中激波管激励声波后水下的第一波形数据和岩心侧面的第二波形数据；对所述第一波形数据进行幅度校正使波形基线归零；对所述第二波形数据和幅度校正后的所述第一波形数据进行幅度一致性处理后进行第一次时间校正，使波形起始时间为波形到达岩心顶端时间；利用斯通利波速度对第一次时间校正后的所述第二波形数据进行第二次时间校正，使所述第二波形数据的多道接收波形中斯通利波到达时间一致；通过第二次时间校正后的所述第二波形数据获取斯通利波幅度。

Description

激波管实验斯通利波幅度提取方法及装置
技术领域
本发明涉及裂缝型碳酸盐岩储层油气勘探领域，针对激波管实验装置对岩心的测量信号，提出了一种激波管实验斯通利波幅度提取方法及装置。
背景技术
在碳酸盐岩、火山岩等裂缝孔洞型储层中，基质渗透率往往很低，对油气流动起主导作用的是各种裂缝，因此，裂缝，特别是裂缝宽度的准确识别及评价对复杂储层油气勘探具有重要意义，这也是目前测井储层评价中面临的主要挑战之一。
沿着井壁传播的斯通利波经过含裂缝储层时，因部分能量反射、部分能量沿裂缝进入储层而使得斯通利波幅度明显衰减，因此可利用斯通利波定量评价裂缝。激波管实验装置能够测量得到纯净的斯通利波信号，为斯通利波衰减规律研究提供了一种可靠的物理实验手段。用激波管实验研究裂缝对斯通利波的影响特征，需对含裂缝岩心进行多次激波管实验，获取多组实验波形数据。然而，实验获取的原始波形中存在基线漂移、幅度不一致性、斯通利波到达时间不一致等多种干扰因素，如何从原始波形中提取准确的斯通利波幅度是一个难题，同时也具有重要意义。
发明内容
本发明目的在于提供一种激波管实验斯通利波幅度提取方法及装置，为利用斯通利波定量评价岩石裂缝渗透率奠定了基础。
为达上述目的，本发明所提供的激波管实验斯通利波幅度提取方法，具体包含：获取激波管实验过程中激波管激励声波后水下的第一波形数据和岩心侧面的第二波形数据；对所述第一波形数据进行幅度校正使波形基线归零；对所述第二波形数据和幅度校正后的所述第一波形数据进行幅度一致性处理后进行第一次时间校正，使波形起始时间为波形到达岩心顶端时间；利用斯通利波速度对第一次时间校正后的所述第二波形数据进行第二次时间校正，使所述第二波形数据的多道接收波形中斯通利波到达时间一致；
通过第二次时间校正后的所述第二波形数据获取斯通利波幅度。
在上述激波管实验斯通利波幅度提取方法中，优选的，对所述第一波形数据进行幅度校正使波形基线归零包含：计算所述第一波形数据中的噪声段的幅度均值，根据所述幅度均值修正波形数据使波形基线归零。
在上述激波管实验斯通利波幅度提取方法中，优选的，对所述第二波形数据和幅度校正后的所述第一波形数据进行幅度一致性处理包含：获取激波管激励声波后所述第一波形数据中的平稳信号，根据所述平稳信号计算获得幅度均值；将所述第一波形数据和所述第二波形数据分别除以所述幅度均值，获得一致性处理后的所述第一波形数据和所述第二波形数据。
在上述激波管实验斯通利波幅度提取方法中，优选的，利用斯通利波速度对第一次时间校正后的所述第二波形数据进行第二次时间校正包含：通过将采集所述第二波形数据的探头与岩心顶端之间的距离除以斯通利波速度，对所述第二波形数据进行再次时间校正。
在上述激波管实验斯通利波幅度提取方法中，优选的，通过第二次时间校正后的所述第二波形数据获取斯通利波幅度包含：通过第二次时间校正后的所述第二波形数据中波包的波峰与波谷之间的差值获得斯通利波幅度。
本发明还提供一种激波管实验斯通利波幅度提取装置，所述装置包含：采集模块、校正模块和分析模块；所述采集模块用于获取激波管实验过程中激波管激励声波后水下的第一波形数据和岩心侧面的第二波形数据；所述校正模块用于对所述第一波形数据进行幅度校正使波形基线归零；对所述第二波形数据和幅度校正后的所述第一波形数据进行幅度一致性处理后进行第一次时间校正，使波形起始时间为波形到达岩心顶端时间；利用斯通利波速度对第一次时间校正后的所述第二波形数据进行第二次时间校正，使所述第二波形数据的多道接收波形中斯通利波到达时间一致；所述分析模块用于通过第二次时间校正后的所述第二波形数据获取斯通利波幅度。
在上述激波管实验斯通利波幅度提取装置中，优选的，所述校正模块包含归零单元，所述归零单元用于计算所述第一波形数据中的噪声段的幅度均值，根据所述幅度均值修正波形数据使波形基线归零。
在上述激波管实验斯通利波幅度提取装置中，优选的，所述校正模块包含调整单元，所述调整单元用于获取激波管激励声波后所述第一波形数据中的平稳信号，根据所述平稳信号计算获得幅度均值；将所述第一波形数据和所述第二波形数据分别除以所述幅度均值，获得一致性处理后的所述第一波形数据和所述第二波形数据。
本发明还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述方法。
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述方法的计算机程序。
本发明的有益技术效果在于：提供了一种激波管实验中斯通利波幅度提取方法及装置，拥有完整的激波管实验数据预处理过程，为利用斯通利波定量评价岩石裂缝渗透率奠定了基础。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，并不构成对本发明的限定。在附图中：
图1A为本发明一实施例所提供的激波管实验装置示意图。
图1B为本发明一实施例所提供的激波管实验中斯通利波幅度提取方法的流程示意图。
图2A至图2B为本发明一实施例所提供的波形基线归零前后P1探头接收波形对比图。
图3A至图3B为本发明一实施例所提供的幅度一致性处理前后P1探头接收波形对比图。
图4A至图4B为本发明一实施例所提供的首次时间校正前后P1探头接收波形对比图。
图5A至图5B为本发明一实施例所提供的首次时间校正前后P2探头接收波形对比图。
图6为本发明一实施例所提供的二次时间校正后P2探头接收波形图。
图7为本发明一实施例所提供的二次时间校正后P2探头接收阵列波形图。
图8为本发明一实施例所提供的斯通利波幅度与岩心上移距离关系图。
图9为本发明一实施例所提供的激波管实验中斯通利波幅度提取装置的结构示意图。
图10为本发明一实施例所提供的电子设备的结构示意图。
具体实施方式
以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达成技术效果的实现过程能充分理解并据以实施。需要说明的是，只要不构成冲突，本发明中的各个实施例及各实施例中的各个特征可以相互结合，所形成的技术方案均在本发明的保护范围之内。
另外，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
请参考图1B所示，本发明所提供的一种本发明所提供的激波管实验斯通利波幅度提取方法，具体包含：S101获取激波管实验过程中激波管激励声波后水下的第一波形数据和岩心侧面的第二波形数据；S102对所述第一波形数据进行幅度校正使波形基线归零；S103对所述第二波形数据和幅度校正后的所述第一波形数据进行幅度一致性处理后进行第一次时间校正，使波形起始时间为波形到达岩心顶端时间；S104利用斯通利波速度对第一次时间校正后的所述第二波形数据进行第二次时间校正，使所述第二波形数据的多道接收波形中斯通利波到达时间一致；S105通过第二次时间校正后的所述第二波形数据获取斯通利波幅度。整体上，本实施例是利用激波管实验研究斯通利波衰减规律进而进行裂缝定量评价的关键，实际工作中，为实施上述实施例可通过如下流程实现：
步骤1：搜集激波管实验数据，包括始终置于水中的P1探头所接收波形数据，以及始终置于岩心侧面的P2探头接收波形数据。
步骤2：对P1探头接收波形实施幅度校正，使波形基线归零；
步骤3：对P1和P2探头接收波形实施幅度一致性处理；
步骤4：对P1和P2探头接收波形实施首次时间校正，使波形起始时间为声波到达岩心顶端时间；
步骤5：利用斯通利波速度对P2探头接收波形实施二次时间校正，使得多道接收波形中斯通利波到达时间一致；
步骤6：从P2探头接收波形中提取斯通利波幅度。
在上述步骤1中对岩心实施一次激波管实验后P1和P2探头均会记录1道波形，然后，岩心在激波管中上移0.5cm，再次实施激波管实验。通常岩心移动15次，从而获取P1和P2探头记录的共计30道波形数据；岩心移动的具体次数根据岩心的长度及裂缝所在位置确定。具体的，所述步骤1中搜集激波管实验数据，包括所有测量次数下P1和P2探头记录的波形。搜集激波管实验数据，包括始终置于水中的P1探头所接收波形数据，以及始终置于岩心侧面的P2探头接收波形数据。具体可参考图1A所示，激波管实际上为一密闭腔，一块含裂缝岩心置于底端，被水层覆盖。水层上为常压空气层，再往上为薄膜封隔的高压空气层。激波管工作原理是薄膜被电流瞬间切开，然后高压空气形成冲击波向下传播，传播到水层中变为声波信号，首先被P1探头接收到；然后再往下传播，进入岩心和激波管壁之间的环形空间内，形成导波：斯通利波。斯通利波会被P2探头接收到。一次激波管实验能够获取一组波形数据，包括P1和P2探头接收的共计两道形。一次试验后通过轴承使激波管底座托着岩心上升一段距离，通常为0.5cm，再次开展激波管实验记录一组波形，重复该步骤直至完成整个实验流程。对于含裂缝岩心，实验目的就是通过上移岩心，能够测量到P2探头经过裂缝前后的多组波形数据，进而定量评价裂缝对斯通利波幅度的影响。本领域相关技术人员当可知，上述实施例仅为本发明的一种可用实例，并不对其做任何限定。
在本发明一实施例中，上述对所述第一波形数据进行幅度校正使波形基线归零包含：计算所述第一波形数据中的噪声段的幅度均值，根据所述幅度均值修正波形数据使波形基线归零。亦即对应的步骤2，激波管顶部激励的声波首先在空气中传播一段时间，然后进入水层中，进而被始终置于水中的P1探头接收到。在声波到达之前，P1接收波形中仅含噪声，因此有效声波到达之前的这段噪声幅度可用于使波形基线归零。具体方法为计算噪声段的幅度均值，然后从整道波形中减去该均值，便达到基线归零的目的。噪声段幅度均值计算公式如下
Figure BDA0002566313100000051
式中ANoise代表噪声段幅度均值；T0代表噪声段起始时间，通常取为零时刻；T1代表有效声波首次到达时间；WAVE(t)P1代表P1探头接收波形中不同时间的振幅。
具体请参考图2A至图2B所示，图2A展示了岩心上移3cm后，激波管实验对应的P1探头接收的原始波形，从中可以观察到0.2ms为声波到达时间，在此之前均为随机噪声信号。以公式(1)计算出噪声段幅度均值，然后实施基线归零校正，校正后波形如图2B所示，从中可观察到波形基线基本已归零。
在本发明一实施例中，对所述第二波形数据和幅度校正后的所述第一波形数据进行幅度一致性处理包含：获取激波管激励声波后所述第一波形数据中的平稳信号，根据所述平稳信号计算获得幅度均值；将所述第一波形数据和所述第二波形数据分别除以所述幅度均值，获得一致性处理后的所述第一波形数据和所述第二波形数据。实际工作中，如上述步骤3所述，经所述步骤2波形基线归零校正后，P1探头接收到一段时间的零值信号。激波管激励的声波传播到P1探头后，波形会有一个阶跃，然后接收一段平稳信号，直到声波向下传播到岩心顶端后形成的反射波再传播到P1探头。计算P1探头接收的这段平稳信号的幅度均值，然后将P1和P2探头的接收波形均除以该幅度均值，便达到了幅度一致性处理目的。幅度均值计算公式如下
Figure BDA0002566313100000061
AAvg式中代表幅度均值。将计算起始时间设置为T1+Δt而不是T1的原因是，首次接收到有效声波后，波形中会出现一段持续几个周期的振荡信号。T2代表P1探头接收的反射波到达时间。若对该块岩心实施N次激波管实验，即岩心上移N次，那么最后一次即第N次实验时岩心顶端与P1探头距离最近，此时P1探头接收到的反射波到达时间T2N最早。因此，需要取T2N作为幅度均值AAvg的计算截止时间。具体的，请参考图3A所示，基线归零校正后P1探头接收波形中可以观察到波形第一次阶跃时间为有效声波到达时间，0.2ms；第二次阶跃时间为声波向下传播到岩心顶端后形成的反射波到达时间，0.39ms。两次阶跃时间的中间时间段即为平稳信号，利用公式(2)计算了该段平稳信号的幅度均值，大约为4mV。P1和P2探头接收波形除以该幅度均值便达到幅度一致性处理的目的。图3B展示了经幅度一致性处理后的P1探头接收波形。
在上述实施例中，步骤4可包含从P1探头接收波形中可清楚观察和记录到声波首次到达时间T1，以及声波传播到岩心顶端反射回来的到达时间T2。时间校正便是将声波到达岩心顶端的时间校正为零时刻，校正公式如下
t＝t-(T1+T2)/2 (3)
需要注意的是，该校正公式同时适用于激波管实验的P1和P2探头接收波形，并且每次激波管实验的校正量均需重新计算。具体的，步骤3中已求得第一次阶跃时间T1(0.2ms)和第二次阶跃时间T2(0.39ms)，利用公式(3)可求得时间校正量0.295ms。图4A和图4B分别展示了时间校正前后的P1探头接收波形，从中可观擦到，时间校正便是将波形向前移动了0.295ms。图5A和图5B分别展示了时间校正前后的P2探头接收波形。
在本发明一实施例中，利用斯通利波速度对第一次时间校正后的所述第二波形数据进行第二次时间校正包含：通过将采集所述第二波形数据的探头与岩心顶端之间的距离除以斯通利波速度，对所述第二波形数据进行再次时间校正。实际工作中，如步骤5所示，在经过所述步骤4时间校正后的P2探头接收波形中，斯通利波到达时间便为P2探头与岩心顶端距离除以斯通利波传播速度，具体公式如下：
TSTi＝Si/VST (4)
式中TSTi代表第i次激波管实验时斯通利波到达时间；Si代表第次激波管实验时P2探头与岩心顶端距离；VST代表斯通利波速度，对于碳酸盐岩、砂岩等岩石，Vst略小于流体声速，可取值为1450m/s。
若将多次激波管实验中斯通利波到达时间校正为一致，则每次P2接收波形对应时间需减去斯通利波到达时间TSTi即可。具体的，请结合前述实施例参考图6所示，当岩心上移3cm后，岩心顶端与P2探头距离经测量可知为0.286m，根据公式(4)计算出此时斯通利波到达时间TSTi为0.197ms，以此时间对该次激波管实验对应的P2探头接收波形实施二次时间校正，校正后波形如图6所示，从中可观察到斯通利波到达时间已校正到零时刻。图7进一步展示了对该岩心实施15次激波管实验对应的经二次时间校正后的P2探头接收波形，可以看到所有接收波形中斯通利波均已校正到位。
在本发明一实施例中，通过第二次时间校正后的所述第二波形数据获取斯通利波幅度包含：通过第二次时间校正后的所述第二波形数据中波包的波峰与波谷之间的差值获得斯通利波幅度。如步骤6所示，实际工作中，P2探头接收波形中斯通利波信号通常显示为一个波包，存在一个波峰和一个波谷，计算两者的差值便为斯通利波幅度。具体请参考图6所示，首先提取斯通利波峰值和谷值，然后计算两者差值即为斯通利波幅度。依据此方法计算出多组波形对应的斯通利波幅度值，它们与岩心上移距离的关系如图8所示，能够看到经过裂缝前斯通利波幅度较高，在经过裂缝之后斯通利波幅度有明显降低，这说明裂缝对斯通利波幅度确实有显著的衰减影响。
请参考图9所示，本发明还提供一种激波管实验斯通利波幅度提取装置，所述装置包含：采集模块、校正模块和分析模块；所述采集模块用于获取激波管实验过程中激波管激励声波后水下的第一波形数据和岩心侧面的第二波形数据；所述校正模块用于对所述第一波形数据进行幅度校正使波形基线归零；对所述第二波形数据和幅度校正后的所述第一波形数据进行幅度一致性处理后进行第一次时间校正，使波形起始时间为波形到达岩心顶端时间；利用斯通利波速度对第一次时间校正后的所述第二波形数据进行第二次时间校正，使所述第二波形数据的多道接收波形中斯通利波到达时间一致；所述分析模块用于通过第二次时间校正后的所述第二波形数据获取斯通利波幅度。
在上述实施例中，所述校正模块可包含归零单元和调整单元，所述归零单元用于计算所述第一波形数据中的噪声段的幅度均值，根据所述幅度均值修正波形数据使波形基线归零。所述调整单元用于获取激波管激励声波后所述第一波形数据中的平稳信号，根据所述平稳信号计算获得幅度均值；将所述第一波形数据和所述第二波形数据分别除以所述幅度均值，获得一致性处理后的所述第一波形数据和所述第二波形数据。
本发明的有益技术效果在于：提供了一种激波管实验中斯通利波幅度提取方法及装置，拥有完整的激波管实验数据预处理过程，为利用斯通利波定量评价岩石裂缝渗透率奠定了基础。
本发明还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述方法。
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述方法的计算机程序。
如图10所示，该电子设备600还可以包括：通信模块110、输入单元120、音频处理单元130、显示器160、电源170。值得注意的是，电子设备600也并不是必须要包括图10中所示的所有部件；此外，电子设备600还可以包括图10中没有示出的部件，可以参考现有技术。
如图10所示，中央处理器100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器100接收输入并控制电子设备600的各个部件的操作。
其中，存储器140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器100可执行该存储器140存储的该程序，以实现信息存储或处理等。
输入单元120向中央处理器100提供输入。该输入单元120例如为按键或触摸输入装置。电源170用于向电子设备600提供电力。显示器160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。
该存储器140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器140还可以是某种其它类型的装置。存储器140包括缓冲存储器141(有时被称为缓冲器)。存储器140可以包括应用/功能存储部142，该应用/功能存储部142用于存储应用程序和功能程序或用于通过中央处理器100执行电子设备600的操作的流程。
存储器140还可以包括数据存储部143，该数据存储部143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器140的驱动程序存储部144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块110即为经由天线111发送和接收信号的发送机/接收机110。通信模块(发送机/接收机)110耦合到中央处理器100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。
基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)110还经由音频处理器130耦合到扬声器131和麦克风132，以经由扬声器131提供音频输出，并接收来自麦克风132的音频输入，从而实现通常的电信功能。音频处理器130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器130还耦合到中央处理器100，从而使得可以通过麦克风132能够在本机上录音，且使得可以通过扬声器131来播放本机上存储的声音。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种激波管实验斯通利波幅度提取方法，其特征在于，所述方法包含：
获取激波管实验过程中激波管激励声波后水下的第一波形数据和岩心侧面的第二波形数据；
对所述第一波形数据进行幅度校正使波形基线归零；
对所述第二波形数据和幅度校正后的所述第一波形数据进行幅度一致性处理后进行第一次时间校正，使波形起始时间为波形到达岩心顶端时间；
利用斯通利波速度对第一次时间校正后的所述第二波形数据进行第二次时间校正，使所述第二波形数据的多道接收波形中斯通利波到达时间一致；
通过第二次时间校正后的所述第二波形数据获取斯通利波幅度。
2.根据权利要求1所述的激波管实验斯通利波幅度提取方法，其特征在于，对所述第一波形数据进行幅度校正使波形基线归零包含：计算所述第一波形数据中的噪声段的幅度均值，根据所述幅度均值修正波形数据使波形基线归零。
3.根据权利要求1所述的激波管实验斯通利波幅度提取方法，其特征在于，对所述第二波形数据和幅度校正后的所述第一波形数据进行幅度一致性处理包含：
获取激波管激励声波后所述第一波形数据中的平稳信号，根据所述平稳信号计算获得幅度均值；
将所述第一波形数据和所述第二波形数据分别除以所述幅度均值，获得一致性处理后的所述第一波形数据和所述第二波形数据。
4.根据权利要求1所述的激波管实验斯通利波幅度提取方法，其特征在于，利用斯通利波速度对第一次时间校正后的所述第二波形数据进行第二次时间校正包含：
通过将采集所述第二波形数据的探头与岩心顶端之间的距离除以斯通利波速度，对所述第二波形数据进行再次时间校正。
5.根据权利要求1所述的激波管实验斯通利波幅度提取方法，其特征在于，通过第二次时间校正后的所述第二波形数据获取斯通利波幅度包含：通过第二次时间校正后的所述第二波形数据中波包的波峰与波谷之间的差值获得斯通利波幅度。
6.一种激波管实验斯通利波幅度提取装置，其特征在于，所述装置包含采集模块、校正模块和分析模块；
所述采集模块用于获取激波管实验过程中激波管激励声波后水下的第一波形数据和岩心侧面的第二波形数据；
所述校正模块用于对所述第一波形数据进行幅度校正使波形基线归零；对所述第二波形数据和幅度校正后的所述第一波形数据进行幅度一致性处理后进行第一次时间校正，使波形起始时间为波形到达岩心顶端时间；利用斯通利波速度对第一次时间校正后的所述第二波形数据进行第二次时间校正，使所述第二波形数据的多道接收波形中斯通利波到达时间一致；
所述分析模块用于通过第二次时间校正后的所述第二波形数据获取斯通利波幅度。
7.根据权利要求6所述的激波管实验斯通利波幅度提取装置，其特征在于，所述校正模块包含归零单元，所述归零单元用于计算所述第一波形数据中的噪声段的幅度均值，根据所述幅度均值修正波形数据使波形基线归零。
8.根据权利要求6所述的激波管实验斯通利波幅度提取装置，其特征在于，所述校正模块包含调整单元，所述调整单元用于获取激波管激励声波后所述第一波形数据中的平稳信号，根据所述平稳信号计算获得幅度均值；将所述第一波形数据和所述第二波形数据分别除以所述幅度均值，获得一致性处理后的所述第一波形数据和所述第二波形数据。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5任一所述方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至5任一所述方法的计算机程序。
CN202010625075.4A 2020-07-02 2020-07-02 激波管实验斯通利波幅度提取方法及装置 Active CN113884580B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010625075.4A CN113884580B (zh) 2020-07-02 2020-07-02 激波管实验斯通利波幅度提取方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010625075.4A CN113884580B (zh) 2020-07-02 2020-07-02 激波管实验斯通利波幅度提取方法及装置

Publications (2)

Publication Number Publication Date
CN113884580A true CN113884580A (zh) 2022-01-04
CN113884580B CN113884580B (zh) 2023-08-22

Family

ID=79012161

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010625075.4A Active CN113884580B (zh) 2020-07-02 2020-07-02 激波管实验斯通利波幅度提取方法及装置

Country Status (1)

Country Link
CN (1) CN113884580B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102222359A (zh) * 2011-05-24 2011-10-19 中国石油天然气股份有限公司 一种岩心三维孔隙结构重构方法
CN104345346A (zh) * 2014-10-30 2015-02-11 中国石油天然气股份有限公司 一种获取裂缝宽度的方法
CN108375785A (zh) * 2018-01-12 2018-08-07 中国石油天然气股份有限公司 裂缝带位置校正方法及装置
CN110173259A (zh) * 2019-05-14 2019-08-27 中石化石油工程技术服务有限公司 一种斯通利波能量补偿校正的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102222359A (zh) * 2011-05-24 2011-10-19 中国石油天然气股份有限公司 一种岩心三维孔隙结构重构方法
CN104345346A (zh) * 2014-10-30 2015-02-11 中国石油天然气股份有限公司 一种获取裂缝宽度的方法
CN108375785A (zh) * 2018-01-12 2018-08-07 中国石油天然气股份有限公司 裂缝带位置校正方法及装置
CN110173259A (zh) * 2019-05-14 2019-08-27 中石化石油工程技术服务有限公司 一种斯通利波能量补偿校正的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
NING LI ET AL.,: ""Shock-induced Stoneley waves in carbonate rock samples"", 《GEOPHYSICS》, vol. 84, no. 5, pages 209 - 216 *
杨博 等: ""基于斯通利波能量衰减的裂缝渗透性评价方法研究"", 《地球物理学进展》, vol. 34, no. 3, pages 1127 - 1131 *

Also Published As

Publication number Publication date
CN113884580B (zh) 2023-08-22

Similar Documents

Publication Publication Date Title
CN108562938B (zh) 一种消除频散效应的方法、装置及系统
CN104345346A (zh) 一种获取裂缝宽度的方法
CN112593922B (zh) 一种阵列声波测井评价固井二界面胶结质量的方法及装置
CN105116448B (zh) 一种转换波方位各向异性校正方法及装置
CN105510880A (zh) 一种基于双差法的微地震震源定位方法
CN110426739B (zh) 一种地质勘探检测方法、装置和存储介质
Cui et al. Assessing the effects of insufficient rebar and missing grout in grouted rock bolts using guided ultrasonic waves
CN202611692U (zh) 在声波测井中消除直达波干扰的系统及声波测井仪
CN103362502B (zh) 在声波测井中消除直达波干扰的方法、系统及声波测井仪
CN103114844A (zh) 水平井声波水泥胶结测井中仪器偏心校正方法
Lei et al. An effective method for laboratory acoustic emission detection and location using template matching
RU2009139867A (ru) Автоматизированная оценка медленности бурового раствора
CN111381292B (zh) 一种预测砂岩含烃储层的测井解释方法与装置
CN104237383B (zh) 一种岩石超声波速度及品质因子同时测量方法及系统
CN107238859B (zh) 计算虚拟井合成地震记录的方法及系统
CN113884580A (zh) 激波管实验斯通利波幅度提取方法及装置
CN108490490B (zh) 裂缝带表征方法、装置、电子设备及计算机存储介质
CN110596747A (zh) 一种可控震源激发节点接收的表层调查方法
CN113672840B (zh) 斯通利波衰减提取方法及装置
CN113077074B (zh) 基于储层预测因子的储层预测方法及装置
CN114814951A (zh) 一种获取快慢横波慢度的方法、装置、设备及存储介质
CN116413787A (zh) 利用das井中地震数据反演地层波阻抗的方法及装置
CN112379407B (zh) 可控震源激发组合台数与覆盖次数的匹配方法及装置
CN112363225B (zh) 横波微测井解释方法及装置
CN112253100B (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