Far detection sound wave three-dimensional imaging method, system, equipment and medium

Abstract

The invention provides a far detection sound wave three-dimensional imaging method, system, equipment and medium, and the method comprises the following steps: converting the longitudinal section data of array sound wave logging at different angles into gray data in sequence, and enabling the logging data at different angles to form a cylinder in sequence in a coaxial manner, each piece of cross section data is stored layer by layer in the vertical direction; performing data interpolation on each layer of cross section according to a circular interpolation method to obtain logging cross section data with higher resolution; converting each piece of cross section gray scale data obtained through interpolation into a gray scale image; reading each layer of cross section grey-scale image in the grey-scale image for three-dimensional imaging to obtain a three-dimensional structure; the opacity and color of the three-dimensional structure are adjusted, and three-dimensional imaging of the far-detection sound waves is completed; according to the method, the sound wave signal data received by the limited array sound wave instrument is fully utilized for interpolation, the real condition of the well body is completely restored under the support of the limited data, and the method has a wide practical application scene.

Description

Technical Field

[0001] The invention belongs to the field of borehole geophysics and relates to a remote detection acoustic wave three-dimensional imaging method, system, equipment and medium. Background Art

[0002] The essence of acoustic logging data visualization is to draw the acoustic signal received by the sampler into a graph in a certain way, so that geological professionals can view and analyze the geological structure in the form of graphics and images, and provide guidance for exploration projects. The application of visualization technology to the processing and interpretation of logging data is called logging data visualization technology. Well logging data visualization involves multiple disciplines such as geophysics, oil and gas exploration, and computer graphics, and is the mutual combination and cross-application of these disciplines. The essence of logging data visualization is to convert abstract and difficult-to-understand array acoustic logging data into specific, intuitive and easy-to-understand graphic image information, thereby helping researchers and users to analyze. At the same time, the resolution of array acoustic logging data can be further improved through the circular interpolation technology of array acoustic logging data. In recent years, due to the development and popularization of computer imaging technology and tools, the logging field has gradually begun to use visualization technology to assist related exploration work, such as the use of computer imaging technology for borehole imaging and reflector imaging, which has provided great help for exploration work.

[0003] Array acoustic logging is an image constructed from the reflected wave data received by the acoustic logging instrument, which represents the longitudinal section of the wellbore. In the past, the visualization method of logging data was usually based on array acoustic logging, and a series of processing and analysis such as denoising and recognition were performed on the image to assist researchers and users; however, two-dimensional images are difficult to reflect the distribution characteristics of reflectors in three-dimensional space, and array acoustic logging images only contain a few angles and cannot restore the condition of the entire well body, which to a certain extent limits the analysis and research of experts; to make up for the above shortcomings, the prior art proposes a three-dimensional imaging method for acoustic logging data, such as Figure 2 As shown in the figure, the existing 3D imaging is to construct a 3D coaxial body based on the longitudinal section image of acoustic logging. Although it reflects the 3D structure and actual condition of the well body to a certain extent, the data volume is small and the detailed information is lacking, and it still cannot completely solve the above problems. Therefore, 3D imaging based on array acoustic logging data has become a difficult problem that needs to be solved urgently. Summary of the invention

[0004] In view of the problems existing in the prior art, the present invention provides a remote detection acoustic three-dimensional imaging method, system, equipment and medium, which improves the existing acoustic logging data visualization method, can improve the resolution of logging data and increase the interactive display of visualization imaging, and restore the actual situation of the well body to the greatest extent.

[0005] The present invention is achieved through the following technical solutions:

[0006] A remote detection acoustic wave three-dimensional imaging method comprises the following steps:

[0007] The longitudinal section data of array acoustic logging at different angles are sequentially converted into grayscale data, the logging data at different angles are sequentially formed into a cylinder in a coaxial manner, and each cross-section data is stored by layer in the vertical direction;

[0008] Interpolating the data of each cross section according to the circular interpolation method to obtain higher resolution logging cross section data;

[0009] Convert the grayscale data of each cross section obtained through interpolation into a grayscale image;

