[go: up one dir, main page]

CN114624772A - Well-seismic calibration method, device, equipment and storage medium for sonic-free curve well - Google Patents

Well-seismic calibration method, device, equipment and storage medium for sonic-free curve well Download PDF

Info

Publication number
CN114624772A
CN114624772A CN202011439612.2A CN202011439612A CN114624772A CN 114624772 A CN114624772 A CN 114624772A CN 202011439612 A CN202011439612 A CN 202011439612A CN 114624772 A CN114624772 A CN 114624772A
Authority
CN
China
Prior art keywords
well
curve
time difference
silent
acoustic
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.)
Pending
Application number
CN202011439612.2A
Other languages
Chinese (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 CN202011439612.2A priority Critical patent/CN114624772A/en
Publication of CN114624772A publication Critical patent/CN114624772A/en
Pending legal-status Critical Current

Links

Images

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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V13/00Manufacturing, calibrating, cleaning, or repairing instruments or devices covered by groups G01V1/00 – G01V11/00

Landscapes

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

Abstract

The application provides a well seismic calibration method, a well seismic calibration device, well seismic calibration equipment and a storage medium for a silent wave curve well. In the method, a first acoustic time difference curve of an adjacent well of a silent wave curve well to be calibrated and a reference curve of the silent wave curve well are obtained. And establishing a response relation equation between the reference curve and the first sound wave time difference curve, and performing multivariate linear fitting processing to obtain a second sound wave time difference curve of the silent wave curve well. And finally, carrying out well seismic calibration on the silent wave curve well according to the second sound wave time difference curve. According to the technical scheme, the well seismic calibration of the silent wave curve well is realized by acquiring the first acoustic time difference curve of the reference curve and the adjacent well and performing multivariate linear fitting processing, and the well seismic calibration efficiency is improved.

Description

无声波曲线井的井震标定方法、装置、设备及存储介质Well-seismic calibration method, device, equipment and storage medium for non-acoustic curve wells

技术领域technical field

本申请涉及测井技术领域,尤其涉及一种无声波曲线井的井震标定方法、装置、设备及存储介质。The present application relates to the technical field of well logging, and in particular, to a well seismic calibration method, device, equipment and storage medium for a non-acoustic curve well.

背景技术Background technique

声波测井数据记录了声波在地层中的传播时间,反映了地层的孔隙特性及波阻抗特性,这与地震数据在反映地层真实属性方面具有很好的一致性。地震资料在横向上具有范围广,连续性好的特点,能够对地层进行横向追踪和预测,具有较高的纵向分辨率,因此将声波测井资料与井旁地震数据联合起来应用到地震勘探中成为了必要且常用的技术手段。Acoustic logging data records the propagation time of acoustic waves in the formation, and reflects the pore characteristics and wave impedance characteristics of the formation, which is in good agreement with seismic data in reflecting the real properties of the formation. The seismic data has the characteristics of wide range and good continuity in the lateral direction, which can track and predict the formation laterally, and has a high vertical resolution. Therefore, the acoustic logging data and the seismic data next to the well are combined and applied to seismic exploration. become a necessary and commonly used technical means.

在井震标定的过程中,经常会遇到缺失声波曲线的情况,通常此时在井震标定及储层反演中就会放弃该井。随着测井技术的发展,对此种情况常用的解决办法是通过岩石物理建模的方法。In the process of well-seismic calibration, it is often encountered that the acoustic curve is missing, and the well is usually abandoned in the well-seismic calibration and reservoir inversion at this time. With the development of logging technology, the common solution to this situation is through the method of petrophysical modeling.

然而,在现有技术中,这种建模方法需要大量的岩心分析、压力-体积-温度分析(Pressure-Volume-Temperature,PVT)、体积组分、含水饱和度、孔隙度等基础资料,且经济成本和时间成本都较高,如果缺少这些资料,该方法将不适用解决无声波时差曲线井的井震标定问题。However, in the prior art, this modeling method requires a lot of basic data such as core analysis, pressure-volume-temperature analysis (PVT), volume composition, water saturation, porosity, etc., and Both the economic cost and the time cost are high. If these data are lacking, this method will not be suitable for solving the well-seismic calibration problem of the well with the silent time difference curve.

发明内容SUMMARY OF THE INVENTION

本申请提供一种无声波曲线井的井震标定方法、装置、设备及存储介质,用以解决既无声波时差曲线又缺少必要基础资料的井,无法进行井震标定的问题。The present application provides a well-seismic calibration method, device, equipment and storage medium for a well without an acoustic wave curve, which is used to solve the problem that the well-seismic calibration cannot be performed for a well without an acoustic time difference curve and lacking necessary basic data.

第一方面,本申请实施例提供一种无声波曲线井的井震标定方法,包括:In a first aspect, an embodiment of the present application provides a well-seismic calibration method for a non-acoustic curve well, including:

获取待标定的无声波曲线井的邻井的第一声波时差曲线,所述邻井与所述无声波曲线井之间的距离小于预设值,且特征相同;Acquiring a first sonic time difference curve of an adjacent well of a non-acoustic curve well to be calibrated, where the distance between the adjacent well and the non-acoustic curve well is less than a preset value and has the same characteristics;

获取所述无声波曲线井的基准曲线;obtaining the reference curve of the silent curve well;

对所述基准曲线和所述第一声波时差曲线之间建立响应关系方程,并进行多元线性拟合处理,得到所述无声波曲线井的第二声波时差曲线;establishing a response relationship equation between the reference curve and the first sonic time difference curve, and performing multivariate linear fitting processing to obtain a second sonic time difference curve of the non-acoustic curve well;

根据所述第二声波时差曲线,对所述无声波曲线井进行井震标定。According to the second sonic time difference curve, the well-seismic calibration is performed on the well without the sonic curve.

在第一方面的一种可能设计中,所述获取待标定的无声波曲线井的邻井的第一声波时差曲线,包括:In a possible design of the first aspect, the acquiring the first acoustic transit time curve of the adjacent well of the non-acoustic curve well to be calibrated includes:

根据所述无声波曲线井的特征,选取与所述无声波曲线井的距离小于所述预设距离且特征一致的测井作为所述邻井;According to the characteristics of the silent curve well, a logging well whose distance from the silent curve well is less than the preset distance and has the same characteristics is selected as the offset well;

从所述邻井中满足预设井眼条件的井段中选择样本层段,所述样本层段的范围为2米至10米;selecting a sample interval from the well intervals in the offset well that meet the preset wellbore conditions, and the sample interval ranges from 2 meters to 10 meters;

根据从所述样本层段中获取的多个样本点,确定所述第一声波时差曲线。The first acoustic transit time curve is determined according to a plurality of sample points obtained from the sample slice.

在该种可能的设计中,所述方法还包括:In this possible design, the method further includes:

从所述样本层段中获取多个样本点,所述多个样本点的个数范围为20个至80个,且所述多个样本点的岩性和流性的类型与所述无声波曲线井的部分井段数据相近。A plurality of sample points are obtained from the sample interval, the number of the plurality of sample points ranges from 20 to 80, and the types of lithology and fluidity of the plurality of sample points are the same as those of the silent wave Some sections of the curve well have similar data.

可选的,所述无声波曲线井和所述邻井的沉积环境和压实程度相同。Optionally, the deposition environment and the degree of compaction of the non-acoustic curve well and the adjacent well are the same.

可选的,所述基准曲线为自然伽马曲线或者/和深电阻率曲线。Optionally, the reference curve is a natural gamma curve or/and a deep resistivity curve.

可选的,所述响应关系方程为:DT=a*GR+b*DEN+c*lg(RD)+N;Optionally, the response relationship equation is: DT=a*GR+b*DEN+c*lg(RD)+N;

其中,DT表示声波时差,GR表示自然伽马,DEN表示密度,RD表示深电阻率,a、b、c和N分别为拟合常数。Among them, DT stands for acoustic transit time, GR stands for natural gamma, DEN stands for density, RD stands for deep resistivity, and a, b, c, and N are fitting constants, respectively.

在第一方面的另一种可能设计中,所述方法还包括:In another possible design of the first aspect, the method further includes:

获取采用所述第二声波时差曲线进行的井震标定结果与实测得到的地震波形之间的差异程度;Obtaining the degree of difference between the well-seismic calibration result performed by using the second acoustic wave time difference curve and the measured seismic waveform;

若所述差异程度大于预设阈值,则重新从所述邻井中选择样本点,获取新的第一声波时差曲线;If the degree of difference is greater than the preset threshold, reselect sample points from the adjacent wells to obtain a new first acoustic wave time difference curve;

根据所述新的第一声波时差曲线和所述基准曲线进行拟合处理,得到新的第二声波时差曲线。The fitting process is performed according to the new first acoustic wave time difference curve and the reference curve to obtain a new second acoustic wave time difference curve.

第二方面,本申请提供一种无声波曲线井的井震标定装置,包括:获取模块和处理模块;In a second aspect, the present application provides a well-seismic calibration device for a non-acoustic curve well, comprising: an acquisition module and a processing module;

所述处理模块,用于获取待标定的无声波曲线井的邻井的第一声波时差曲线,所述邻井与所述无声波曲线井之间的距离小于预设值,且特征相同;The processing module is used to obtain the first sonic time difference curve of the adjacent well of the non-acoustic curve well to be calibrated, the distance between the adjacent well and the non-acoustic curve well is less than a preset value, and the characteristics are the same;

所述获取模块,用于获取所述无声波曲线井的基准曲线;the acquisition module, configured to acquire the reference curve of the soundless curve well;

所述处理模块,还用于对所述基准曲线和所述第一声波时差曲线之间建立响应关系方程,并进行多元线性拟合处理,得到所述无声波曲线井的第二声波时差曲线;The processing module is further configured to establish a response relationship equation between the reference curve and the first sonic time difference curve, and perform multivariate linear fitting processing to obtain the second sonic time difference curve of the non-acoustic curve well ;

所述处理模块,还用于根据所述第二声波时差曲线,对所述无声波曲线井进行井震标定。The processing module is further configured to perform well-seismic calibration on the well without the acoustic wave curve according to the second acoustic wave time difference curve.

在第二方面的一种可能设计中,所述处理模块,用于获取待标定的无声波曲线井的邻井的第一声波时差曲线,包括:In a possible design of the second aspect, the processing module, configured to obtain the first acoustic transit time curve of the adjacent well of the non-acoustic curve well to be calibrated, includes:

根据所述无声波曲线井的特征,选取与所述无声波曲线井的距离小于所述预设距离且特征一致的测井作为所述邻井;According to the characteristics of the silent curve well, a logging well whose distance from the silent curve well is less than the preset distance and has the same characteristics is selected as the offset well;

从所述邻井中满足预设井眼条件的井段中选择样本层段,所述样本层段的范围为2米至10米;selecting a sample interval from the well intervals in the offset well that meet the preset wellbore conditions, and the sample interval ranges from 2 meters to 10 meters;

根据从所述样本层段中获取的多个样本点,确定所述第一声波时差曲线。The first acoustic transit time curve is determined according to a plurality of sample points obtained from the sample slice.

在该种可能的设计中,所述处理模块还包括从所述样本层段中获取多个样本点,所述多个样本点的个数范围为20个至80个,且所述多个样本点的岩性和流性的类型与所述无声波曲线井的部分井段数据相近。In this possible design, the processing module further includes acquiring a plurality of sample points from the sample layer, the number of the plurality of sample points ranges from 20 to 80, and the plurality of samples The types of lithology and fluidity of the points are similar to the data of some sections of the silent curve well.

可选的,所述无声波曲线井和所述邻井的沉积环境和压实程度相同。Optionally, the deposition environment and the degree of compaction of the non-acoustic curve well and the adjacent well are the same.

可选的,所述基准曲线为自然伽马曲线或者/和深电阻率曲线。Optionally, the reference curve is a natural gamma curve or/and a deep resistivity curve.

可选的,所述响应关系方程为:DT=a*GR+b*DEN+c*lg(RD)+N;Optionally, the response relationship equation is: DT=a*GR+b*DEN+c*lg(RD)+N;

其中,DT表示声波时差,GR表示自然伽马,DEN表示密度,RD表示深电阻率,a、b、c和N分别为拟合常数。Among them, DT stands for acoustic transit time, GR stands for natural gamma, DEN stands for density, RD stands for deep resistivity, and a, b, c, and N are fitting constants, respectively.

在第二方面的另一种可能设计中,所述处理模块,还用于:In another possible design of the second aspect, the processing module is further configured to:

获取采用所述第二声波时差曲线进行的井震标定结果与实测得到的地震波形之间的差异程度;Obtaining the degree of difference between the well-seismic calibration result performed by using the second acoustic wave time difference curve and the measured seismic waveform;

若所述差异程度大于预设阈值,则重新从所述邻井中选择样本点,获取新的第一声波时差曲线;If the degree of difference is greater than the preset threshold, reselect sample points from the adjacent wells to obtain a new first acoustic wave time difference curve;

根据所述新的第一声波时差曲线和所述基准曲线进行拟合处理,得到新的第二声波时差曲线。The fitting process is performed according to the new first acoustic wave time difference curve and the reference curve to obtain a new second acoustic wave time difference curve.

第三方面,本申请实施例提供一种计算机设备,包括:处理器、存储器和收发器;In a third aspect, an embodiment of the present application provides a computer device, including: a processor, a memory, and a transceiver;

所述存储器用于存储可在处理器上运行的计算机程序指令;the memory for storing computer program instructions executable on the processor;

所述处理器执行所述计算机程序指令时实现如上述第一方面以及各可能设计提供的方法。The processor, when executing the computer program instructions, implements the methods provided by the first aspect and possible designs described above.

第四方面,本申请实施例可提供一种计算机可读存储介质,所述计算机可读存储介质中存储有计算机执行指令,当所述计算机执行指令被处理器执行时用于实现第一方面以及各可能设计提供的方法。In a fourth aspect, embodiments of the present application may provide a computer-readable storage medium, where computer-executable instructions are stored in the computer-readable storage medium, and when the computer-executable instructions are executed by a processor, are used to implement the first aspect and methods provided by each possible design.

第五方面,本申请实施例提供一种程序,当该程序被处理器执行时,用于执行如第一方面所述的方法。In a fifth aspect, an embodiment of the present application provides a program for executing the method according to the first aspect when the program is executed by a processor.

第六方面,本申请实施例提供一种计算机程序产品,包括程序指令,程序指令用于实现如第一方面所述的方法。In a sixth aspect, embodiments of the present application provide a computer program product, including program instructions, where the program instructions are used to implement the method described in the first aspect.

本申请实施例提供的无声波曲线井的井震标定方法、装置、设备及存储介质,通过获取待标定的无声波曲线井的邻井的第一声波时差曲线,以及无声波曲线井的基准曲线。对基准曲线和第一声波时差曲线之间建立响应关系方程,并进行多元线性拟合处理,得到无声波曲线井的第二声波时差曲线。最后根据第二声波时差曲线,对无声波曲线井进行井震标定。该技术方案中,通过获取基准曲线与邻井的第一声波时差曲线,并进行多元线性拟合处理,实现了对无声波曲线井的井震标定,并提高了井震标定的效率。The well-seismic calibration method, device, device, and storage medium for a non-acoustic curve well provided by the embodiments of the present application are obtained by acquiring the first acoustic time difference curve of the adjacent well of the non-acoustic curve well to be calibrated, and the benchmark of the non-acoustic curve well curve. A response relationship equation is established between the reference curve and the first acoustic time difference curve, and multivariate linear fitting is performed to obtain the second acoustic time difference curve of the well without acoustic wave curve. Finally, according to the second sonic time difference curve, the well seismic calibration is carried out for the well without sonic curve. In this technical scheme, by obtaining the reference curve and the first acoustic time difference curve of the adjacent well, and performing multivariate linear fitting processing, the well seismic calibration of the well without the acoustic wave curve is realized, and the efficiency of the well seismic calibration is improved.

附图说明Description of drawings

图1为本申请实施例提供的无声波曲线井的井震标定方法实施例一的流程图;1 is a flow chart of Embodiment 1 of a well-seismic calibration method for a non-acoustic curve well provided by an embodiment of the present application;

图2为本申请实施例提供的无声波曲线井的井震标定方法实施例二的流程图;2 is a flow chart of Embodiment 2 of the well-seismic calibration method for a non-acoustic curve well provided by an embodiment of the present application;

图3为本申请实施例提供的无声波曲线井的井震标定方法实施例三的流程图;3 is a flowchart of Embodiment 3 of the well-seismic calibration method for a non-acoustic curve well provided by an embodiment of the present application;

图4为本申请实施例提供的无声波曲线井的井震标定方法的总流程图;Fig. 4 is the general flow chart of the well-seismic calibration method of the non-acoustic curve well provided by the embodiment of the present application;

图5为本申请实施例提供的无声波曲线井的井震标定装置的结构示意图;5 is a schematic structural diagram of a well seismic calibration device for a non-acoustic curve well provided in an embodiment of the present application;

图6为本申请实施例提供的计算机设备的结构示意图。FIG. 6 is a schematic structural diagram of a computer device provided by an embodiment of the present application.

具体实施方式Detailed ways

为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。In order to make the purposes, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments It is a part of the embodiments of the present application, but not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those of ordinary skill in the art without creative work fall within the protection scope of the present application.

在介绍本申请的实施例之前,首先对本申请的背景技术进行解释说明:Before introducing the embodiments of the present application, the background technology of the present application will be explained first:

井震标定是连接测井与地震信息的桥梁,精细的地震标定是进行层位解释和油藏描述的基础和前提。而影响井震标结果的主要原因有:测井资料、地震资料、子波以及其他因素。Well-seismic calibration is a bridge connecting logging and seismic information, and fine seismic calibration is the basis and premise of horizon interpretation and reservoir description. The main factors that affect the results of well seismic markers are: logging data, seismic data, wavelets and other factors.

井震标定结果不准确可能会导致储层构造深度误差、储层预测不准,以及油藏认识偏移。在井震标定时,常用的方法是通过岩石物理建模方法解决无声波时差曲线井的井震标定问题,具体的,利用岩心分析、压力-体积-温度分析(Pressure-Volume-Temperature,PVT)、体积组分、含水饱和度、孔隙度等基础资料来进行岩石物理建模。但是这种方法经济成本和时间成本都较高,如果缺少这些资料,该方法将不适用解决无声波时差曲线井的井震标定问题。另外,这种方法对于技术人员的地震解释经验具有依赖性,导致分析结果具有较强的主观性。Inaccurate well-seismic calibration results may lead to reservoir structural depth errors, inaccurate reservoir predictions, and shifts in reservoir knowledge. In the well-seismic calibration, the commonly used method is to solve the well-seismic calibration problem of the non-acoustic transit time curve well through the rock physics modeling method. Specifically, the core analysis, pressure-volume-temperature analysis (Pressure-Volume-Temperature, PVT) , volume composition, water saturation, porosity and other basic data to carry out petrophysical modeling. However, the economic cost and time cost of this method are high. If these data are lacking, this method will not be suitable for solving the well-seismic calibration problem of the wells with non-acoustic transit time curve. In addition, this method is dependent on the seismic interpretation experience of technicians, resulting in strong subjectivity of analysis results.

具体的,当缺失声波曲线时,以塔里木油田为例,第一、塔里木油田井深通常在5000-7000米,取1筒心+PVT+岩心分析,费用约100万左右;第二、岩石物理建模由于参数众多、过程冗长,从测井评价得到Vclay、

Figure BDA0002829675060000051
Sw等参数,建立体积模型,根据工区储层、流体类型选择合适的岩石物理模型、测试调整各类参数,到最终正演出缺失的声波曲线,整个过程基本在2-3天以上;第三、传统采用岩石物理建模的方法,对于没有录取岩心、PVT等基础资料的单井,就无法进行井震标定,并且在储层反演中,也只能放弃该井。Specifically, when the sound wave curve is missing, taking Tarim Oilfield as an example, first, the well depth of Tarim Oilfield is usually 5000-7000 meters, taking 1 tube core + PVT + core analysis, the cost is about 1 million; second, petrophysical modeling Due to the numerous parameters and the lengthy process, the Vclay,
Figure BDA0002829675060000051
Sw and other parameters, establish a volume model, select an appropriate petrophysical model according to the reservoir and fluid type in the work area, test and adjust various parameters, and finally perform the missing sound wave curve, the whole process basically takes more than 2-3 days; third, The traditional petrophysical modeling method cannot perform well-seismic calibration for a single well without basic data such as core and PVT, and in the reservoir inversion, the well can only be abandoned.

针对上述技术问题,本申请的发明构思如下:由于缺少必要的基础资料,可以在需要进行井震标定的无声波曲线井的附近寻找具有相同特征的临井,然后对临井进行特征分析,结合需要进行井震标定的无声波曲线井的特征,建立一定的关系,便可以对该无声波曲线井做出准确地井震标定。In view of the above-mentioned technical problems, the inventive concept of the present application is as follows: due to the lack of necessary basic data, an adjacent well with the same characteristics can be found in the vicinity of the non-acoustic curve well that needs to perform well seismic calibration, and then the characteristics of the adjacent well can be analyzed, combined with According to the characteristics of the non-acoustic curve wells that need to be calibrated by well-seismic, and establishing a certain relationship, the accurate well-seismic calibration of the non-acoustic curve wells can be made.

下面,通过具体实施例对本申请的技术方案进行详细说明。Hereinafter, the technical solutions of the present application will be described in detail through specific embodiments.

需要说明的是,下面这几个具体的实施例可以相互结合,对于相同或相似的概念或过程可能在某些实施例中不再赘述。It should be noted that the following specific embodiments may be combined with each other, and the same or similar concepts or processes may not be repeated in some embodiments.

图1为本申请实施例提供的无声波曲线井的井震标定方法实施例一的流程图。如图1所示,该井震标定方法可以包括如下步骤:FIG. 1 is a flowchart of Embodiment 1 of a well-seismic calibration method for a non-acoustic curve well provided by an embodiment of the present application. As shown in Figure 1, the well seismic calibration method may include the following steps:

步骤11、获取待标定的无声波曲线井的邻井的第一声波时差曲线。Step 11: Acquire the first acoustic time difference curve of the adjacent well of the non-acoustic curve well to be calibrated.

在本步骤中,当需要对待标定的无声波曲线井进行井震标定时(即未进行声波时差测井的井),需要在找到相邻的其他井(即邻井),对邻井进行处理,获取邻井的声波时差曲线(即第一声波时差曲线)。In this step, when it is necessary to perform well-seismic calibration for the non-acoustic curve well to be calibrated (that is, the well without sonic time difference logging), it is necessary to find other adjacent wells (that is, adjacent wells), and process the adjacent wells , to obtain the acoustic transit time curve of the adjacent well (ie, the first acoustic transit time curve).

可选的,邻井与无声波曲线井之间的距离小于预设值,且特征相同。Optionally, the distance between the offset well and the silent curve well is less than a preset value and has the same characteristics.

具体的,为了使得邻井与无声波曲线井之间的地质条件等因素尽量接近,以减小井震标定的误差,需要将邻井到无声波曲线井的距离小于预设值,该预设值可以是人为规定的距离值。Specifically, in order to make the geological conditions and other factors between the offset well and the silent curve well as close as possible to reduce the error of well seismic calibration, the distance from the offset well to the silent curve well needs to be smaller than the preset value. The value can be an artificially specified distance value.

可选的,邻井与无声波曲线井的特征相同可以是:沉积环境和压实程度相同。Optionally, the characteristics of the offset well and the silent curve well may be the same: the deposition environment and the degree of compaction are the same.

步骤12、获取无声波曲线井的基准曲线。Step 12, obtaining the reference curve of the silent curve well.

其中,基准曲线可以是自然伽马曲线或者/和深电阻率曲线。主要是因为自然伽马曲线和深电阻率曲线是石油探测过程中受井眼条件影响较小的两个曲线。The reference curve may be a natural gamma curve or/and a deep resistivity curve. Mainly because the natural gamma curve and the deep resistivity curve are the two curves that are less affected by wellbore conditions in the process of oil exploration.

可选的,自然伽马曲线是由自然伽马测井技术获取的基准曲线,基于地球物理探测手段,了解被测的无声波曲线井中不同井段的导电、传热、声学、放射性等物理性质,用以确定岩层的岩性、物性、含油气性、含盐卤矿化度与地球物理特性的关系。Optionally, the natural gamma curve is a reference curve obtained by the natural gamma logging technology. Based on the geophysical detection method, the physical properties such as electrical conductivity, heat transfer, acoustics, and radioactivity of different well sections in the well measured without sound wave curve can be understood. , to determine the relationship between the lithology, physical properties, oil and gas properties, salt-bearing brine salinity and geophysical properties of rock formations.

具体的,由于自然界岩石中所含的铀、钍及其衰变产物和钾的放射性同位素能引起地层的自然放射性。这些放射性元素在没有任何外界激发的情况下,可以释放出射线。这些具有放射性的元素在衰变过程中能发射α粒子、β粒子和γ射线。α粒子和β粒子的穿透能力很差,不能用于测井;而γ射线是一种具有很强穿透能力的高能电磁波,它能在钻井过程中被探测到。因此随钻自然伽马测量系统可以通过测量地层的自然伽马变化来反映地层的岩性。Specifically, the radioactive isotopes of uranium, thorium and their decay products and potassium contained in natural rocks can cause natural radioactivity in the formation. These radioactive elements emit radiation without any external excitation. These radioactive elements emit alpha, beta, and gamma rays during decay. Alpha particles and beta particles have poor penetrating ability and cannot be used for well logging; while gamma rays are high-energy electromagnetic waves with strong penetrating ability, which can be detected during drilling. Therefore, the natural gamma measurement while drilling system can reflect the lithology of the formation by measuring the natural gamma change of the formation.

可选的,深电阻率曲线是由深电阻率测探法(resistivity sounding,电测深法)获取的基准曲线。可选的,石油和水的电阻率相差很大,同样的储集层,含油时比含水时电阻率要高。Optionally, the deep resistivity curve is a reference curve obtained by deep resistivity sounding (resistivity sounding, electrical sounding method). Optionally, the resistivity of oil and water is very different, and the same reservoir has higher resistivity when it contains oil than when it contains water.

步骤13、对基准曲线和第一声波时差曲线之间建立响应关系方程,并进行多元线性拟合处理,得到无声波曲线井的第二声波时差曲线。Step 13 , establishing a response relationship equation between the reference curve and the first sonic time difference curve, and performing multivariate linear fitting processing to obtain a second sonic time difference curve of a well without an acoustic wave curve.

在本步骤中,基于步骤11和步骤12获取到的第一声波时差曲线和基准曲线,建立响应关系方程,进行多元线性拟合处理,得到无声波曲线井的声波时差曲线,即第二声波时差曲线。In this step, based on the first sonic time difference curve and the reference curve obtained in steps 11 and 12, a response relationship equation is established, and multivariate linear fitting is performed to obtain the sonic time difference curve of the non-acoustic curve well, that is, the second acoustic wave Time difference curve.

在一种可能的实现中,为了更加准确的进行井震标定,可以在基准曲线中加入密度曲线,参与到多元线性拟合处理中来,即响应关系方程可以是:DT=a*GR+b*DEN+c*lg(RD)+N。In a possible implementation, in order to perform well-seismic calibration more accurately, a density curve can be added to the reference curve to participate in the multivariate linear fitting process, that is, the response relationship equation can be: DT=a*GR+b *DEN+c*lg(RD)+N.

其中,DT表示声波时差,GR表示自然伽马,DEN表示密度,RD表示深电阻率,a、b、c和N分别为拟合常数。Among them, DT stands for acoustic transit time, GR stands for natural gamma, DEN stands for density, RD stands for deep resistivity, and a, b, c, and N are fitting constants, respectively.

可选的,密度曲线是用伽马源发射的伽马射线照射地层,根据康普顿效应测量地层体积密度,受井眼垮塌影响较大,但由于该曲线具有线性特征,很容易用自然伽马曲线校正,所以校正之后可与自然伽马曲线或者/和深电阻率曲线同时作为基准曲线参与到声波时差曲线的拟合处理中来。Optionally, the density curve irradiates the formation with gamma rays emitted by a gamma source, and measures the formation bulk density according to the Compton effect, which is greatly affected by borehole collapse. Horse curve correction, so after the correction, it can be used as a reference curve to participate in the fitting process of the acoustic time difference curve at the same time as the natural gamma curve or/and the deep resistivity curve.

步骤14、根据第二声波时差曲线,对无声波曲线井进行井震标定。Step 14 , according to the second sonic time difference curve, perform well seismic calibration for the well without sonic curve.

在本步骤中,根据多元线性拟合处理得到的第二声波时差曲线,对待标定的无声波曲线井进行井震标定。In this step, according to the second acoustic wave time difference curve obtained by the multivariate linear fitting process, the well seismic calibration is performed on the well without the acoustic wave curve to be calibrated.

本申请实施例提供的声波曲线井的井震标定方法,通过获取待标定的无声波曲线井的邻井的第一声波时差曲线以及无声波曲线井的基准曲线,对基准曲线和第一声波时差曲线之间建立响应关系方程,并进行多元线性拟合处理,得到无声波曲线井的第二声波时差曲线。最后根据第二声波时差曲线,对无声波曲线井进行井震标定。该技术方案中,通过待标定的无声波曲线井基准曲线与邻井的第一声波时差曲线,并进行多元线性拟合处理,实现了对无声波曲线井的井震标定,并提高了井震标定的效率。In the well-seismic calibration method for a sonic curve well provided by the embodiment of the present application, by obtaining the first sonic time difference curve of the adjacent well of the sonic curve well to be calibrated and the reference curve of the sonic curve well, the reference curve and the first acoustic curve are obtained. The response relationship equation is established between the wave time difference curves, and multivariate linear fitting is performed to obtain the second acoustic wave time difference curve of the well without acoustic wave curve. Finally, according to the second sonic time difference curve, the well seismic calibration is carried out for the well without sonic curve. In this technical scheme, by performing multivariate linear fitting processing on the reference curve of the non-acoustic curve well to be calibrated and the first acoustic time difference curve of the adjacent well, the well seismic calibration of the non-acoustic curve well is realized, and the wellbore is improved. Efficiency of vibration calibration.

在上述实施例的基础上,图2为本申请实施例提供的无声波曲线井的井震标定方法实施例二的流程图。如图2所示,上述步骤11可以通过如下步骤实现:On the basis of the foregoing embodiment, FIG. 2 is a flowchart of Embodiment 2 of the well seismic calibration method for a non-acoustic curve well provided by the embodiment of the present application. As shown in Figure 2, the above step 11 can be implemented by the following steps:

步骤21、根据无声波曲线井的特征,选取与无声波曲线井的距离小于预设距离且特征一致的测井作为邻井。Step 21 , according to the characteristics of the silent curve well, select the logging well whose distance from the silent curve well is less than the preset distance and has the same characteristics as the offset well.

在本步骤中,对于邻井的选取,应该以无声波曲线井的特征作为参考,选取特征一致(即沉积环境和压实程度相同)的测井。除了特征一致,为了保证地质条件等因素尽量接近,选取的邻井还应该与无声波曲线井的距离小于预设距离。In this step, for the selection of offset wells, the characteristics of the non-acoustic curve wells should be used as a reference, and the logging wells with the same characteristics (ie, the same depositional environment and the same degree of compaction) should be selected. In addition to the same characteristics, in order to ensure that the geological conditions and other factors are as close as possible, the distance between the selected offset well and the silent curve well should be smaller than the preset distance.

可选的,在实际选取邻井时,根据勘探开发不同阶段,该预设距离可以是几百米,也可以是几十千米,甚至更大。Optionally, when an offset well is actually selected, according to different stages of exploration and development, the preset distance may be several hundred meters, or several tens of kilometers, or even greater.

步骤22、从邻井中满足预设井眼条件的井段中选择样本层段。Step 22: Select a sample interval from the well intervals in the offset well that satisfy the preset wellbore condition.

在本步骤中,从邻井中选择满足预设井眼条件的井段中选择样本层段,井段的范围为距离样本层段纵向±50m的距离。其中,样本层段的范围可以是2米至10米。In this step, a sample interval is selected from the well intervals in the offset wells that meet the preset wellbore conditions, and the range of the well interval is a distance of ±50m longitudinally from the sample interval. Wherein, the range of the sample interval may be 2 meters to 10 meters.

可选的,预设井眼条件可以是该井段中的声波时差曲线不失真。Optionally, the preset wellbore condition may be that the acoustic transit time curve in the well section is not distorted.

步骤23、根据从样本层段中获取的多个样本点,确定第一声波时差曲线。Step 23: Determine a first acoustic wave time difference curve according to a plurality of sample points obtained from the sample slice.

在本步骤中,从上述步骤确定的样本层段中获取多个样本点,并根据这些样本点,确定出第一声波时差曲线。In this step, a plurality of sample points are obtained from the sample layers determined in the above steps, and according to these sample points, a first acoustic wave time difference curve is determined.

其中,多个样本点的个数范围为20个至80个,且多个样本点的岩性和流性的类型与无声波曲线井的部分井段数据相近。Among them, the number of multiple sample points ranges from 20 to 80, and the types of lithology and fluidity of the multiple sample points are similar to the data of part of the well section of the soundless curve well.

可选的,为了进一步保证井震标定的准确性,样本点的个数范围少于20个或多于80个都会导致拟合处理的数据不准确,影响标定结果。Optionally, in order to further ensure the accuracy of the well-seismic calibration, if the number of sample points is less than 20 or more than 80, the fitting data will be inaccurate and the calibration result will be affected.

本申请实施例提供的声波曲线井的井震标定方法,根据无声波曲线井的特征,选取与无声波曲线井的距离小于预设距离且特征一致的测井作为邻井,并从邻井中满足预设井眼条件的井段中选择样本层段,根据从样本层段中获取的多个样本点,确定第一声波时差曲线。该技术方案中,通过对邻井中多个样本点的选取,实现了更加准确的确定第一声波时差曲线,为后续的井震标定提高了可靠的基础。In the well-seismic calibration method for a sonic curve well provided by the embodiment of the present application, according to the characteristics of the silent curve well, a logging well whose distance from the silent curve well is less than a preset distance and has the same characteristics is selected as an offset well, and the offset wells satisfy the requirements from the offset well. A sample interval is selected from a well interval with preset wellbore conditions, and a first acoustic time difference curve is determined according to a plurality of sample points obtained from the sample interval. In this technical solution, by selecting multiple sample points in the offset well, a more accurate determination of the first acoustic time difference curve is realized, which provides a reliable basis for subsequent well-seismic calibration.

图3为本申请实施例提供的无声波曲线井的井震标定方法实施例三的流程图。如图3所示,该井震标定方法还可以包括如下步骤:FIG. 3 is a flowchart of Embodiment 3 of the well-seismic calibration method for a non-acoustic curve well provided by the embodiment of the present application. As shown in Figure 3, the well seismic calibration method may further include the following steps:

步骤31、获取采用第二声波时差曲线进行的井震标定结果与实测得到的地震波形之间的差异程度。Step 31: Obtain the degree of difference between the well-seismic calibration result using the second acoustic wave time difference curve and the seismic waveform obtained by actual measurement.

在本步骤中,当采用第二声波时差曲线进行井震标定得到井震标定结果后,判断井震标定结果是否准确。In this step, after the well-seismic calibration result is obtained by using the second acoustic wave transit time curve to obtain the well-seismic calibration result, it is judged whether the well-seismic calibration result is accurate.

在一种可能的实现中,利用实测得到的地震波形与井震标定结果进行对比,确定地震波形与井震标定结果之间的差异程度。In a possible implementation, the measured seismic waveforms are compared with the well-seismic calibration results to determine the degree of difference between the seismic waveforms and the well-seismic calibration results.

步骤32、若差异程度大于预设阈值,则重新从邻井中选择样本点,获取新的第一声波时差曲线。Step 32: If the degree of difference is greater than the preset threshold, re-select sample points from the adjacent wells to obtain a new first acoustic wave time difference curve.

在本步骤中,设定预设阈值,用于上述步骤中地震波形与井震标定结果之间的差异程度进行对比,若差异程度大于预设阈值,则认为该井震标定结果不准确。In this step, a preset threshold is set for comparing the degree of difference between the seismic waveform and the well-seismic calibration result in the above steps. If the difference is greater than the preset threshold, the well-seismic calibration result is considered to be inaccurate.

可选的,存在差异程度大于预设阈值的原因可以是样本点的选取不符合无声波曲线井的井段,在一种可能的实现中,需要重新从邻井中选择样本点,获取新的第一声波时差曲线,即重复执行上述步骤11。Optionally, the reason why the degree of difference is greater than the preset threshold may be that the selection of sample points does not conform to the well section of the silent curve well. For the sound wave time difference curve, repeat step 11 above.

可选的,若差异程度不大于预设阈值,则认为地震波形与井震标定结果符合程度高,即采用第二声波时差曲线进行的井震标定结果为准确的。Optionally, if the degree of difference is not greater than the preset threshold, it is considered that the seismic waveform and the well-seismic calibration result have a high degree of conformity, that is, the well-seismic calibration result using the second acoustic wave time difference curve is accurate.

步骤33、根据新的第一声波时差曲线和基准曲线进行拟合处理,得到新的第二声波时差曲线。Step 33: Perform fitting processing according to the new first acoustic wave time difference curve and the reference curve to obtain a new second acoustic wave time difference curve.

在本步骤中,获取到新的第一声波时差曲线之后,参照上述步骤13,与基准曲线进行拟合处理,得到新的第二声波时差曲线。In this step, after the new first acoustic wave time difference curve is acquired, referring to the above-mentioned step 13, a fitting process is performed with the reference curve to obtain a new second acoustic wave time difference curve.

在一种可能的实现中,根据新的第二声波时差曲线确定井震标定结果,重复步骤31的操作,直到地震波形与井震标定结果之间的差异程度不大于预设阈值,则认为该井震标定结果准确。In a possible implementation, the well-seismic calibration result is determined according to the new second acoustic wave time difference curve, and the operation of step 31 is repeated until the degree of difference between the seismic waveform and the well-seismic calibration result is not greater than a preset threshold, then it is considered that the The well-seismic calibration results are accurate.

本申请实施例提供的声波曲线井的井震标定方法,通过获取采用第二声波时差曲线进行的井震标定结果与实测得到的地震波形之间的差异程度。若差异程度大于预设阈值,则重新从邻井中选择样本点,获取新的第一声波时差曲线,基于该新的第一声波时差曲线,与基准曲线进行拟合处理,得到新的第二声波时差曲线。该技术方案中,通过第二声波时差曲线确定的井震标定结果与地震波形之间的差异程度,验证井震标定结果的准确性,若不准确,可重新确定第一声波时差曲线,避免了确定过程中样本点的选取可能存在的误差,进一步地提高了井震标定结果的准确性。The well-seismic calibration method for a sonic curve well provided by the embodiment of the present application obtains the degree of difference between the well-seismic calibration result using the second sonic time difference curve and the measured seismic waveform. If the degree of difference is greater than the preset threshold, re-select sample points from the adjacent wells, obtain a new first acoustic time difference curve, and perform fitting processing with the reference curve based on the new first acoustic time difference curve to obtain a new first acoustic time difference curve. Two sound wave time difference curve. In this technical solution, the degree of difference between the well-seismic calibration result and the seismic waveform determined by the second sonic time difference curve is used to verify the accuracy of the well-seismic calibration result. If it is not accurate, the first sonic time difference curve can be re-determined to avoid The possible errors in the selection of sample points in the determination process are eliminated, and the accuracy of the well-seismic calibration results is further improved.

在上述实施例的基础上,图4为本申请实施例提供的无声波曲线井的井震标定方法的总流程图。如图4所示,该方法包括如下步骤:On the basis of the above-mentioned embodiment, FIG. 4 is a general flow chart of the well-seismic calibration method of the non-acoustic curve well provided by the embodiment of the present application. As shown in Figure 4, the method includes the following steps:

第1步、获取邻井的多个样本点。Step 1. Obtain multiple sample points of the offset well.

第2步、获取无声波曲线井的基准曲线。Step 2: Obtain the reference curve of the silent curve well.

第3步、建立响应关系方程。The third step is to establish the response relationship equation.

第4步、多元线性拟合处理得到第二声波时差曲线。Step 4, multivariate linear fitting processing to obtain the second acoustic wave time difference curve.

第5步、对无声波曲线井进行井震标定。Step 5: Perform well-seismic calibration for the silent wave curve well.

第6步、确定井震标定结果与实测得到的地震波形之间的差异程度小于预设阈值,若是,则执行第7步;若否,则重新执行第1步。Step 6: Determine that the degree of difference between the well-seismic calibration result and the measured seismic waveform is less than a preset threshold, if yes, perform Step 7; if not, perform Step 1 again.

第7步、结束。Step 7, end.

本申请实施例提供的无声波曲线井的井震标定方法,根据邻井的多个样本点以及无声波曲线井的基准曲线,建立响应关系方程,并经过多元线性拟合处理得到第二声波时差曲线,根据第二声波时差曲线,进行井震标定。最后通过确定井震标定结果与实测得到的地震波形之间的差异程度,确定井震标定结果的准确性,该技术方案,实现了对无声波曲线井进行井震标定,并提高了井震标定的效率。In the well-seismic calibration method for a non-acoustic curve well provided by the embodiment of the present application, a response relationship equation is established according to multiple sample points of an adjacent well and a reference curve of a non-acoustic curve well, and a second acoustic time difference is obtained through multivariate linear fitting processing. Curve, according to the second acoustic wave time difference curve, the well-seismic calibration is carried out. Finally, by determining the degree of difference between the well-seismic calibration results and the measured seismic waveforms, the accuracy of the well-seismic calibration results is determined. This technical solution realizes the well-seismic calibration of the curve-free wells and improves the well-seismic calibration. s efficiency.

图5为本申请实施例提供的无声波曲线井的井震标定装置的结构示意图。如图5所示,该井震标定装置包括:获取模块51和处理模块52。FIG. 5 is a schematic structural diagram of a well-seismic calibration device for a non-acoustic curve well provided in an embodiment of the present application. As shown in FIG. 5 , the well-seismic calibration device includes: an acquisition module 51 and a processing module 52 .

处理模块52,用于获取待标定的无声波曲线井的邻井的第一声波时差曲线,邻井与无声波曲线井之间的距离小于预设值,且特征相同;The processing module 52 is configured to acquire the first acoustic time difference curve of the adjacent well of the non-acoustic curve well to be calibrated, the distance between the adjacent well and the non-acoustic curve well is less than a preset value, and the characteristics are the same;

获取模块51,用于获取无声波曲线井的基准曲线;an acquisition module 51, configured to acquire a reference curve of a soundless curve well;

处理模块52,还用于对基准曲线和第一声波时差曲线之间建立响应关系方程,并进行多元线性拟合处理,得到无声波曲线井的第二声波时差曲线;The processing module 52 is further configured to establish a response relationship equation between the reference curve and the first acoustic wave time difference curve, and perform multivariate linear fitting processing to obtain a second acoustic wave time difference curve of a well without an acoustic wave curve;

处理模块52,还用于根据第二声波时差曲线,对无声波曲线井进行井震标定。The processing module 52 is further configured to perform well-seismic calibration for a well without an acoustic wave curve according to the second acoustic wave time difference curve.

在本申请实施例一种可能设计中,处理模块52,用于获取待标定的无声波曲线井的邻井的第一声波时差曲线,包括:In a possible design of the embodiment of the present application, the processing module 52 is configured to obtain the first acoustic time difference curve of the adjacent well of the non-acoustic curve well to be calibrated, including:

根据无声波曲线井的特征,选取与无声波曲线井的距离小于预设距离且特征一致的测井作为邻井;According to the characteristics of the silent curve wells, the logging wells whose distance from the silent curve wells are less than the preset distance and have the same characteristics are selected as offset wells;

从邻井中满足预设井眼条件的井段中选择样本层段,样本层段的范围为2米至10米;The sample interval is selected from the well intervals in the offset well that meet the preset wellbore conditions, and the range of the sample interval is 2 meters to 10 meters;

根据从样本层段中获取的多个样本点,确定第一声波时差曲线。A first acoustic wave transit time curve is determined according to a plurality of sample points obtained from the sample slice.

在该种可能的设计中,处理模块52还包括从样本层段中获取多个样本点,多个样本点的个数范围为20个至80个,且多个样本点的岩性和流性的类型与无声波曲线井的部分井段数据相近。In this possible design, the processing module 52 further includes acquiring a plurality of sample points from the sample interval, the number of the plurality of sample points ranges from 20 to 80, and the lithology and fluidity of the plurality of sample points are The type is similar to the data of some sections of the silent curve well.

可选的,无声波曲线井和邻井的沉积环境和压实程度相同。Optionally, the depositional environment and the degree of compaction are the same for the silent curve well and the offset well.

可选的,基准曲线为自然伽马曲线或者/和深电阻率曲线。Optionally, the reference curve is a natural gamma curve or/and a deep resistivity curve.

可选的,响应关系方程为:DT=a*GR+b*DEN+c*lg(RD)+N;Optionally, the response relation equation is: DT=a*GR+b*DEN+c*lg(RD)+N;

其中,DT表示声波时差,GR表示自然伽马,DEN表示密度,RD表示深电阻率,a、b、c和N分别为拟合常数。Among them, DT stands for acoustic transit time, GR stands for natural gamma, DEN stands for density, RD stands for deep resistivity, and a, b, c, and N are fitting constants, respectively.

在本申请实施例另一种可能设计中,处理模块52,还用于:In another possible design of the embodiment of the present application, the processing module 52 is further configured to:

获取采用第二声波时差曲线进行的井震标定结果与实测得到的地震波形之间的差异程度;Obtain the degree of difference between the well-seismic calibration result using the second acoustic wave time difference curve and the measured seismic waveform;

若差异程度大于预设阈值,则重新从邻井中选择样本点,获取新的第一声波时差曲线;If the difference degree is greater than the preset threshold, re-select sample points from the offset wells to obtain a new first acoustic wave time difference curve;

根据新的第一声波时差曲线和基准曲线进行拟合处理,得到新的第二声波时差曲线。The fitting process is performed according to the new first acoustic wave time difference curve and the reference curve to obtain a new second acoustic wave time difference curve.

本申请实施例提供的无声波曲线井的井震标定装置,可用于执行上述实施例中的无声波曲线井的井震标定方法,其实现原理和技术效果类似,在此不再赘述。The well-seismic calibration device for a non-acoustic curve well provided by the embodiment of the present application can be used to execute the well-seismic calibration method for a non-acoustic curve well in the above-mentioned embodiment, and its realization principle and technical effect are similar, and are not repeated here.

需要说明的是,应理解以上装置的各个模块的划分仅仅是一种逻辑功能的划分,实际实现时可以全部或部分集成到一个物理实体上,也可以物理上分开。且这些模块可以全部以软件通过处理元件调用的形式实现;也可以全部以硬件的形式实现;还可以部分模块通过处理元件调用软件的形式实现,部分模块通过硬件的形式实现。例如,确定模块可以为单独设立的处理元件,也可以集成在上述装置的某一个芯片中实现,此外,也可以以程序代码的形式存储于上述装置的存储器中,由上述装置的某一个处理元件调用并执行以上确定模块的功能。其它模块的实现与之类似。此外这些模块全部或部分可以集成在一起,也可以独立实现。这里的处理元件可以是一种集成电路,具有信号的处理能力。在实现过程中,上述方法的各步骤或以上各个模块可以通过处理器元件中的硬件的集成逻辑电路或者软件形式的指令完成。It should be noted that it should be understood that the division of each module of the above apparatus is only a division of logical functions, and in actual implementation, all or part of it may be integrated into a physical entity, or may be physically separated. And these modules can all be implemented in the form of software calling through processing elements; they can also all be implemented in hardware; some modules can also be implemented in the form of calling software through processing elements, and some modules can be implemented in hardware. For example, the determination module may be a separately established processing element, or may be integrated into a certain chip of the above-mentioned device to be implemented, in addition, it may also be stored in the memory of the above-mentioned device in the form of program code, and a certain processing element of the above-mentioned device may Call and execute the function of the above determined module. The implementation of other modules is similar. In addition, all or part of these modules can be integrated together, and can also be implemented independently. The processing element here may be an integrated circuit with signal processing capability. In the implementation process, each step of the above-mentioned method or each of the above-mentioned modules can be completed by an integrated logic circuit of hardware in the processor element or an instruction in the form of software.

图6为本申请实施例提供的计算机设备的结构示意图。如图6所示,该计算机设备可以包括:处理器61、存储器62、收发器63。FIG. 6 is a schematic structural diagram of a computer device provided by an embodiment of the present application. As shown in FIG. 6 , the computer device may include: a processor 61 , a memory 62 , and a transceiver 63 .

其中,处理器61执行存储器62存储的计算机执行指令,使得处理器61执行上述实施例中的方案。处理器61可以是通用处理器,包括中央处理器CPU、网络处理器(networkprocessor,NP)等;还可以是数字信号处理器DSP、专用集成电路ASIC、现场可编程门阵列FPGA或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。The processor 61 executes the computer-executed instructions stored in the memory 62, so that the processor 61 executes the solutions in the foregoing embodiments. The processor 61 can be a general-purpose processor, including a central processing unit CPU, a network processor (NP), etc.; it can also be a digital signal processor DSP, an application-specific integrated circuit ASIC, a field programmable gate array FPGA or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.

收发器63用于和其他设备进行通信。可选的,在硬件实现上,上述图5所示实施例中的获取模块51对应于本实施例中的收发器63,该收发器63构成通信接口。The transceiver 63 is used to communicate with other devices. Optionally, in terms of hardware implementation, the obtaining module 51 in the embodiment shown in FIG. 5 above corresponds to the transceiver 63 in this embodiment, and the transceiver 63 constitutes a communication interface.

可选的,该计算机设备的上述各个器件之间可以通过系统总线连接。Optionally, the above-mentioned components of the computer equipment may be connected through a system bus.

本申请实施例提供的计算机设备,可用于执行上述实施例中的方案,其实现原理和技术效果类似,在此不再赘述。The computer equipment provided in the embodiments of the present application can be used to execute the solutions in the foregoing embodiments, and the implementation principles and technical effects thereof are similar, and details are not described herein again.

本申请实施例还提供一种运行指令的芯片,该芯片用于执行上述实施例中的方案。The embodiment of the present application further provides a chip for running instructions, and the chip is used to execute the solution in the above-mentioned embodiment.

本申请实施例还提供一种计算机可读存储介质,该计算机可读存储介质中存储有计算机指令,当该计算机指令在计算机上运行时,使得计算机执行上述实施例的方案。Embodiments of the present application further provide a computer-readable storage medium, where computer instructions are stored in the computer-readable storage medium, and when the computer instructions are executed on a computer, the computer can execute the solutions of the foregoing embodiments.

本申请实施例还提供一种计算机程序产品,该计算机程序产品包括计算机程序,其存储在计算机可读存储介质中,至少一个处理器可以从计算机可读存储介质读取计算机程序,至少一个处理器执行计算机程序时可实现上述实施例中的方案。Embodiments of the present application further provide a computer program product, where the computer program product includes a computer program, which is stored in a computer-readable storage medium, at least one processor can read the computer program from the computer-readable storage medium, and at least one processor The solutions in the above-described embodiments can be implemented when a computer program is executed.

最后应说明的是:以上各实施例仅用以说明本申请的技术方案,而非对其限制;尽管参照前述各实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本申请各实施例技术方案的范围。Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, but not to limit them; although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: The technical solutions described in the foregoing embodiments can still be modified, or some or all of the technical features thereof can be equivalently replaced; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the technical solutions of the embodiments of the present application. scope.

Claims (10)

1. A well-seismic calibration method for a silent wave curve well is characterized by comprising the following steps:
acquiring a first acoustic time difference curve of an adjacent well of a silent wave curve well to be calibrated, wherein the distance between the adjacent well and the silent wave curve well is smaller than a preset value, and the characteristics are the same;
acquiring a reference curve of the sonic-free curve well;
establishing a response relation equation between the reference curve and the first sound wave time difference curve, and performing multivariate linear fitting processing to obtain a second sound wave time difference curve of the silent wave curve well;
and carrying out well seismic calibration on the sonic-free curve well according to the second sound wave time difference curve.
2. The method of claim 1, wherein said obtaining a first sonic moveout profile of an adjacent well to a silent curve well to be calibrated comprises:
selecting a logging with the distance from the silent wave curve well being smaller than the preset distance and consistent in characteristics as the adjacent well according to the characteristics of the silent wave curve well;
selecting a sample interval from the well sections meeting the preset well bore conditions in the adjacent wells, wherein the range of the sample interval is 2-10 meters;
determining the first acoustic moveout curve from a plurality of sample points taken from the sample interval.
3. The method of claim 2, further comprising:
obtaining a plurality of sample points from the sample interval, wherein the number of the plurality of sample points ranges from 20 to 80, and the lithology and the fluidity of the plurality of sample points are similar to the partial well interval data of the silent wave curve well.
4. The method of any one of claims 1 to 3, wherein the deposition environment and compaction level of the anechoic curve well and the adjacent well are the same.
5. The method according to any one of claims 1 to 3, wherein the reference curve is a natural gamma curve or/and a deep resistivity curve.
6. The method of any one of claims 1 to 3, wherein the response relation equation is: DT ═ a × GR + b × DEN + c × lg (rd) + N;
where DT denotes acoustic moveout, GR denotes natural gamma, DEN denotes density, RD denotes deep resistivity, and a, b, c, and N are fitting constants, respectively.
7. The method according to any one of claims 1 to 3, further comprising:
acquiring the difference degree between the well seismic calibration result carried out by adopting the second sound wave time difference curve and the seismic waveform obtained by actual measurement;
if the difference degree is larger than a preset threshold value, selecting a sample point from the adjacent well again to obtain a new first acoustic time difference curve;
and fitting the new first sound wave time difference curve and the reference curve to obtain a new second sound wave time difference curve.
8. A well shake calibration device of a silent wave curve well is characterized by comprising: the device comprises an acquisition module and a processing module;
the processing module is used for acquiring a first acoustic time difference curve of an adjacent well of a silent wave curve well to be calibrated, wherein the distance between the adjacent well and the silent wave curve well is smaller than a preset value, and the characteristics of the adjacent well and the silent wave curve well are the same;
the acquisition module is used for acquiring a reference curve of the silent wave curve well;
the processing module is further configured to establish a response relation equation between the reference curve and the first acoustic time difference curve, and perform multivariate linear fitting processing to obtain a second acoustic time difference curve of the silent wave curve well;
and the processing module is also used for carrying out well seismic calibration on the sonic-free curve well according to the second sound wave time difference curve.
9. A computer device, comprising: a processor, a memory, and a transceiver;
the memory for storing computer program instructions executable on the processor;
the processor, when executing the computer program instructions, implements the method of any of claims 1-7 above.
10. A computer-readable storage medium having computer-executable instructions stored thereon, which when executed by a processor, perform the method of any one of claims 1-7.
CN202011439612.2A 2020-12-10 2020-12-10 Well-seismic calibration method, device, equipment and storage medium for sonic-free curve well Pending CN114624772A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011439612.2A CN114624772A (en) 2020-12-10 2020-12-10 Well-seismic calibration method, device, equipment and storage medium for sonic-free curve well

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011439612.2A CN114624772A (en) 2020-12-10 2020-12-10 Well-seismic calibration method, device, equipment and storage medium for sonic-free curve well

Publications (1)

Publication Number Publication Date
CN114624772A true CN114624772A (en) 2022-06-14

Family

ID=81896127

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011439612.2A Pending CN114624772A (en) 2020-12-10 2020-12-10 Well-seismic calibration method, device, equipment and storage medium for sonic-free curve well

Country Status (1)

Country Link
CN (1) CN114624772A (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104502996A (en) * 2014-12-23 2015-04-08 中国石油天然气集团公司 Density curve correction method and system
CN104656142A (en) * 2013-11-19 2015-05-27 中国石油天然气集团公司 Seismic horizon calibration method utilizing vertical seismic profiling (VSP) and well-logging combination
CN105759312A (en) * 2016-03-01 2016-07-13 中国海洋石油总公司 Reservoir characteristic method well-seismic calibration method
CN109932749A (en) * 2019-03-05 2019-06-25 中国石油天然气股份有限公司 Well-seismic calibration method and device
WO2020106287A1 (en) * 2018-11-21 2020-05-28 Halliburton Energy Services, Inc. Enhanced anisotropy analysis with multi-component dipole sonic data
CN111650646A (en) * 2020-05-26 2020-09-11 中国石油天然气集团有限公司 Acoustic curve correction method for seismic processing

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104656142A (en) * 2013-11-19 2015-05-27 中国石油天然气集团公司 Seismic horizon calibration method utilizing vertical seismic profiling (VSP) and well-logging combination
CN104502996A (en) * 2014-12-23 2015-04-08 中国石油天然气集团公司 Density curve correction method and system
CN105759312A (en) * 2016-03-01 2016-07-13 中国海洋石油总公司 Reservoir characteristic method well-seismic calibration method
WO2020106287A1 (en) * 2018-11-21 2020-05-28 Halliburton Energy Services, Inc. Enhanced anisotropy analysis with multi-component dipole sonic data
CN109932749A (en) * 2019-03-05 2019-06-25 中国石油天然气股份有限公司 Well-seismic calibration method and device
CN111650646A (en) * 2020-05-26 2020-09-11 中国石油天然气集团有限公司 Acoustic curve correction method for seismic processing

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
周韬 等: "拟声波曲线重构技术在储层预测中的应用", 石油地质与工程, vol. 26, no. 04, 25 July 2012 (2012-07-25), pages 52 - 54 *
张秀丽: "井约束地震反演技术在储层预测中的应用――以大庆油田杏56区块为例", 石油地球物理勘探, vol. 44, no. 1, 30 November 2009 (2009-11-30), pages 145 - 149 *
肖毅 等: "合成记录制作中井眼影响因素分析及对策", 煤田地质与勘探, vol. 37, no. 03, 25 June 2009 (2009-06-25), pages 64 - 67 *

Similar Documents

Publication Publication Date Title
US11802985B2 (en) Method and system for analyzing filling for karst reservoir based on spectrum decomposition and machine learning
US10877171B2 (en) Wavelet estimation for four-dimensional characterization of subsurface properties based on dynamic simulation
US9892366B2 (en) Facies definition using unsupervised classification procedures
EP2513674A2 (en) Workflow for petrophysical and geophysical formation evaluation of wireline and lwd log data
CN106761666B (en) A method and device for four-probe scattering gamma logging and nonlinear data inversion
CN107688037A (en) It is a kind of that the method for determining Rock in Well grading curve is distributed using nuclear magnetic resonance log T2
CN109283596A (en) A method for interpreting physical properties of carbonate reservoirs
US20210311223A1 (en) System and method for seismic inversion
CN110095811B (en) Method and device for constructing and processing velocity model of paste rock stratum
CN107450108A (en) The determination method and apparatus in dessert area
CN109061729B (en) A method for reconstructing gas-bearing sensitivity curves of high-temperature and high-pressure gas reservoirs
CN113435066B (en) Logging interpretation reservoir evaluation method based on digital core technology
CN107238859A (en) Calculate the method and system of virtual log synthetic seismogram
CN111502647B (en) Method and device for determining drilling geological environment factors and storage medium
CN114624772A (en) Well-seismic calibration method, device, equipment and storage medium for sonic-free curve well
CN113806674B (en) Quantification method and device for longitudinal scale of ancient river channel, electronic equipment and storage medium
CN114488304A (en) Lithology identification method and device and electronic equipment
Donald et al. Digital transformation of borehole sonic services in the oil and gas industry
CN111856571A (en) Method, device and equipment for converting velocity in logging frequency band to velocity in seismic frequency band
US20240272326A1 (en) Identification and quantification of subsurface geobodies
CN118501951B (en) Wellbore sonic imaging method, device, equipment and storage medium for offshore drilling
CN118091786B (en) Method and device for correcting linear density of cracks based on underground logging data
CN114428293B (en) Reservoir parameter random inversion method and device, electronic equipment and storage medium
Grilli et al. Sonic Log Quality Control: A Mandatory Workflow to Ensure Reliable Results in Different Disciplines
CN114428294A (en) Reservoir parameter statistical inversion method and device, electronic equipment and storage medium

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