[go: up one dir, main page]

CN115994338B - A method, device, electronic device and medium for fusion of well-seismic information - Google Patents

A method, device, electronic device and medium for fusion of well-seismic information

Info

Publication number
CN115994338B
CN115994338B CN202111210323.XA CN202111210323A CN115994338B CN 115994338 B CN115994338 B CN 115994338B CN 202111210323 A CN202111210323 A CN 202111210323A CN 115994338 B CN115994338 B CN 115994338B
Authority
CN
China
Prior art keywords
structure tensor
field
velocity model
interpolation
logging data
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
Application number
CN202111210323.XA
Other languages
Chinese (zh)
Other versions
CN115994338A (en
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.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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 China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN202111210323.XA priority Critical patent/CN115994338B/en
Publication of CN115994338A publication Critical patent/CN115994338A/en
Application granted granted Critical
Publication of CN115994338B publication Critical patent/CN115994338B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The application discloses a well earthquake information fusion method, a well earthquake information fusion device, electronic equipment and a medium. The method comprises the steps of obtaining acoustic time difference logging data and an offset profile, carrying out median filtering and Gaussian filtering on the acoustic time difference logging data, extracting a structure tensor and a coherent body from the offset profile, calculating a structure tensor field by using the structure tensor and the coherent body, carrying out interpolation processing on the logging data by combining the structure tensor field by using a mixed neighborhood method to obtain an interpolation speed model and a diffusion time field, calculating a weighting coefficient by using the diffusion time field, and fusing the chromatography speed model and the interpolation speed model to obtain a fusion speed model. The application improves the precision and resolution of the chromatographic velocity model by using a small amount of logging data so as to efficiently obtain more accurate structural imaging results and guide the oil-gas exploration and development.

Description

Well earthquake information fusion method and device, electronic equipment and medium
Technical Field
The invention relates to the technical field of oil-gas seismic exploration, in particular to a well seismic information fusion method, a well seismic information fusion device, electronic equipment and a medium.
Background
In the field of oil-gas seismic exploration, a velocity modeling technology is the basis of subsequent seismic wave migration imaging and reservoir inversion, ray tomography inversion is used as the current mainstream velocity modeling technology, a relatively accurate low-frequency velocity model can be obtained, but a high wave number component of the velocity model is difficult to obtain, a wave tomography inversion method is generated, a more accurate velocity model can be obtained, in practical application, the problems of inversion efficiency, dependence on an initial model, data quality and the like still exist, the method is difficult to be widely applied, a set of methods between a conventional ray theory and a wave equation theory are developed in the industry, gaussian beam tomography realizes beam tomography of linear approximation of a wave equation by filtering, compared with the traditional ray tomography, the Gaussian beam tomography increases the coverage of rays, reduces the condition number of a nuclear matrix, reduces the disease state of a chromatographic equation set, and can obtain a more accurate velocity model, meanwhile, in a region to be researched, a small amount of well measurement data often exist, the acoustic time difference can accurately reflect the velocity information, and often serve as constraint conditions to be added into the chromatographic inversion system, the method can be limited by the constraint method, the method can be fast in the condition of the constraint method, but the well drilling accuracy is limited by the constraint on the condition of the constraint method can be greatly reduced, the accuracy of the well drilling information can be greatly reduced, and the method can not be limited by the constraint on the real-time information.
Aiming at the problems, the well earthquake information fusion method is expected to be provided, the accuracy and the resolution of a chromatographic velocity model are improved by using a small amount of logging data, so that more accurate structural imaging results are obtained efficiently, and the oil and gas exploration and development are guided.
The information disclosed in the background section of the invention is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.
Disclosure of Invention
The invention provides a well earthquake information fusion method, a well earthquake information fusion device, electronic equipment and a well earthquake information fusion medium, which at least solve the technical problems that the accuracy and the resolution of a logging data chromatographic velocity model are low, and an accurate structural imaging result cannot be obtained efficiently.
In a first aspect, an embodiment of the present disclosure provides a method for fusion of well seismic information, including:
acquiring acoustic time difference logging data and an offset profile;
Performing median filtering and Gaussian filtering on the acoustic time difference logging data;
extracting a structure tensor and a coherent body from the offset profile, and calculating a structure tensor field by using the structure tensor and the coherent body;
Performing interpolation processing on logging data by using a mixed neighborhood method and combining a structure tensor field to obtain an interpolation speed model and a diffusion time field;
and calculating a weighting coefficient by using the diffusion time field, and fusing the chromatographic velocity model and the interpolation velocity model to obtain a fusion velocity model.
As a specific implementation manner of the embodiments of the present disclosure, the calculation formula of the structural tensor field is as follows:
Where S is the structure tensor extracted from the offset profile, c is the coherence volume extracted from the offset profile, and D is the structure tensor field.
As a specific implementation manner of the embodiment of the disclosure, the interpolation formula is as follows:
τ(x)=0,X∈χ
wherein τ is a diffusion time field, D is a structure tensor field, χ is an acoustic time difference logging position area, p is an acoustic time difference logging value with the shortest diffusion time, and q is an interpolation result.
As a specific implementation manner of the embodiment of the present disclosure, the fused calculation formula is as follows:
vmix=l·vinterp+(1-l)·vtomo
Wherein T max is a time threshold, τ is a diffusion time field, l is a weighting coefficient, v interp is an interpolation velocity model, and v tomo is a chromatography velocity model.
In a second aspect, an embodiment of the present disclosure further provides a device for fusion of borehole seismic information, including:
the acquisition module is used for acquiring acoustic time difference logging data and an offset profile;
the processing module is used for carrying out median filtering and Gaussian filtering on the acoustic time difference logging data;
the calculation module is used for extracting a structure tensor and a coherent body from the offset profile and calculating a structure tensor field by using the structure tensor and the coherent body;
The interpolation processing module is used for carrying out interpolation processing on the logging data by utilizing a mixed neighborhood method and combining the structure tensor field to obtain an interpolation speed model and a diffusion time field;
and the fusion module is used for calculating a weighting coefficient by using the diffusion time field, fusing the chromatographic velocity model and the interpolation velocity model, and obtaining a fusion velocity model.
As a specific implementation manner of the embodiments of the present disclosure, the calculation formula of the structural tensor field is as follows:
Where S is the structure tensor extracted from the offset profile, c is the coherence volume extracted from the offset profile, and D is the structure tensor field.
As a specific implementation manner of the embodiment of the disclosure, the interpolation formula is as follows:
τ(x)=0,x∈χ
wherein τ is a diffusion time field, D is a structure tensor field, x is an acoustic time difference logging position area, p is an acoustic time difference logging value with the shortest diffusion time, and q is an interpolation result.
As a specific implementation manner of the embodiment of the present disclosure, the fused calculation formula is as follows:
vmix=l·vinterp+(1-l)·vtomo
Wherein T max is a time threshold, τ is a diffusion time field, l is a weighting coefficient, v interp is an interpolation velocity model, and v tomo is a chromatography velocity model.
In a third aspect, embodiments of the present disclosure further provide an electronic device, including:
At least one processor, and
A memory communicatively coupled to the at least one processor, wherein,
The memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of well intervention information fusion as described above.
In a fourth aspect, the presently disclosed embodiments also provide a non-transitory computer readable storage medium storing computer instructions for causing a computer to perform a well-shock information fusion method as described above.
The invention has the beneficial effects that:
The method utilizes a mixed neighborhood method to interpolate a small amount of acoustic time difference logging data, adopts interpolation precision parameters to weight the multi-source velocity model, recovers high-frequency information of the velocity model to a certain extent, improves the precision and resolution of the chromatographic velocity model, realizes accurate imaging, and is simple and direct, and has higher efficiency.
According to the invention, the mixed neighborhood method is used for carrying out interpolation processing on logging data in combination with the structure tensor field, the diffusion time field is used for calculating the weighting coefficient, the chromatographic velocity model and the interpolation velocity model are fused, the fusion velocity model is obtained, the high-frequency information of the chromatographic velocity model is recovered to a certain extent, the resolution ratio of the velocity model is improved, a foundation is laid for subsequent high-precision imaging, and the method is simple and direct, and can realize rapid modeling imaging.
The method and apparatus of the present invention have other features and advantages which will be apparent from or are set forth in detail in the accompanying drawings and the following detailed description, which are incorporated herein, and which together serve to explain certain principles of the present invention.
Drawings
The foregoing and other objects, features and advantages of the invention will be apparent from the following more particular descriptions of exemplary embodiments of the invention as illustrated in the accompanying drawings wherein like reference numbers generally represent like parts throughout the exemplary embodiments of the invention.
FIG. 1 is a flow chart of a method of fusion of well seismic information according to embodiment 1;
FIG. 2 is a diagram of marmousi standard models employed in example 1;
FIG. 3 is a Gaussian beam tomography velocity field diagram of example 1;
FIG. 4 is a velocity model diagram of the mixed domain interpolation in example 1;
FIG. 5 is a graph of the velocity model after fusion in example 1;
FIG. 6 is a graph showing the comparison of the true velocity values, the chromatographic velocity model and the fusion velocity model in example 1;
fig. 7 is a block diagram of a well seismic information fusion apparatus according to embodiment 2 of the present invention.
Detailed Description
Preferred embodiments of the present invention will be described in more detail below. While the preferred embodiments of the present invention are described below, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein.
The invention provides a well earthquake information fusion method, which comprises the following steps:
acquiring acoustic time difference logging data and an offset profile;
Performing median filtering and Gaussian filtering on the acoustic time difference logging data;
extracting a structure tensor and a coherent body from the offset profile, and calculating a structure tensor field by using the structure tensor and the coherent body;
Performing interpolation processing on logging data by using a mixed neighborhood method and combining a structure tensor field to obtain an interpolation speed model and a diffusion time field;
and calculating a weighting coefficient by using the diffusion time field, and fusing the chromatographic velocity model and the interpolation velocity model to obtain a fusion velocity model.
In one example, the structural tensor field is calculated as follows:
Where S is the structure tensor extracted from the offset profile, c is the coherence volume extracted from the offset profile, and D is the structure tensor field.
In one example, the interpolation formula is as follows:
τ(x)=0,x∈χ
wherein τ is a diffusion time field, D is a structure tensor field, χ is an acoustic time difference logging position area, p is an acoustic time difference logging value with the shortest diffusion time, and q is an interpolation result.
In one example, the fused calculation formula is as follows:
vmix=l·vinterp+(1-l)·vtomo
Wherein T max is a time threshold, τ is a diffusion time field, l is a weighting coefficient, v interp is an interpolation velocity model, and v tomo is a chromatography velocity model.
The invention also provides a well earthquake information fusion device, which comprises:
the acquisition module is used for acquiring acoustic time difference logging data and an offset profile;
the processing module is used for carrying out median filtering and Gaussian filtering on the acoustic time difference logging data;
the calculation module is used for extracting a structure tensor and a coherent body from the offset profile and calculating a structure tensor field by using the structure tensor and the coherent body;
The interpolation processing module is used for carrying out interpolation processing on the logging data by utilizing a mixed neighborhood method and combining the structure tensor field to obtain an interpolation speed model and a diffusion time field;
and the fusion module is used for calculating a weighting coefficient by using the diffusion time field, fusing the chromatographic velocity model and the interpolation velocity model, and obtaining a fusion velocity model.
In one example, the structural tensor field is calculated as follows:
Where S is the structure tensor extracted from the offset profile, c is the coherence volume extracted from the offset profile, and D is the structure tensor field.
In one example, the interpolation formula is as follows:
τ(x)=0,x∈χ
wherein τ is a diffusion time field, D is a structure tensor field, χ is an acoustic time difference logging position area, p is an acoustic time difference logging value with the shortest diffusion time, and q is an interpolation result.
In one example, the fused calculation formula is as follows:
vmix=l·vinterp+(1-l)·vtomo
Wherein T max is a time threshold, τ is a diffusion time field, l is a weighting coefficient, v interp is an interpolation velocity model, and v tomo is a chromatography velocity model.
The present invention also provides an electronic device including:
and a memory communicatively coupled to the at least one processor, wherein,
The memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of well intervention information fusion as described above.
The present invention also provides a non-transitory computer readable storage medium storing computer instructions for causing a computer to perform a well shock information fusion method as described above.
In order to facilitate understanding of the solution and the effects of the embodiments of the present invention, four specific application examples are given below. It will be understood by those of ordinary skill in the art that the examples are for ease of understanding only and that any particular details thereof are not intended to limit the present invention in any way.
Example 1
FIG. 1 shows a flow chart of the steps of a method of fusion of well seismic information, as shown in FIG. 1, according to one embodiment of the invention, comprising:
s01, acquiring acoustic time difference logging data and an offset profile;
S02, carrying out median filtering and Gaussian filtering on the acoustic time difference logging data;
S03, extracting a structure tensor and a coherent body from the offset profile, and calculating a structure tensor field by using the structure tensor and the coherent body;
S04, carrying out interpolation processing on logging data by utilizing a mixed neighborhood method and combining a structure tensor field to obtain an interpolation speed model and a diffusion time field;
And S05, calculating a weighting coefficient by using the diffusion time field, and fusing the chromatographic velocity model and the interpolation velocity model to obtain a fusion velocity model.
Performing interpolation processing on logging data by combining a structure tensor field by using a mixed neighborhood method, calculating a weighting coefficient by using a diffusion time field, fusing a chromatographic velocity model and an interpolation velocity model to obtain a fusion velocity model, the method has the advantages that the high-frequency information of the chromatographic velocity model is recovered to a certain extent, the resolution of the velocity model is improved, a foundation is laid for subsequent high-precision imaging, and the method is simple and direct and can realize rapid modeling imaging.
Fig. 2 is a marmousi standard model diagram adopted in the embodiment of the present invention, a gaussian beam chromatography velocity model in the embodiment of the present invention of fig. 3 is selected as an example for verification, a geophysical standard model marmousi model is selected for verification, a velocity model is obtained by using reflected wave gaussian beam chromatography, real velocity values at positions x=200m and x=1100 m are taken as logging data, a structure tensor and a coherent body are extracted from an offset section, the logging data is interpolated by using a hybrid neighborhood method, and the interpolated velocity model is obtained, as shown in fig. 4. The interpolated velocity model has a higher frequency than the tomographic velocity model, with velocity values near the well location closer to the true velocity model, and velocity values further from the well location offset from the true value. Setting the time threshold as 2000ms, calculating a weighting coefficient by using a diffusion time field, fusing the two velocity models to obtain a fused velocity model, taking X=500m and X=1500m data, comparing the real velocity value, the chromatographic velocity model and the fused velocity model, and displaying the result that the fused velocity model is closer to the real velocity value, as shown in fig. 6.
Example 2
FIG. 7 illustrates a borehole seismic information fusion apparatus in accordance with one embodiment of the invention.
As shown in fig. 7, the well-shock information fusion device includes:
the acquisition module is used for acquiring acoustic time difference logging data and an offset profile;
the processing module is used for carrying out median filtering and Gaussian filtering on the acoustic time difference logging data;
the calculation module is used for extracting a structure tensor and a coherent body from the offset profile and calculating a structure tensor field by using the structure tensor and the coherent body;
The interpolation processing module is used for carrying out interpolation processing on the logging data by utilizing a mixed neighborhood method and combining the structure tensor field to obtain an interpolation speed model and a diffusion time field;
and the fusion module is used for calculating a weighting coefficient by using the diffusion time field, fusing the chromatographic velocity model and the interpolation velocity model, and obtaining a fusion velocity model.
As a specific implementation manner of the embodiments of the present disclosure, the calculation formula of the structural tensor field is as follows:
Where S is the structure tensor extracted from the offset profile, c is the coherence volume extracted from the offset profile, and D is the structure tensor field.
As a specific implementation manner of the embodiment of the disclosure, the interpolation formula is as follows:
τ(x)=0,x∈χ
wherein τ is a diffusion time field, D is a structure tensor field, χ is an acoustic time difference logging position area, p is an acoustic time difference logging value with the shortest diffusion time, and q is an interpolation result.
As a specific implementation manner of the embodiment of the present disclosure, the fused calculation formula is as follows:
vmix=l·vinterp+(1-l)·vtomo
Wherein T max is a time threshold, τ is a diffusion time field, l is a weighting coefficient, v interp is an interpolation velocity model, and v tomo is a chromatography velocity model.
Example 3
The present disclosure provides an electronic device comprising at least one processor, and a memory communicatively coupled to the at least one processor, wherein,
The memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of seismic information fusion described above.
An electronic device according to an embodiment of the present disclosure includes a memory and a processor.
The memory is for storing non-transitory computer readable instructions. In particular, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and/or non-volatile memory. The volatile memory may include, for example, random Access Memory (RAM) and/or cache memory (cache), and the like. The non-volatile memory may include, for example, read Only Memory (ROM), hard disk, flash memory, and the like.
The processor may be a Central Processing Unit (CPU) or other form of processing unit having data processing and/or instruction execution capabilities, and may control other components in the electronic device to perform the desired functions. In one embodiment of the present disclosure, the processor is configured to execute the computer readable instructions stored in the memory.
It should be understood by those skilled in the art that, in order to solve the technical problem of how to obtain a good user experience effect, the present embodiment may also include well-known structures such as a communication bus, an interface, and the like, and these well-known structures are also included in the protection scope of the present disclosure.
The detailed description of the present embodiment may refer to the corresponding description in the foregoing embodiments, and will not be repeated herein.
Example 4
Embodiments of the present disclosure provide a non-transitory computer readable storage medium storing computer instructions for causing a computer to perform a well-shock information fusion method as described above.
The non-transitory computer readable storage medium according to embodiments of the present disclosure, when executed by a processor, performs all or part of the steps of the methods of embodiments of the present disclosure described above.
Such computer readable storage media include, but are not limited to, optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., tape or removable hard disk), media with built-in rewritable non-volatile memory (e.g., memory card), and media with built-in ROM (e.g., ROM cartridge).
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention has been given for the purpose of illustrating the benefits of embodiments of the invention only and is not intended to limit embodiments of the invention to any examples given.
The foregoing description of embodiments of the invention has been presented for purposes of illustration and description, and is not intended to be exhaustive or limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the various embodiments described.