[0010] Read the grayscale images of the cross sections of each layer in the grayscale image to perform three-dimensional imaging and obtain a three-dimensional structure;

[0011] The opacity and color of the three-dimensional structure are adjusted to complete the three-dimensional imaging of the far-reaching sound waves.

[0012] Preferably, the circular interpolation method is to obtain interpolation data between two data by fusing linear interpolation and applying the linear interpolation principle to arc interpolation to obtain cross-sectional data with higher resolution.

[0013] Preferably, the circular interpolation method is:

[0014]

[0015]

[0016]

[0017] In the formula, Q 11 , Q 12 , Q 21 With Q 22 They represent four points on two adjacent radii where data exists on the cross section, where Q 11 With Q 12 , Q 21 With Q 22 are on the same radius, and Q 11 With Q 21 , Q 12 With Q 22 The distances from the center of the cross section are the same, R1 and R2 represent Q 11 With Q 12 , Q 21 With Q 22 The midpoint of the arc In the previous point, F represents the data value at a certain point.

[0018] Preferably, the process of converting the grayscale data of each cross section obtained by interpolation into a grayscale image is:

[0019]

[0020] In the formula, data i,j represents the well logging data value obtained by interpolation at the midpoint [i, j] of the cross section, max(data) represents the maximum value of the acoustic logging data value of the entire wellbore after interpolation, gamma_data i,j Represents the ratio of the acoustic logging data value at point [i, j] to the maximum value, used for the grayscale value in the grayscale image.

[0021] Preferably, the process of reading the grayscale images of each layer of cross section in the grayscale image to perform three-dimensional imaging and obtain the three-dimensional structure is: using the idea of ​​ray casting volume rendering, and drawing each layer of grayscale image into the three-dimensional well body structure through the VTK tool.

[0022] Preferably, the ray casting is a direct numerical calculation method of volume rendering integral, which calculates the propagation of light by integrating the interaction of light along the observation light, and can map scalar data to optical properties for simulating original volume data. The specific formula is as follows:

[0023]

[0024]

[0025] Among them, C i and A i represent the color value and opacity sampled on the i-th body texture, respectively. The color value and opacity sampled on the body texture are the data contained in the voxel. and They represent the accumulated color value and opacity respectively. After the volume rendering is completed, the 3D imaging angle and size are adjusted synchronously.

[0026] Preferably, the process of adjusting the opacity and color of the three-dimensional structure is: adjusting the opacity and color of the three-dimensional imaging by adding an opacity transfer function and a color transfer function.

[0027] A remote detection acoustic wave three-dimensional imaging system, comprising:

[0028] A preprocessing module is used to sequentially convert longitudinal section data of array acoustic logging at different angles into grayscale data, sequentially form cylinders with logging data at different angles in a coaxial manner, and store each cross-section data by layer in the vertical direction;

[0029] A processing module, used for performing data interpolation on the cross sections of each layer according to a circular interpolation method to obtain logging cross section data with higher resolution;

[0030] A conversion module, used for converting the grayscale data of each cross section obtained through interpolation into a grayscale image;

[0031] A three-dimensional imaging module is used to read the grayscale images of the cross sections of each layer in the grayscale image to perform three-dimensional imaging and obtain a three-dimensional structure;

[0032] The output module is used to adjust the opacity and color of the three-dimensional structure to complete the three-dimensional imaging of the far-detection sound waves.

[0033] A computer device comprises a memory, a processor and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, the steps of a remote detection acoustic wave three-dimensional imaging method are implemented.

[0034] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the steps of a remote detection acoustic wave three-dimensional imaging method.

[0035] Compared with the prior art, the present invention has the following beneficial technical effects:

