CN117333575B - Shaft distributed optical fiber acoustic vibration monitoring data imaging method and processing terminal - Google Patents
Shaft distributed optical fiber acoustic vibration monitoring data imaging method and processing terminal Download PDFInfo
- Publication number
- CN117333575B CN117333575B CN202310742868.8A CN202310742868A CN117333575B CN 117333575 B CN117333575 B CN 117333575B CN 202310742868 A CN202310742868 A CN 202310742868A CN 117333575 B CN117333575 B CN 117333575B
- Authority
- CN
- China
- Prior art keywords
- data
- amplitude
- acoustic vibration
- sample
- segy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 42
- 238000003384 imaging method Methods 0.000 title claims abstract description 18
- 239000013307 optical fiber Substances 0.000 title claims abstract description 18
- 238000012545 processing Methods 0.000 title claims abstract description 16
- 239000011159 matrix material Substances 0.000 claims abstract description 23
- 230000008859 change Effects 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 238000005070 sampling Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 6
- 238000010586 diagram Methods 0.000 claims description 5
- 238000004519 manufacturing process Methods 0.000 description 15
- 239000012530 fluid Substances 0.000 description 11
- 238000013508 migration Methods 0.000 description 9
- 230000005012 migration Effects 0.000 description 9
- 230000000694 effects Effects 0.000 description 8
- 230000004044 response Effects 0.000 description 7
- 230000007774 longterm Effects 0.000 description 3
- 238000013500 data storage Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 1
- 101000827703 Homo sapiens Polyphosphoinositide phosphatase Proteins 0.000 description 1
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 1
- 102100023591 Polyphosphoinositide phosphatase Human genes 0.000 description 1
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
- G01H9/004—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/10—File systems; File servers
- G06F16/16—File or folder operations, e.g. details of user interfaces specifically adapted to file systems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/10—File systems; File servers
- G06F16/17—Details of further file system functions
- G06F16/1727—Details of free space management performed by the file system
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/10—File systems; File servers
- G06F16/17—Details of further file system functions
- G06F16/174—Redundancy elimination performed by the file system
- G06F16/1744—Redundancy elimination performed by the file system using compression, e.g. sparse files
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Human Computer Interaction (AREA)
- Fluid Mechanics (AREA)
- Environmental & Geological Engineering (AREA)
- Geophysics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域Technical Field
本发明涉及油气田开发井筒监测技术领域,尤其涉及一种井筒分布式光纤声振动监测数据成像方法及处理终端。The invention relates to the technical field of wellbore monitoring for oil and gas field development, and in particular to a wellbore distributed optical fiber acoustic vibration monitoring data imaging method and a processing terminal.
背景技术Background technique
分布式光纤声振动监测系统(DAS)现广泛应用于油气管道非法开挖监测、石油石化安防、边境线安防、光纤断电监测等,具有传输距离远、定位准、可连续分布式测量等优点。现在油气生产过程中的井筒生产状态监测也有应用DAS,根据井身结构将传感光缆通过油管送入井筒内,与井筒流体完全耦合,可以实现对井中是否出砂、流体如何运移以及井筒生产工具的活动状态进行监测。Distributed fiber acoustic vibration monitoring system (DAS) is now widely used in monitoring illegal excavation of oil and gas pipelines, petroleum and petrochemical security, border security, fiber power failure monitoring, etc. It has the advantages of long transmission distance, accurate positioning, and continuous distributed measurement. DAS is also used in the wellbore production status monitoring during the oil and gas production process. According to the wellbore structure, the sensor optical cable is sent into the wellbore through the oil pipe and fully coupled with the wellbore fluid. It can monitor whether sand is produced in the well, how the fluid moves, and the activity status of the wellbore production tools.
DAS对监测井中压裂活动等特殊现象较为敏感,因为压裂产生的声振动能量强,与周围背景声振动响应存在明显差异,在DAS监测结果上可以有较为直观的反映。但井筒多相流体运移的声振动响应规律并不清晰,其引起的声振动能量相对井下压裂等这类活动要低很多,DAS现场监测的原始数据无法直观判断井下的流体运移状态。而井筒多相流的运移状态又是指导现场生产的重要指标之一,有着十分重要的地位,现阶段,分布式光纤声振动监测系统在井筒多相流监测方面具有无可替代的地位,但是市场上比较成熟的DAS都为国外相关厂商生产,对于监测数据的处理和应用往往也受制于国外服务商的限制,有的需要到生产结果后再进行处理,时效性低。DAS is more sensitive to special phenomena such as fracturing activities in monitoring wells, because the acoustic vibration energy generated by fracturing is strong and there is a significant difference with the surrounding background acoustic vibration response, which can be reflected more intuitively in the DAS monitoring results. However, the acoustic vibration response law of wellbore multiphase fluid migration is not clear, and the acoustic vibration energy caused by it is much lower than that of downhole fracturing and other such activities. The original data of DAS on-site monitoring cannot intuitively judge the fluid migration state in the well. The migration state of wellbore multiphase flow is one of the important indicators to guide on-site production, and it has a very important position. At this stage, the distributed optical fiber acoustic vibration monitoring system has an irreplaceable position in wellbore multiphase flow monitoring, but the more mature DAS on the market are all produced by relevant foreign manufacturers. The processing and application of monitoring data are often subject to the restrictions of foreign service providers. Some need to be processed after the production results, which has low timeliness.
在海洋常规油气及非常规油气开发过程中,通常会在井筒内布设DAS,实现井筒中声振动能量的变化,进一步指示井下生产活动,常规的DAS监测可用于评价井筒中瞬时或者短周期内的工况变化,如井筒压裂等,但对于井筒内流体运移等长时间尺度非稳态现象的监测,由于数据量大,且数据信号中井筒工业噪声干扰严重,尚缺少实时显示的应用。In the process of conventional and unconventional offshore oil and gas development, DAS is usually deployed in the wellbore to realize the change of acoustic vibration energy in the wellbore and further indicate the downhole production activities. Conventional DAS monitoring can be used to evaluate the instantaneous or short-term operating condition changes in the wellbore, such as wellbore fracturing, but for the monitoring of long-term non-steady-state phenomena such as fluid migration in the wellbore, there is still a lack of real-time display applications due to the large amount of data and serious interference from wellbore industrial noise in the data signal.
发明内容Summary of the invention
针对现有技术中存在的不足,本发明的目的之一在于建立一种井筒分布式光纤声振动监测数据成像方法及处理终端,所指的监测数据是井筒内利用分布式光纤进行被动监测的数据,目的是服务于井筒流体运移等长时间尺度生产活动的刻画,可对海量的DAS数据进行针对性处理,处理后的DAS数据可实现井筒微振动现象的连续性显示,同时对背景噪声进行压制,能够实现海量监测数据的高信噪比快速成像,从而避免了现有DAS监测数据在井下生产监测功能单一、只能对短周期的生产活动进行解释的缺点。In view of the deficiencies in the prior art, one of the purposes of the present invention is to establish a wellbore distributed optical fiber acoustic vibration monitoring data imaging method and processing terminal. The monitoring data referred to is the data of passive monitoring in the wellbore using distributed optical fiber. The purpose is to serve the characterization of long-term production activities such as wellbore fluid migration. Massive DAS data can be processed in a targeted manner. The processed DAS data can realize the continuous display of wellbore micro-vibration phenomena, while suppressing background noise, and can achieve high signal-to-noise ratio rapid imaging of massive monitoring data, thereby avoiding the shortcomings of existing DAS monitoring data in underground production monitoring with a single function and only being able to explain short-term production activities.
为实现以上目的,本发明采取的技术方案:In order to achieve the above objectives, the technical solution adopted by the present invention is:
本发明提供了一种井筒分布式光纤声振动监测数据成像方法,包括以下步骤:The present invention provides a wellbore distributed optical fiber acoustic vibration monitoring data imaging method, comprising the following steps:
步骤一、对DAS监测到的k个segy文件数据进行预处理,以获取每一个segy文件的数据矩阵,所述数据矩阵的大小为N×M;Step 1: pre-process the k segy file data monitored by the DAS to obtain a data matrix of each segy file, where the size of the data matrix is N×M;
其中,M为样本数据道的道数,N=t/dt为每道样本数据道的样点总数,t为数据记录时间,dt为采样间隔,k为正整数,DAS为分布式光纤声振动监测系统;Wherein, M is the number of sample data channels, N=t/dt is the total number of sample points in each sample data channel, t is the data recording time, dt is the sampling interval, k is a positive integer, and DAS is a distributed optical fiber acoustic vibration monitoring system;
步骤二、所述步骤一获取的每一个segy文件的数据矩阵,取矩阵样本数据每一道中t秒内的N个样点数据,对其振幅做均方根振幅计算并记录下每一个振幅值,将每一个segy文件中M个振幅值进行集合,以形成每一条振幅能量变化曲线;Step 2: for each segy file data matrix obtained in step 1, take N sample point data within t seconds in each channel of matrix sample data, calculate the root mean square amplitude of the amplitude and record each amplitude value, and aggregate the M amplitude values in each segy file to form each amplitude energy change curve;
步骤三、所述步骤二获得的k条振幅能量变化曲线,按segy文件形成的时间顺序进行叠加排列,以形成长时间尺度下的井筒内声振动能量变化剖面图。Step 3: The k amplitude energy change curves obtained in step 2 are superimposed and arranged in the time sequence of the formation of the segy file to form a cross-sectional diagram of the acoustic vibration energy change in the wellbore on a long time scale.
进一步的,所述步骤二中所述的均方根振幅计算,其计算公式为:Furthermore, the root mean square amplitude calculation in step 2 is calculated as follows:
其中,AM为第m道的均方根振幅值,n为每道中样点数据的序列,Amn为第m道t秒内第n个样点的振幅值,M为样本数据道的道数,N=t/dt为每道样本数据道的样点总数,t为数据记录时间,dt为采样间隔。Wherein, AM is the RMS amplitude value of the mth channel, n is the sequence of sample data in each channel, A mn is the amplitude value of the nth sample within t seconds of the mth channel, M is the number of sample data channels, N = t/dt is the total number of sample points in each sample data channel, t is the data recording time, and dt is the sampling interval.
进一步的,所述预处理,包括:读取原始数据中标准segy格式的道头信息,基于所述道头信息,识别出道头中记录的采样间隔dt、数据记录时间t、样本数据道的道数M,确定每道样本数据道的样点总数为N=t/dt,以获取每一个segy文件的数据矩阵,所述数据矩阵的大小为N×M。Furthermore, the preprocessing includes: reading the header information in the standard segy format in the original data, and based on the header information, identifying the sampling interval dt, data recording time t, and the number of sample data tracks M recorded in the header, and determining that the total number of sample points in each sample data track is N=t/dt, so as to obtain the data matrix of each segy file, and the size of the data matrix is N×M.
本发明提供了一种实时处理终端,包括:The present invention provides a real-time processing terminal, comprising:
存储器,用于实时读写存储的程序指令;Memory, used to read and write stored program instructions in real time;
处理器,用于运行所述存储器存储的程序指令;A processor, configured to execute program instructions stored in the memory;
其中,所述处理器通过并行算法,以执行基于一种井筒分布式光纤声振动监测数据成像方法的步骤。The processor uses a parallel algorithm to execute steps based on a wellbore distributed optical fiber acoustic vibration monitoring data imaging method.
本发明的有益效果为:The beneficial effects of the present invention are:
(1)保护数据安全。(1) Protect data security.
(2)提高现场工作效率,支撑现场工程决策。由于segy格式数据文件占用存储空间大,生产平台现场往往按照30s或更短周期为一个振动状态剖面显示,短时间内井筒内声振动响应并不能直观反映流体运移的状态,同时长周期的井筒DAS监测数据是海量数据,但其生产结束后的处理工作繁重,使用本发明可以基于并行算法对海量数据在不影响监测质量的前提下,刻画长周期的井筒声振动响应结果,压缩数据量,优化数据存储空间,进一步提升数据处理效率,进而有效记录井筒流体振动响应变化。(2) Improve on-site work efficiency and support on-site engineering decision-making. Since the segy format data file occupies a large storage space, the production platform site is often displayed as a vibration state profile according to a period of 30s or less. The acoustic vibration response in the wellbore in a short period of time cannot intuitively reflect the state of fluid migration. At the same time, the long-period wellbore DAS monitoring data is massive data, but the processing work after the production is completed is arduous. The use of the present invention can be based on a parallel algorithm to describe the long-period wellbore acoustic vibration response results without affecting the monitoring quality, compress the data volume, optimize the data storage space, further improve the data processing efficiency, and effectively record the changes in the wellbore fluid vibration response.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1为井筒分布式光纤声振动监测数据成像方法的流程图;FIG1 is a flow chart of a method for imaging wellbore distributed optical fiber acoustic vibration monitoring data;
图2为单个文件数据均方根振幅计算及长时间尺度下数据成像方法图;Fig. 2 is a diagram showing the calculation of the root mean square amplitude of single file data and the data imaging method at a long time scale;
图3为图(a)短周期内原始声振动剖面实例与图(b)处理后2个月内声振动剖面实例的对比图;FIG3 is a comparison diagram of an original acoustic vibration profile example within a short period in FIG(a) and an acoustic vibration profile example within 2 months after processing in FIG(b);
图4为是本申请实施例提供的处理终端的示意图。FIG. 4 is a schematic diagram of a processing terminal provided in an embodiment of the present application.
具体实施方式Detailed ways
为使本发明的目的、技术方案及优点更加清楚、明确,下面结合附图和具体实施方式对本发明的内容做进一步详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部内容。In order to make the purpose, technical solution and advantages of the present invention clearer and more specific, the content of the present invention is further described in detail below in conjunction with the accompanying drawings and specific implementation methods. It is understood that the specific embodiments described herein are only used to explain the present invention, rather than to limit the present invention. It should also be noted that, for the convenience of description, only the parts related to the present invention are shown in the accompanying drawings, rather than all the contents.
如图1至图3所示,针对井筒多相流运移引起的微振动响应,结合DAS监测的数据,本发明提供了一种井筒分布式光纤声振动监测数据成像方法,包括以下步骤:As shown in FIGS. 1 to 3 , in view of the micro-vibration response caused by the migration of multiphase flow in the wellbore, combined with the data monitored by DAS, the present invention provides a wellbore distributed optical fiber acoustic vibration monitoring data imaging method, comprising the following steps:
本发明提供了一种井筒分布式光纤声振动监测数据成像方法,包括以下步骤:The present invention provides a wellbore distributed optical fiber acoustic vibration monitoring data imaging method, comprising the following steps:
步骤一、对DAS监测到的k个segy文件数据进行预处理,以获取每一个segy文件的数据矩阵,所述数据矩阵的大小为N×M,矩阵横向为道数递增方向,纵向为时间递增方向;Step 1: pre-process the k segy file data monitored by DAS to obtain a data matrix of each segy file, wherein the size of the data matrix is N×M, the horizontal direction of the matrix is the increasing direction of the number of channels, and the vertical direction is the increasing direction of time;
其中,M为样本数据道的道数,N=t/dt为每道样本数据道的样点总数,t为数据记录时间,dt为采样间隔,k为正整数,DAS为分布式光纤声振动监测系统;Wherein, M is the number of sample data channels, N=t/dt is the total number of sample points in each sample data channel, t is the data recording time, dt is the sampling interval, k is a positive integer, and DAS is a distributed optical fiber acoustic vibration monitoring system;
步骤二、所述步骤一获取的每一个segy文件的数据矩阵,取每一个矩阵样本数据中每一道t秒内的N个样点数据,对这些样点的振幅做均方根计算并记录下为一个新的振幅值,将每一个segy文件中M个计算后的振幅值进行集合,以形成一条新的振幅能量变化曲线;Step 2: for each segy file data matrix obtained in step 1, take N sample point data within t seconds in each matrix sample data, perform RMS calculation on the amplitude of these sample points and record them as a new amplitude value, and aggregate the M calculated amplitude values in each segy file to form a new amplitude energy change curve;
步骤三、所述步骤二获得的每一条振幅能量变化曲线,按k个segy文件形成的时间顺序进行叠加排列,实现海量数据的压缩和优化,以形成长时间尺度下的井筒声振动能量变化剖面图,从而实现井筒长周期声振动数据的成像;同时可实现井下多相振动状态实时显示,使井筒DAS监测到的数据得到更高效的利用,进而实现对井筒内生产工况的解释。Step three, each amplitude energy change curve obtained in the step two is superimposed and arranged in the time sequence formed by k segy files, so as to achieve compression and optimization of massive data, so as to form a wellbore acoustic vibration energy change profile on a long time scale, thereby realizing imaging of the long-period acoustic vibration data of the wellbore; at the same time, the real-time display of the multi-phase vibration state of the wellbore can be realized, so that the data monitored by the wellbore DAS can be used more efficiently, thereby realizing the interpretation of the production conditions in the wellbore.
本发明对海量的DAS数据进行针对性处理,处理后的DAS数据可实时连续性显示井筒微振动现象,由于系统性能引起的背景噪声强度往往低于实际信号的强度,且具有随机性,本发明通过对一批真实且长周期内有一定规律的实际信号进行叠加后,随机的背景噪声也相应的得到压制,提升了数据整体信噪比,通过查明井筒流体运移响应特征,避免了现有DAS监测数据主要用于井筒压裂监测的缺点,开辟了DAS系统监测方法新的应用领域。The present invention processes massive amounts of DAS data in a targeted manner. The processed DAS data can continuously display micro-vibration phenomena in the wellbore in real time. The intensity of background noise caused by system performance is often lower than the intensity of actual signals and is random. After superimposing a batch of real signals with certain regularity in a long period, the present invention suppresses the random background noise accordingly, thereby improving the overall signal-to-noise ratio of the data. By identifying the response characteristics of wellbore fluid migration, the shortcoming of existing DAS monitoring data being mainly used for wellbore fracturing monitoring is avoided, thus opening up a new application field for the DAS system monitoring method.
优选地,所述步骤二中所述的均方根振幅计算,其计算公式为:Preferably, the calculation formula for the root mean square amplitude calculation in step 2 is:
其中,AM为第m道的均方根振幅值,n为每道中样点数据的序列,Amn为第m道t秒内第n个样点的振幅值,M为样本数据道的道数,N=t/dt为每道样本数据道的样点总数,t为数据记录时间,dt为采样间隔。Wherein, AM is the RMS amplitude value of the mth channel, n is the sequence of sample data in each channel, A mn is the amplitude value of the nth sample within t seconds of the mth channel, M is the number of sample data channels, N = t/dt is the total number of sample points in each sample data channel, t is the data recording time, and dt is the sampling interval.
优选地,所述预处理,包括:读取原始数据中标准segy格式的道头信息,基于所述道头信息,识别出道头中记录的采样间隔dt、数据记录时间t、数据样本道的道数M,确定每道样本数据道的样点总数为N=t/dt,以获取每一个segy文件的数据矩阵,所述数据矩阵的大小为N×M。Preferably, the preprocessing includes: reading the header information in the standard segy format in the original data, and based on the header information, identifying the sampling interval dt, data recording time t, and the number of data sample tracks M recorded in the header, and determining that the total number of sample points in each sample data track is N=t/dt, so as to obtain the data matrix of each segy file, wherein the size of the data matrix is N×M.
如图4所示,本发明提供了一种实时处理终端,包括:As shown in FIG4 , the present invention provides a real-time processing terminal, comprising:
存储器,用于实时读写存储的程序指令;Memory, used to read and write stored program instructions in real time;
处理器,用于运行所述存储器存储的程序指令;A processor, configured to execute program instructions stored in the memory;
其中,所述处理器通过并行算法,以执行基于一种井筒分布式光纤声振动监测数据成像方法的步骤,例如图1、图2所示,能够实现海量监测数据的高信噪比快速成像。优选地,所述程序指令可以被分割成一个或多个模块/单元,所述一个或者多个模块/单元被存储在所述存储器中,并由所述处理器执行完成本发明。所述一个或多个模块/单元可以是能够完成特定功能的一系列程序指令段,该指令段用于描述所述程序指令在所述处理终端中的执行过程。优选地,所述处理终端可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述处理终端可包括,但不仅限于,处理器、存储器。Wherein, the processor uses a parallel algorithm to execute the steps based on a wellbore distributed optical fiber acoustic vibration monitoring data imaging method, such as shown in Figures 1 and 2, which can achieve high signal-to-noise ratio rapid imaging of massive monitoring data. Preferably, the program instructions can be divided into one or more modules/units, and the one or more modules/units are stored in the memory and executed by the processor to complete the present invention. The one or more modules/units can be a series of program instruction segments that can perform specific functions, and the instruction segments are used to describe the execution process of the program instructions in the processing terminal. Preferably, the processing terminal can be a computing device such as a desktop computer, a notebook, a PDA, and a cloud server. The processing terminal may include, but is not limited to, a processor and a memory.
本发明通过对海量监测数据的按时间进行划分,并在短周期下的数据中进行均方根平均,在不影响数据质量的前提大,有效压缩了数据体量,实现对长时间尺度下DAS数据的连续显示,此外还实现对系统随机干扰噪声的压制,其特点是噪声的振幅小,但不具有规律性,其压制方法就是上述提及的均方根叠加,使得随机噪声进一步被压制了;其数据计算处理方法简单,计算效率高,可直接作为现场数据处理成像软件中的独立模块使用,解决原始DAS数据对现场流体运移等相关长周期尺度的生产工况指示不清的问题,为现场常规油气及非常规油气开采活动提供了重要的数据指导,进一步明确了井筒内的生产流体的活动状态。The present invention divides the massive monitoring data by time and performs root mean square averaging on the data in a short period, thereby effectively compressing the data volume without affecting the data quality, and realizing continuous display of DAS data in a long time scale. In addition, the present invention also realizes the suppression of random interference noise of the system, which is characterized by a small amplitude of the noise but no regularity. The suppression method is the root mean square superposition mentioned above, so that the random noise is further suppressed; the data calculation and processing method is simple and has high calculation efficiency, and can be directly used as an independent module in the field data processing and imaging software, so as to solve the problem that the original DAS data cannot clearly indicate the production conditions of the long period scale related to the field fluid migration, etc., and provides important data guidance for the field conventional oil and gas and unconventional oil and gas production activities, and further clarifies the activity state of the production fluid in the wellbore.
需要说明的是,这里最终形成的声振动能量变化剖面时间样点数与采集形成的segy文件数量一致,虽然数据显示的时间间隔变大,但对于以月或年为周期的长时间连续监测数据来说,其采样时间间隔单位以毫秒变更为秒的影响可忽略不计,此外,该处理方法大大压缩了数据的存储空间,且不会影响整体工况的显示。It should be noted that the number of time sample points of the acoustic vibration energy change profile finally formed here is consistent with the number of segy files formed by the acquisition. Although the time interval of the data display becomes larger, for long-term continuous monitoring data with a period of months or years, the impact of changing the sampling time interval unit from milliseconds to seconds is negligible. In addition, this processing method greatly compresses the data storage space and will not affect the display of the overall working conditions.
上述实施例只是为了说明本发明的技术构思及特点,其目的是在于让本领域内的普通技术人员能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡是根据本发明内容的实质所做出的等效的变化或修饰,都应涵盖在本发明的保护范围内。The above embodiments are only for illustrating the technical concept and features of the present invention, and their purpose is to enable ordinary technicians in the field to understand the content of the present invention and implement it accordingly, and they cannot be used to limit the protection scope of the present invention. Any equivalent changes or modifications made based on the essence of the content of the present invention should be included in the protection scope of the present invention.
Claims (3)
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310742868.8A CN117333575B (en) | 2023-06-21 | 2023-06-21 | Shaft distributed optical fiber acoustic vibration monitoring data imaging method and processing terminal |
PCT/CN2023/129895 WO2024259874A1 (en) | 2023-06-21 | 2023-11-06 | Distributed acoustic sensing data imaging method for shaft, and processing terminal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310742868.8A CN117333575B (en) | 2023-06-21 | 2023-06-21 | Shaft distributed optical fiber acoustic vibration monitoring data imaging method and processing terminal |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117333575A CN117333575A (en) | 2024-01-02 |
CN117333575B true CN117333575B (en) | 2024-05-31 |
Family
ID=89289112
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310742868.8A Active CN117333575B (en) | 2023-06-21 | 2023-06-21 | Shaft distributed optical fiber acoustic vibration monitoring data imaging method and processing terminal |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN117333575B (en) |
WO (1) | WO2024259874A1 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN119290128B (en) * | 2024-10-15 | 2025-06-20 | 华北科技学院(中国煤矿安全技术培训中心) | Method for extracting velocity spectrum from DAS data, readable storage medium and electronic device |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106199702A (en) * | 2016-07-19 | 2016-12-07 | 中国地质大学(北京) | Geological data storage tape falls the anomalous amplitude drawing method that magnetic powder produces |
CN106610505A (en) * | 2016-12-29 | 2017-05-03 | 中国石油大学(华东) | Logging and seismic data matching method based on combination of DTW (Dynamic Time Warping) and ABA (Amplitude Balance Adjustment) |
CN113219539A (en) * | 2020-02-05 | 2021-08-06 | 中国石油天然气集团有限公司 | Method and device for removing borehole wave interference in optical fiber acoustic sensing seismic data |
CN113962304A (en) * | 2021-10-22 | 2022-01-21 | 中海石油(中国)有限公司 | Multilayer simultaneous intelligent seismic horizon interpretation method and system |
CN114993269A (en) * | 2022-04-13 | 2022-09-02 | 南方海洋科学与工程广东省实验室(广州) | A Method of Extracting Track Water Depth Using SEGY Data |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106204497B (en) * | 2016-07-20 | 2018-12-25 | 长安大学 | A kind of pavement crack extraction algorithm based on smooth smoothed curve and matched curve |
US11467308B2 (en) * | 2018-05-21 | 2022-10-11 | West Virginia University | Fibro: a fiber optic data processing software for unconventional reservoirs |
-
2023
- 2023-06-21 CN CN202310742868.8A patent/CN117333575B/en active Active
- 2023-11-06 WO PCT/CN2023/129895 patent/WO2024259874A1/en active Application Filing
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106199702A (en) * | 2016-07-19 | 2016-12-07 | 中国地质大学(北京) | Geological data storage tape falls the anomalous amplitude drawing method that magnetic powder produces |
CN106610505A (en) * | 2016-12-29 | 2017-05-03 | 中国石油大学(华东) | Logging and seismic data matching method based on combination of DTW (Dynamic Time Warping) and ABA (Amplitude Balance Adjustment) |
CN113219539A (en) * | 2020-02-05 | 2021-08-06 | 中国石油天然气集团有限公司 | Method and device for removing borehole wave interference in optical fiber acoustic sensing seismic data |
CN113962304A (en) * | 2021-10-22 | 2022-01-21 | 中海石油(中国)有限公司 | Multilayer simultaneous intelligent seismic horizon interpretation method and system |
CN114993269A (en) * | 2022-04-13 | 2022-09-02 | 南方海洋科学与工程广东省实验室(广州) | A Method of Extracting Track Water Depth Using SEGY Data |
Non-Patent Citations (3)
Title |
---|
Ground motion response to an ML 4.3 earthquakr using co-located distributed acoustic sensing and seismometer arrays;Herberf F. Wang .etal;Geophysical Journal International;20180317;全文 * |
Spark框架下地震属性处理方法研究;朱丽萍 等;计算机与数字工程;20180831;第46卷(第08期);全文 * |
三维速度场建立方法;韩宏伟 等;油气地质与采收率;20100131(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN117333575A (en) | 2024-01-02 |
WO2024259874A1 (en) | 2024-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN117333575B (en) | Shaft distributed optical fiber acoustic vibration monitoring data imaging method and processing terminal | |
EP3028038A1 (en) | Vibration condition monitoring system and associated methods | |
Li et al. | Adaptive identification of time-varying cable tension based on improved variational mode decomposition | |
CN118656586A (en) | Vulnerability analysis and toughness assessment method and system for substation structures under strong earthquakes | |
Elichev et al. | Understanding well events with machine learning | |
Yang et al. | Vibration prediction and analysis of the main beam of the TBM based on a multiple linear regression model | |
CN119558098A (en) | A simulation method and system for predicting mechanical parameters of rock mass structural surface at base-cover interface | |
Greenwood | Vibration monitoring and mitigation-an integrated measurement system | |
CN118859336A (en) | A dam safety monitoring method and system | |
Wang et al. | Advances in vibration analysis and modeling of large rotating mechanical equipment in mining arena: A review | |
CN112801033A (en) | AlexNet network-based construction disturbance and leakage identification method along long oil and gas pipeline | |
CN117927218A (en) | Liquid production profile interpretation method and device based on distributed optical fibers | |
Lim et al. | Diagnosis for loose blades in gas turbines using wavelet analysis | |
Clemens | Sucker rod pumping system optimisation based on a novel finite elements analysis | |
Shi-Yi et al. | Analyzing transient pressure from permanent downhole gauges (PDG) using wavelet method | |
Goswami et al. | Data-driven sensor selection for industrial gearbox fault diagnosis using principal component analysis | |
Zhang et al. | Full-Azimuth Quantitative Identification of Near-Fault Pulse-Like Ground Motions Based on M&P Wavelet | |
Molotkov et al. | ESP Monitoring and Diagnosis Enables Operational Efficiency | |
Ma et al. | Research on acoustic emission precursor characteristics of coal sample unloading failure based on discrete wavelet analysis | |
Suzuki et al. | Automatic detection of pressure-buildup intervals from permanent downhole pressure data using filter convolution | |
Vekved et al. | Optimizing Rod Pump Performance Using High Frequency, High Resolution Pressure Data at the Wellhead | |
Li et al. | Research on Pipeline Damage Identification Based on Optical Fiber Sensing and Deep Learning | |
Ning et al. | Research on dynamic analysis method of pier safety based on time-frequency domain combination | |
CN118569748A (en) | Underground pipe network modeling and power grid infrastructure engineering dynamic early warning system | |
Tserng et al. | Measuring the Natural Frequencies of Bridges by Wireless Sensors to Monitor Foundation Scouring |
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 |