Claims (6)

1.一种井震信息融合方法,其特征在于,包括:1. A method for fusing well and seismic information, comprising: 获取声波时差测井数据和偏移剖面;Acquire sonic time-of-day logging data and migration profiles; 对所述声波时差测井数据进行中值滤波、高斯滤波处理;performing median filtering and Gaussian filtering on the acoustic time difference logging data; 对所述偏移剖面提取结构张量和相干体,利用结构张量和相干体计算结构张量场;extracting a structure tensor and a coherence volume from the offset section, and calculating a structure tensor field using the structure tensor and the coherence volume; 利用混合邻域法,结合结构张量场对测井数据进行插值处理,获取插值速度模型和扩散时间场;The hybrid neighborhood method is used to interpolate the well logging data in combination with the structure tensor field to obtain the interpolation velocity model and diffusion time field. 利用所述扩散时间场计算出加权系数,融合层析速度模型和插值速度模型,获取融合速度模型;Calculating a weighting coefficient using the diffusion time field, fusing a tomographic velocity model and an interpolation velocity model to obtain a fused velocity model; 其中,所述结构张量场的计算公式如下:The calculation formula of the structure tensor field is as follows: , 其中, S为偏移剖面中提取的结构张量,c为偏移剖面中提取的相干体,D为结构张量场;Where S is the structure tensor extracted from the migration section, c is the coherence volume extracted from the migration section, and D is the structure tensor field; 其中,所述融合速度的计算公式如下:The calculation formula of the fusion speed is as follows: , 其中, Tmax为时间阈值,τ为扩散时间场,l为加权系数,vinterp为插值速度模型,vtomo为层析速度模型,vmix为融合速度。Where T max is the time threshold, τ is the diffusion time field, l is the weighting coefficient, v interp is the interpolation velocity model, v tomo is the tomographic velocity model, and v mix is the fusion velocity. 2.根据权利要求1所述的井震信息融合方法,其特征在于,所述插值公式如下:2. The well-seismic information fusion method according to claim 1, wherein the interpolation formula is as follows: , 其中,τ为扩散时间场,D为结构张量场,χ表示声波时差测井位置区域,p为扩散时间最短的声波时差测井值,q为插值结果。Where τ is the diffusion time field, D is the structure tensor field, χ represents the sonic transit time logging location area, p is the sonic transit time logging value with the shortest diffusion time, and q is the interpolation result. 3.一种井震信息融合装置,其特征在于,3. A well-seismic information fusion device, characterized in that: 获取模块,获取声波时差测井数据和偏移剖面;Acquisition module, which acquires acoustic time difference logging data and migration profiles; 处理模块,对所述声波时差测井数据进行中值滤波、高斯滤波处理;A processing module, performing median filtering and Gaussian filtering on the acoustic time difference logging data; 计算模块,对所述偏移剖面提取结构张量和相干体,利用结构张量和相干体计算结构张量场;a calculation module, extracting a structure tensor and a coherence volume from the offset section, and calculating a structure tensor field using the structure tensor and the coherence volume; 插值处理模块,利用混合邻域法,结合结构张量场对测井数据进行插值处理,获取插值速度模型和扩散时间场;The interpolation processing module uses the hybrid neighborhood method and the structural tensor field to interpolate the logging data to obtain the interpolation velocity model and diffusion time field; 融合模块,利用所述扩散时间场计算出加权系数,融合层析速度模型和插值速度模型,获取融合速度模型;A fusion module calculates a weighting coefficient using the diffusion time field, fuses the tomographic velocity model and the interpolation velocity model, and obtains a fused velocity model; 其中,所述结构张量场的计算公式如下:The calculation formula of the structure tensor field is as follows: , 其中, S为偏移剖面中提取的结构张量,c为偏移剖面中提取的相干体,D为结构张量场;Where S is the structure tensor extracted from the migration section, c is the coherence volume extracted from the migration section, and D is the structure tensor field; 其中,所述融合速度的计算公式如下:The calculation formula of the fusion speed is as follows: , 其中,Tmax为时间阈值,τ为扩散时间场,l为加权系数,vinterp为插值速度模型,vtomo为层析速度模型,vmix为融合速度。Where T max is the time threshold, τ is the diffusion time field, l is the weighting coefficient, v interp is the interpolation velocity model, v tomo is the tomographic velocity model, and v mix is the fusion velocity. 4.根据权利要求3所述的井震信息融合装置,其特征在于,所述插值公式如下:4. The well-seismic information fusion device according to claim 3, wherein the interpolation formula is as follows: , 其中,τ为扩散时间场,D为结构张量场,χ表示声波时差测井位置区域,p为扩散时间最短的声波时差测井值,q为插值结果。Where τ is the diffusion time field, D is the structure tensor field, χ represents the sonic transit time logging location area, p is the sonic transit time logging value with the shortest diffusion time, and q is the interpolation result. 5.一种电子设备,其特征在于,所述电子设备包括:5. An electronic device, characterized in that the electronic device comprises: 至少一个处理器;以及,at least one processor; and, 与所述至少一个处理器通信连接的存储器;其中,a memory communicatively connected to the at least one processor; wherein, 所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行权利要求1或2所述的井震信息融合方法。The memory stores instructions that can be executed by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute the well-seismic information fusion method according to claim 1 or 2. 6.一种非暂态计算机可读存储介质,其特征在于,所述非暂态计算机可读存储介质存储计算机指令,该计算机指令用于使计算机执行权利要求1或2所述的井震信息融合方法。6. A non-transitory computer-readable storage medium, characterized in that the non-transitory computer-readable storage medium stores computer instructions, which are used to enable a computer to execute the well-seismic information fusion method according to claim 1 or 2.
CN202111210323.XA 2021-10-18 2021-10-18 A method, device, electronic device and medium for fusion of well-seismic information Active CN115994338B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111210323.XA CN115994338B (en) 2021-10-18 2021-10-18 A method, device, electronic device and medium for fusion of well-seismic information

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111210323.XA CN115994338B (en) 2021-10-18 2021-10-18 A method, device, electronic device and medium for fusion of well-seismic information