[0036] The present invention provides a remote detection acoustic wave three-dimensional imaging method, system, equipment and medium, comprising the following steps: sequentially converting longitudinal section data of acoustic wave well logging at different angles into grayscale data, sequentially forming a cylinder with the well logging data at different angles in a coaxial manner, and storing each cross-sectional data by layer in the vertical direction; performing data interpolation on the cross sections of each layer according to a circular interpolation method to obtain well logging cross-sectional data with higher resolution; converting the grayscale data of each cross section obtained by interpolation into a grayscale image; reading the grayscale image of each cross section in the grayscale image to perform three-dimensional imaging to obtain a three-dimensional structure; and performing three-dimensional imaging on the three-dimensional structure. The opacity and color are adjusted to complete the three-dimensional imaging of the long-range detection sound waves; the present application makes full use of the limited sound wave signal data received by the array acoustic wave instrument for interpolation to obtain data with higher resolution. With the support of limited data, the real condition of the well body is completely restored, and it has a wide range of practical application scenarios; compared with the previous two-dimensional plane imaging and three-dimensional section coaxial body imaging, the three-dimensional imaging based on the array acoustic wave logging data of the present application has clearer information on details, higher overall resolution, and is closer to the real well body condition, which greatly reduces the difficulty of logging data analysis and improves the accuracy of analysis. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] Figure 1 This is a flow chart of a remote detection acoustic wave three-dimensional imaging method in an embodiment of the present invention;

[0038] Figure 2 A comparison diagram of a three-dimensional volume rendering result and a three-dimensional coaxial volume result in the prior art;

[0039] Figure 3 This is a flowchart of a remote detection acoustic wave three-dimensional imaging method in an embodiment of the present invention;

[0040] Figure 4 A schematic diagram of converting longitudinal section data of acoustic logging into cross-sectional data in an embodiment of the present invention;

[0041] Figure 5 is a schematic diagram of a circular interpolation algorithm according to an embodiment of the present invention;

[0042] Figure 6 Schematic diagram of a volume rendering method based on the ray casting principle in an embodiment of the present invention. DETAILED DESCRIPTION

[0043] The present invention is further described in detail below in conjunction with specific embodiments, which are intended to explain the present invention rather than to limit it.

[0044] In order to enable those skilled in the art to better understand the scheme of the present invention, the technical scheme in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

[0045] It should be noted that the terms "first", "second", etc. in the specification and claims of the present invention and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units that are clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices.

[0046] The present invention provides a remote detection acoustic wave three-dimensional imaging method, such as Figure 1 As shown, the following steps are included:

[0047] The longitudinal section data of array acoustic logging at different angles are sequentially converted into grayscale data, the logging data at different angles are sequentially formed into a cylinder in a coaxial manner, and each cross-section data is stored by layer in the vertical direction;

[0048] Interpolating the data of each cross section according to the circular interpolation method to obtain higher resolution logging cross section data;

[0049] Convert the grayscale data of each cross section obtained through interpolation into a grayscale image;

[0050] Read the grayscale images of the cross sections of each layer in the grayscale image to perform three-dimensional imaging and obtain a three-dimensional structure;

[0051] The opacity and color of the three-dimensional structure are adjusted to complete the three-dimensional imaging of the far-reaching sound waves.

[0052] Preferably, the circular interpolation method is to obtain interpolation data between two data by fusing linear interpolation and applying the linear interpolation principle to arc interpolation, so as to obtain cross-sectional data with higher resolution; further, the circular interpolation method is:

[0053]

[0054]

[0055]

[0056] In the formula, Q 11 , Q 12 , Q 21 With Q 22 They represent four points on two adjacent radii where data exists on the cross section, where Q 11 With Q 12 , Q 21 With Q 22 are on the same radius, and Q 11 With Q 21 , Q 12 With Q 22 The distances from the center of the cross section are the same, R1 and R2 represent Q 11 With Q 12 , Q 21 With Q 22 The midpoint of the arc In the previous point, F represents the data value at a certain point.

[0057] Preferably, the process of converting the grayscale data of each cross section obtained by interpolation into a grayscale image is:

[0058]

[0059] In the formula, data i,j represents the well logging data value obtained by interpolation at the midpoint [i, j] of the cross section, max(data) represents the maximum value of the acoustic logging data value of the entire wellbore after interpolation, gamma_data i,j Represents the ratio of the acoustic logging data value at point [i, j] to the maximum value, used for the grayscale value in the grayscale image.

[0060] Preferably, the process of reading the grayscale images of each layer of cross section in the grayscale image to perform three-dimensional imaging and obtain the three-dimensional structure is: using the idea of ​​ray casting volume rendering, and drawing each layer of grayscale image into the three-dimensional well body structure through the VTK tool.

[0061] It should be noted that the idea of ​​ray casting volume rendering is mainly to emit a ray of light from each pixel of the image along a fixed direction, usually the line of sight, and the ray of light passes through the entire image sequence. During the passage, the image sequence is sampled to obtain color information; at the same time, the color value is accumulated according to the light absorption model until the light passes through the entire image sequence. The final color value is the color of the rendered image. This idea is similar to the process in which the human eye projects the light received in the three-dimensional space onto the retina to form a two-dimensional image.

[0062] It should be further explained that VTK (Visualization Toolkit) is an open source software system for 3D computer graphics, image processing and visualization. VTK is designed and implemented based on object-oriented principles. Its kernel is built in C++ and contains about 250,000 lines of code, more than 2,000 classes, and several conversion interfaces. Therefore, VTK can also be used freely in various languages ​​such as Java, Tcl / Tk and Python.

[0063] Preferably, the ray casting is a direct numerical calculation method of volume rendering integral, which calculates the propagation of light by integrating the interaction of light along the observation light, and can map scalar data to optical properties for simulating original volume data. The specific formula is as follows:

[0064]

[0065]

[0066] Among them, C i and A i represent the color value and opacity sampled on the i-th body texture, respectively. The color value and opacity sampled on the body texture are the data contained in the voxel. and They represent the accumulated color value and opacity respectively. After the volume rendering is completed, the 3D imaging angle and size are adjusted synchronously.

[0067] Preferably, the process of adjusting the opacity and color of the three-dimensional structure is: adjusting the opacity and color of the three-dimensional imaging by adding an opacity transfer function and a color transfer function.

[0068] It should be noted that the opacity transfer function is a concept in VTK, which defines the mapping relationship of the opacity of pixels inside an object. Specifically, the opacity transfer function can map the grayscale value of the sampling point in the light projection process to different opacity values ​​to determine the final color value.

[0069] The opacity transfer function in VTK can usually be implemented through the vtkPiecewiseFunction class. This class is a piecewise linear scalar mapping function, whose piecewise linear units can be added, modified, and deleted to meet different needs; if you need to set the opacity transfer function, you first need to create a vtkPiecewiseFunction object, and then add each segment by calling the AddSegment() function; four segments are added through the AddSegment() function, each segment consists of two points, the first point is the starting value of the segment, that is, the gray value corresponding to the light sampling point, and the second point is the opacity value of the segment. The last statement applies the opacity transfer function to the volumeProperty object to determine the opacity of its internal pixels.

[0070] It should be further explained that the color transfer function is a concept in VTK, which is used to map a scalar value to a color value. This color value can be an RGB value or an HSV value. The color transfer function in VTK can usually be implemented through the vtkColorTransferFunction class, which is a scalar mapping function. Its mapping relationship can be added, modified, and deleted to meet different needs. If you need to set a color transfer function, you first need to create a vtkColorTransferFunction object, and then add each breakpoint by calling the AddRGBPoint() or AddHSVPoint() function. For example, color rFunc is a color transfer function. Five breakpoints are added through the AddRGBPoint() function. Each breakpoint consists of four points. The first point is the starting value of the segment, that is, the grayscale value corresponding to the light sampling point. The next three points are the values ​​of the three RGB color channels respectively. The last statement applies the color transfer function to the volumeProperty object to determine the color of its internal pixels. If the components of the multivariate data are independent of each other, the corresponding opacity, gradient opacity, and color transfer function can be set for each component separately. It should be noted that the setting of the opacity transfer function is closely related to the specific algorithm and scenario of volume rendering. In practical applications, it needs to be set according to specific needs.

[0071] Example