Publications (2)

Publication Number Publication Date
CN115994338A CN115994338A (en) 2023-04-21
CN115994338B true CN115994338B (en) 2025-09-05

Family

ID=85990661

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111210323.XA Active CN115994338B (en) 2021-10-18 2021-10-18 A method, device, electronic device and medium for fusion of well-seismic information

Country Status (1)

Country Link
CN (1) CN115994338B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200115A (en) * 2014-09-12 2014-12-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Geostatistics simulation based full-formation velocity modeling method
CN112017289A (en) * 2020-08-31 2020-12-01 电子科技大学 Well-seismic combined initial lithology model construction method based on deep learning

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5729451A (en) * 1995-12-01 1998-03-17 Coleman Research Corporation Apparatus and method for fusing diverse data
US8131469B2 (en) * 2008-10-03 2012-03-06 Baker Hughes Incorporated Data acquisition and processing for invasion profile and gas zone analysis with NMR dual or multiple interecho spacing time logs
CN106703787B (en) * 2017-02-27 2019-11-29 中国石油大学(北京) A kind of the well track parameter calculation method and device of radially horizontal well
CN109884700B (en) * 2019-03-20 2021-02-26 中国石油化工股份有限公司 Multi-information fusion seismic velocity modeling method
CN110333551B (en) * 2019-07-26 2020-09-25 长江大学 Dolostone reservoir prediction method and system based on well-seismic combination and storage medium
CN111502647B (en) * 2020-03-27 2021-02-02 中国石油化工股份有限公司石油工程技术研究院 Method and device for determining drilling geological environment factors and storage medium
EP4133310B1 (en) * 2020-04-07 2025-11-12 Chevron U.S.A. Inc. System and method for seismic inversion
CN112433247B (en) * 2020-11-17 2022-09-02 中国石油化工股份有限公司 While-drilling adjusting method and device for position of stratum to be drilled

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200115A (en) * 2014-09-12 2014-12-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Geostatistics simulation based full-formation velocity modeling method
CN112017289A (en) * 2020-08-31 2020-12-01 电子科技大学 Well-seismic combined initial lithology model construction method based on deep learning