[0072] A method for long-range detection three-dimensional imaging, the flow chart of the method for long-range detection three-dimensional imaging of the present invention is as follows Figure 3 As shown, the method comprises:

[0073] S1, data preprocessing, sequentially converting the longitudinal section data of array acoustic logging at different angles into grayscale data, sequentially forming cylinders with the logging data at different angles in a coaxial manner, and storing each cross-sectional data layer by layer in the vertical direction. The process of obtaining cross-sectional data by this method is as follows Figure 4 As shown in the figure, the left figure shows that the array acoustic wave data is initially in the form of a longitudinal section of the well body, the middle image shows the coaxial three-dimensional structure reconstructed according to the image sequence and interval angle, and the right figure shows that the constructed three-dimensional structure is layered in the vertical dimension to obtain the cross-sectional data of the array acoustic wave logging. At this time, the cross-sectional data is very sparse, and the data points only exist at the intersection of the solid lines, which cannot completely restore the actual well body conditions.

[0074] S2, data interpolation, according to the innovative circular interpolation algorithm, the data of each layer cross-section data obtained in S1 is interpolated to obtain higher resolution logging cross-section data. The circular interpolation algorithm is designed based on the linear interpolation algorithm, and its specific formula is as follows:

[0075]

[0076]

[0077]

[0078] In the above formula, Q 11 , Q 12 , Q 21 With Q 22 They represent four points on two adjacent radii where data exists on the cross section, where Q 11 With Q 12 , Q 21 With Q 22 are on the same radius, and Q 11 With Q 21 , Q 12 With Q 22 The distances from the center of the cross section are the same, R1 and R2 represent Q 11 With Q 12 , Q 21 With Q 22 The midpoint of the arc In the previous point, F represents the data value at a certain point. By fusing linear interpolation and applying the linear interpolation principle to arc interpolation, the interpolation data between the two data is calculated, thereby obtaining cross-sectional data with higher resolution.

[0079] The core of circular interpolation is to use formulas (1), (2) and (3) to interpolate curves between existing data points according to the ratio of the line segment length or arc length between the interpolation point and the original data point. The circular interpolation method is mainly based on the arc length ratio interpolation, which is essentially the same as the linear interpolation idea. Figure 5 As shown in the figure, the left figure is the original cross-sectional data without interpolation. Points R1 and R2 are interpolated by formula (1), and then point P is interpolated by formula (3). By repeating the above process, the data is converted from the left figure to the form of the intermediate image. At this time, the data has been increased by a circle and a diameter. The right figure is a schematic diagram of the result of repeating the interpolation process until the data resolution is sufficient to restore the entire cross section.

[0080] S3, image conversion, maps the grayscale data of each layer of cross-section obtained by interpolation in S2 with the pixel points in the image, and uses the data value as the grayscale value of the pixel point, thereby converting the multi-layer cross-section data into the same number of cross-section grayscale images. The specific implementation is shown in the following formula:

[0081]

[0082] In the formula, data i,j represents the well logging data value obtained by interpolation at the midpoint [i, j] of the cross section, and max(data) represents the maximum value of the acoustic logging data value of the entire wellbore after interpolation.

[0083] gamma_data i,j Represents the ratio of the acoustic logging data value at point [i, j] to the maximum value, used for the grayscale value in the grayscale image.

[0084] S4, 3D volume rendering. Read the grayscale images of each cross-section in S3 to perform 3D imaging, and use the mouse and keyboard to rotate and scale the constructed 3D structure. Figure 6 As shown in the figure, the 3D imaging process adopts the idea of ​​ray casting volume rendering, which is essentially to send a line of sight to the object through the viewpoint, map the optical properties with the scalar based on the light influence received by the line of sight during the transmission process, and then draw the result on the plane used for 3D imaging, so as to achieve the 3D effect on the plane. Based on the principle of ray casting, the grayscale images of each layer are drawn into the 3D well structure through the VTK tool. The 3D imaging obtained by drawing can greatly restore the details of the well condition. The excellent light processing ability enables the 3D imaging to better restore the actual well by adjusting the opacity and color. At the same time, the corresponding keyboard and mouse response functions are used to transmit its actions to the 3D imaging process in real time, thereby triggering the corresponding bound functions.

[0085] S5, adjust the three-dimensional imaging, adjust the opacity and color of the three-dimensional structure through the opacity control and color control, so that the three-dimensional imaging results conform to the actual well body conditions or meet specific needs. The same principle as the keyboard and mouse, this process is achieved by adding a transparency and color adjustment module to the drawing method. The transparency and color adjustment module will reflect the transparency and color parameters of the current three-dimensional imaging in real time in a visual form, and can be manually adjusted to pass the adjusted parameters into the volume rendering process, thereby redrawing the three-dimensional imaging.

[0086] In summary, the long-range detection three-dimensional imaging method of the present invention is suitable for array acoustic wave data analysis in the field of well logging exploration. It does not require the provision of a large amount of well logging data. The sparse data is complemented by a circular interpolation algorithm, and three-dimensional imaging is performed by a ray casting volume rendering method, which completely restores the true condition of the well body. It can reduce the difficulty of well logging workers in analyzing acoustic wave array data, improve the precision and accuracy of the analysis, and has great application value.

[0087] The present invention provides a remote detection acoustic wave three-dimensional imaging system, comprising:

[0088] A preprocessing module is used to sequentially convert longitudinal section data of array acoustic logging at different angles into grayscale data, sequentially form cylinders with logging data at different angles in a coaxial manner, and store each cross-section data by layer in the vertical direction;

[0089] A processing module, used for performing data interpolation on the cross sections of each layer according to a circular interpolation method to obtain logging cross section data with higher resolution;

[0090] A conversion module, used for converting the grayscale data of each cross section obtained through interpolation into a grayscale image;

[0091] A three-dimensional imaging module is used to read the grayscale images of the cross sections of each layer in the grayscale image to perform three-dimensional imaging and obtain a three-dimensional structure;

[0092] The output module is used to adjust the opacity and color of the three-dimensional structure to complete the three-dimensional imaging of the far-detection sound waves.

[0093] In another embodiment of the present invention, a computer device is provided, the computer device including a processor and a memory, the memory is used to store a computer program, the computer program includes program instructions, and the processor is used to execute the program instructions stored in the computer storage medium. The processor may be a central processing unit (CPU), or other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field-programmable gate arrays (FPGA) or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc., which are the computing core and control core of the terminal, which are suitable for implementing one or more instructions, and are specifically suitable for loading and executing one or more instructions in the computer storage medium to implement the corresponding method flow or corresponding function; the processor described in the embodiment of the present invention can be used for the operation of a remote detection acoustic wave three-dimensional imaging method.

[0094] In another embodiment of the present invention, the present invention also provides a storage medium, specifically a computer-readable storage medium (Memory), which is a memory device in a computer device for storing programs and data. It can be understood that the computer-readable storage medium here can include both built-in storage media in the computer device and, of course, extended storage media supported by the computer device. The computer-readable storage medium provides a storage space, which stores the operating system of the terminal. In addition, one or more instructions suitable for being loaded and executed by the processor are also stored in the storage space, and these instructions can be one or more computer programs (including program codes). It should be noted that the computer-readable storage medium here can be a high-speed RAM memory or a non-volatile memory, such as at least one disk memory. The processor can load and execute one or more instructions stored in the computer-readable storage medium to implement the corresponding steps of a remote detection acoustic wave three-dimensional imaging method in the above embodiment.

[0095] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0096] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0097] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0098] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0099] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of the present invention.

Claims

1. A remote detection acoustic wave three-dimensional imaging method, characterized in that: The following steps are involved: The longitudinal section data of array acoustic logging at different angles are sequentially converted into grayscale data, the logging data at different angles are sequentially formed into a cylinder in a coaxial manner, and each cross-section data is stored by layer in the vertical direction; Interpolating the data of each cross section according to the circular interpolation method to obtain higher resolution logging cross section data; Convert the grayscale data of each cross section obtained through interpolation into a grayscale image; Read the grayscale images of the cross sections of each layer in the grayscale image to perform three-dimensional imaging and obtain a three-dimensional structure; The opacity and color of the three-dimensional structure are adjusted to complete the three-dimensional imaging of the far-reaching sound waves.