Also Published As

Publication number Publication date
CN115994338A (en) 2023-04-21

Similar Documents

Publication Publication Date Title
Kontoroupi et al. Online noise identification for joint state and parameter estimation of nonlinear systems
CN111596366A (en) Wave impedance inversion method based on seismic signal optimization processing
US9588245B2 (en) Efficient wavefield extrapolation in anisotropic media
CN113820745A (en) Seismic velocity modeling method, device, electronic apparatus, and medium
CN118210024A (en) Seismic structure interpretation method and device based on multi-source multi-task deep learning
CN107479091A (en) A kind of method for extracting reverse-time migration angle gathers
CN115994338B (en) A method, device, electronic device and medium for fusion of well-seismic information
CN111502647B (en) Method and device for determining drilling geological environment factors and storage medium
CN112363243A (en) Method, device, computer equipment and medium for predicting high-quality reservoir of shale oil
EP3510423A1 (en) Methods for efficient wavefield solutions
CN112308206B (en) Magnetotelluric-seismic wave first arrival time joint inversion method based on deep learning
CN115993645B (en) Migration velocity modeling method, device, electronic device and medium
CN112630840A (en) Random inversion method based on statistical characteristic parameters and processor
CN115166822A (en) Method, device and electronic device for predicting reservoir elastic parameters
CN119375936B (en) Stratum dip angle determining method and device based on local linear constraint
KR101969718B1 (en) Peridynamic method having added mirroring node and numerical analysis apparatus using the same
CN119667776B (en) Diffraction wave impedance inversion method and device, electronic equipment and medium
CN116859452A (en) VTI medium ray tracing method and device, electronic equipment and medium
CN113887034A (en) Electromagnetic wave logging while drilling cross-dimensional simulation method suitable for 2D formation model
CN120122167A (en) Method, device and medium for constructing low-frequency background model guided by seismic data structure
CN113219534B (en) Pre-stack depth migration speed quality control method and device, medium and electronic equipment
US9971050B2 (en) Generalized internal multiple imaging
CN116009081B (en) Full waveform inversion method, device, electronic equipment and medium based on error function
CN118191929A (en) Seismic velocity model boundary fusion calculation method and device, electronic equipment and medium
CN119620208A (en) Method, equipment and medium for interpreting geological radar intelligent data by integrating multi-source data

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