2. The method for three-dimensional imaging of remote sound waves according to claim 1, characterized in that: The circular interpolation method is to obtain interpolation data between two data by fusing linear interpolation and applying the linear interpolation principle to arc interpolation, thereby obtaining cross-sectional data with higher resolution.

3. The method for three-dimensional imaging of remote acoustic waves according to claim 1, characterized in that: The circular interpolation method is: In the formula, Q 11 , Q 12 , Q 21 With Q 22 They represent four points on two adjacent radii where data exists on the cross section, where Q 11 With Q 12 , Q 21 With Q 22 are on the same radius, and Q 11 With Q 21 , Q 12 With Q 22 The distances from the center of the cross section are the same, R1 and R2 represent Q 11 With Q 12 , Q 21 With Q 22 The midpoint of the arc In the previous point, F represents the data value at a certain point.

4. The method for three-dimensional imaging of remote sound waves according to claim 1, characterized in that: The process of converting the grayscale data of each cross section obtained by interpolation into a grayscale image is as follows: In the formula, data i,j represents the well logging data value obtained by interpolation at the midpoint [i, j] of the cross section, max(data) represents the maximum value of the acoustic logging data value of the entire wellbore after interpolation, gamma_data i,j Represents the ratio of the acoustic logging data value at point [i, j] to the maximum value, used for the grayscale value in the grayscale image.

5. The method for three-dimensional imaging of remote sound waves according to claim 1, characterized in that: The process of reading the grayscale images of each layer of the cross section in the grayscale image to perform three-dimensional imaging and obtain the three-dimensional structure is: using the idea of ​​ray casting volume rendering, and drawing each layer of grayscale image into the three-dimensional well body structure through the VTK tool.

6. A remote detection acoustic wave three-dimensional imaging method according to claim 5, characterized in that: The ray casting is a direct numerical calculation method for volume rendering integrals. It calculates the propagation of light by integrating the interaction of light along the observation light, and can map scalar data to optical properties to simulate the original volume data. The specific formula is as follows: Among them, C i and A i represent the color value and opacity sampled on the i-th body texture, respectively. The color value and opacity sampled on the body texture are the data contained in the voxel. and They represent the accumulated color value and opacity respectively. After the volume rendering is completed, the 3D imaging angle and size are adjusted synchronously.

7. The method for three-dimensional imaging of remote sound waves according to claim 1, characterized in that: The process of adjusting the opacity and color of the three-dimensional structure is: adjusting the opacity and color of the three-dimensional imaging by adding an opacity transfer function and a color transfer function.

8. A remote detection acoustic wave three-dimensional imaging system, characterized in that: A remote detection acoustic wave three-dimensional imaging method according to any one of claims 1 to 7, comprising: A preprocessing module is used to sequentially convert longitudinal section data of array acoustic logging at different angles into grayscale data, sequentially form cylinders with logging data at different angles in a coaxial manner, and store each cross-section data by layer in the vertical direction; A processing module, used for performing data interpolation on the cross sections of each layer according to a circular interpolation method to obtain logging cross section data with higher resolution; A conversion module, used for converting the grayscale data of each cross section obtained through interpolation into a grayscale image; A three-dimensional imaging module is used to read the grayscale images of the cross sections of each layer in the grayscale image to perform three-dimensional imaging and obtain a three-dimensional structure; The output module is used to adjust the opacity and color of the three-dimensional structure to complete the three-dimensional imaging of the far-detection sound waves.

9. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the steps of the remote detection acoustic wave three-dimensional imaging method as claimed in any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium storing a computer program, characterized in that: When the computer program is executed by a processor, the steps of a remote detection acoustic wave three-dimensional imaging method as claimed in any one of claims 1 to 7 are implemented.

Images