Near-surface and mid-deep layer Q field fusion method, system and device based on mask body and medium

Abstract

The invention belongs to the technical field of oil exploration information processing, and particularly discloses a near-surface and mid-deep layer Q field fusion method, system, equipment and medium based on a mask body, and the method comprises the steps: generating a near-surface Q field F1 and a mid-deep layer Q field F2, carrying out the weighted subtraction of the two Q fields to form a data body F3, and carrying out the fusion of the near-surface Q field F1 and the mid-deep layer Q field F2; carrying out abnormal value removal processing on the data body F3 to obtain a data body F3 '', carrying out mask transformation on the data body F3'' to respectively obtain mask bodies F1 / 0 and F0 / 1, and taking a data boundary where 0 and 1 values in the mask bodies F1 / 0 and F0 / 1 are located as an ideal fusion surface; and multiplying the points corresponding to the near-surface Q field F1 and the middle-deep layer Q field F2 by the mask bodies F1 / 0 and F0 / 1 to obtain a data body F < upper > and a data body F < lower > which are complementary along the ideal fusion surface, and finally adding the data body F < upper > and the data body F < lower >. The near-surface Q field and the middle-deep layer Q field can be intelligently and efficiently fused, the fusion precision is high, and the efficiency is high. The method is suitable for near-surface and middle-deep layer Q field fusion and full-series Q field establishment.

Description

Technical Field

[0001] The present invention relates to the technical field of oil exploration information processing, and in particular to a mask-based near-surface and mid-deep Q-field fusion method, system, equipment and medium. Background Art

[0002] In the field of oil exploration, the quality factor Q is an important parameter that describes the absorption and attenuation characteristics of seismic waves by the medium. A reasonable and detailed description of the Q field is crucial to improving the resolution of seismic data. At present, the Q processing technology series based on "Q modeling, Q compensation, and Q migration" has been widely used and has achieved good results, especially in basins with complex near-surface conditions, such as the Ordos Basin. Due to the serious absorption and attenuation of data by its thick loess plateau, the main frequency of the data is low and the frequency band is narrow, which seriously restricts the prediction and lithology interpretation of thin-layer sand bodies. The application of Q processing technology can better solve the high-frequency loss, phase dispersion and wavelet distortion caused by complex near-surface absorption and attenuation in the basin.

[0003] The near-surface Q field and mid-deep Q field of Q processing technology are obtained by calculating and fitting micro-logging data and VSP data, respectively. In the current application of Q processing technology, it has developed from "near-surface Q compensation, mid-deep Q migration" to "full-layer Q migration". "Near-surface Q compensation, mid-deep Q migration" is to independently model the two Q fields, which are applied to the energy compensation stage and the migration stage respectively; "full-layer Q migration" is to transitionally merge the two Q fields in some way to construct a full-layer Q field with a Q value distribution model covering the entire formation (from shallow to deep), and finally build an overall model and evaluation, and apply it to the migration stage. Compared with the former, the weak reflections in the profile of "full-layer Q migration" are clearer, and the main frequency and bandwidth have been further improved.

[0004] There are two main methods for fusing the near-surface Q field with the mid-deep Q field to obtain the full-layer Q field. One is to fuse the near-surface Q field with the mid-deep Q field along the top surface of the high-speed layer, and the other is to fuse along the zero-difference interface. The method of fusing along the top surface of the high-speed layer does not require picking the fusion surface, but the coupling of the two Q fields at the transition position is poor, and a large smoothing process is required after fusion, which is not conducive to the fine description and evaluation of the overall Q field, resulting in low Q migration accuracy for the full-layer system; the method of fusing along the zero-difference interface has a better coupling effect in the transition zone, but it is necessary to manually pick the zero-difference interface, which is labor-intensive and time-consuming, and the picking accuracy varies due to different manual picking experiences, especially when the grid density of three-dimensional seismic is large, the picking workload increases sharply, resulting in high experimental modeling costs, which is difficult to meet industrial production needs and is not conducive to the comprehensive promotion of full-layer Q field application technology. Summary of the invention

[0005] In order to solve the above-mentioned deficiencies in the prior art, the present invention aims to provide a near-surface and mid-deep Q-field fusion method, system, equipment and medium based on a mask body, and to improve the modeling method and process of the two Q-field fusion, so as to solve the problems of low accuracy and low efficiency in establishing the full-layer Q field in the prior art, and to realize the intelligent and efficient fusion of the near-surface and mid-deep Q fields.

[0006] To achieve the above purpose, the technical method adopted by the present invention is as follows: a near-surface and mid-deep Q field fusion method based on a mask body, comprising the following steps: S1, generate near-surface Q field F1 and mid-deep Q field F2; S2, the near-surface Q field F1 and the mid-deep Q field F2 are weightedly subtracted to form a data volume F3 containing a zero-value interface; S3, remove outliers from the data body F3 to obtain data body F3”; S4, mask the data volume F3" to obtain the near-surface mask volume F 1 / 0 , medium-deep mask body F 0 / 1 ; S5, near-surface Q field F1, mid-depth Q field F2 respectively multiply the mask volume F 1 / 0 、F 0 / 1 Then add them together to complete the fusion of the two Q fields.

[0007] As a limitation, the weighted subtraction formula in S2 is: F3=F1-αF2, where the value range of the weight coefficient α is 0.5 to 1.5.

[0008] As a limitation, the method for removing outliers in S3 is the sliding average method.

[0009] As a further limitation, the moving average method adopts a five-point moving average method, and the calculation steps of the five-point moving average method are: S31, marking the layer data where the zero-value interface in the data volume F3 is located as (x, y, Depth / Time), where Depth is the depth domain data type, and Time is the time domain data type; S32. Calculate the boundary values ​​excluding the two ends of the data. The processing formula is as follows:

[0010] Among them, A (x, y, z n ) represents the zero-value interface sample point value in the data volume F3, A”(x, y, z n ) represents the zero-value interface sample value in the processed data body F3”; the data length of the sample value is k, where 3≤n≤k-2; S33. Calculate the boundary values ​​at both ends of the data. The processing formula is as follows:

[0011] Thus, the unique data volume F3 of the zero-value interface on the vertical single channel is obtained.

[0012] As a limitation, the mask body F in S4 1 / 0 The logic function GE (>=) or LE (<=) is used to perform logical judgment on the data body F3", and the numbers greater than 0 are taken as 1, and the numbers less than 0 are taken as 0; the mask body F 0 / 1 The logical function LOGICALNOT is used to judge it and get the same as F 1 / 0 Complementary mask body F 0 / 1 .

[0013] As another limitation, the mask body F 0 / 1 It uses the logic function GE (>=) or LE (<=) to perform logical judgment on the data body F3", taking the number less than 0 as 1 and the number greater than 0 as 0.

[0014] As a limitation, the near-surface Q field F1 in S1 is obtained by inverting the near-surface velocity field using the cannon first arrival data and fitting it with the dual-well micro-logging data; the mid-deep Q field F2 is obtained by using VSP well information and horizon interpolation.

[0015] The present invention also provides a near-surface and mid-deep Q-field fusion system based on a mask body, comprising: Q field data generation module, used to generate near-surface Q field F1 and mid-deep Q field F2; The data volume F3 acquisition module is used to obtain the zero-value interface including the transition area between the near-surface Q field F1 and the mid-deep Q field F2; The data body F3 processing module is used to remove the abnormal values ​​of the zero-value interface on the vertical single channel, so that the zero-value interface is unique on the vertical single channel, and obtain the data body F3"; Data body F3" mask transformation module, used to obtain mask body F 1 / 0 、F 0 / 1 , mask body F 1 / 0 、F 0 / 1 The data boundary where the 0 and 1 values ​​are located is the ideal fusion surface; The Q field fusion module is used to fuse the two Q fields on the ideal fusion surface, and the point multiplication mask volume F corresponding to the near-surface Q field F1 and the mid-depth Q field F2 is 1 / 0 、F 0 / 1 , we get a data volume F that only contains the corresponding parts and is complementary along the ideal fusion surface 上 and data body F 下 ; The data body F 上 With data body F 下 By adding them together, the fusion of the two Q fields is realized, completing the establishment of the Q field of the entire system.

[0016] The present invention also provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor calls the computer program in the memory to execute any of the above-mentioned mask-based near-surface and mid-deep Q field fusion methods.

[0017] The present invention also provides a computer-readable storage medium, in which a computer program is stored. When the computer program is executed by a processor, it is used to implement any of the above-mentioned mask-based near-surface and mid-deep Q field fusion methods.

[0018] Due to the adoption of the above technical solution, the present invention has the following beneficial effects compared with the prior art: In the fusion of the surface Q field and the mid-deep Q field in the present invention, the concept of zero difference interface fusion is followed, and the two Q fields are weighted subtracted to form a data body F3, in which the zero value interface corresponding to the manually picked best fusion surface is included. At this time, the mask transformation technology in the computer field is used to extend the mask concept to the field of three-dimensional attribute body modeling of seismic data, and the data body to be picked is transformed to form a mask body. The boundary value of the mask body is consistent with the zero difference surface to be picked, so as to ensure that the zero difference interface fusion method has a good coupling effect in the transition zone; when the two Q fields are fused, the two complementary masks are directly combined. The body is fused with the two Q-field data by dot multiplication and then addition. The whole process does not require manual picking of fusion surface data, that is, the traditional "Q field + fusion surface" fusion idea of ​​near-surface and mid-deep Q field fusion is transformed into the "Q field + mask body" idea. The whole fusion process does not involve the "fusion surface", and the "fusion surface" is implied on the data boundary of the mask body, eliminating the work process of manual picking of the fusion surface. While avoiding manual picking errors, it minimizes the cost of manual picking, is beneficial to production and quality control needs, and achieves the expected effect of intelligent and efficient fusion of Q field.

[0019] In summary, the present invention is suitable for the intelligent and efficient integration of near-surface and mid-deep Q fields, meeting the fast-paced and fine exploration needs of three-dimensional production. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] Figure 1 This is a simplified flowchart of a near-surface and mid-deep Q-field fusion method based on a mask in Example 1 of the present invention; Figure 2 is a diagram showing the generation result of the near-surface Q field F1 in step S1 in embodiment 1 of the present invention; Figure 3 This is a diagram showing the generation result of the middle-depth Q field F2 in step S1 in Example 1 of the present invention; Figure 4 is a method flow chart of steps S2 to S5 in Example 1 of the present invention; Figure 5 It is a single-channel display diagram of the near-surface Q field F1 in a certain seismic channel data in Example 1 of the present invention; Figure 6 In Example 1 of the present invention, Figure 5 The corresponding single-channel display of the mid-deep Q field F2 in the seismic data; Figure 7 It is a single channel display diagram of the data body F3 in Example 1 of the present invention and the data body F3" obtained by using the five-point sliding average method; Figure 8 The near-surface Q field mask body F obtained in step S4 of embodiment 1 of the present invention 1 / 0 Schematic diagram of the shape, where: (a) is the mask body F 1 / 0 Single channel schematic diagram, (b) is the mask body F 1 / 0 A schematic diagram of one of the cross sections, (c) is the mask body F 1 / 0 3D schematic diagram; Fig. 9 The middle-depth Q field mask body F obtained in step S4 of embodiment 1 of the present invention 0 / 1 Schematic diagram of the shape, where: (a) is the mask body F 0 / 1 Single channel schematic diagram, (b) is the mask body F 0 / 1 A schematic diagram of one of the cross sections; (c) is the mask body F 0 / 1 3D schematic diagram; Fig.10 It is a three-dimensional schematic diagram of the Q field of the entire layer system after fusion in Example 1 of the present invention; Fig.11 The full-layer Q prestack time migration profile in Example 1 of the present invention, wherein: (a) is a conventional prestack time migration profile (without applying the Q field fusion technology), (b) is a full-layer Q prestack time migration profile obtained by fusing the Q field using a manual picking method, and (c) is a full-layer Q prestack time migration profile obtained by applying the fusion method of the present invention; Fig.12 This is a structural block diagram of a near-surface and mid-deep Q-field fusion system based on a mask body according to Example 2 of the present invention. DETAILED DESCRIPTION

[0021] The preferred embodiments of the present invention are described below in conjunction with the accompanying drawings. It should be understood that the preferred embodiments described herein are only used to illustrate and explain the present invention and do not constitute a limitation of the present invention. Example 1

[0022] A mask-based near-surface and mid-depth Q-field fusion method, such as Figure 1The figure shows a simplified flow chart of the method. Specifically, in this embodiment, the three-dimensional seismic data processing of a certain block in a basin in the central region of China is taken as an example. The area of ​​the work area is about 300km 2 The surface type belongs to the typical loess plateau mountain, and the near-surface height conditions are complex. The specific process of using the fusion method of the present invention is as follows: S1. Generate near-surface Q field F1 and mid-deep Q field F2 based on the seismic data of the work area.

[0023] In this embodiment, the seismic data in the Geo-East seismic data processing system or other existing processing systems are used for data processing. Figure 2 As shown in Figure 1, the obtained dual-well micro-logging data are fitted with the Qv relationship, and the near-surface velocity field is inverted with the cannon first arrival data, and finally the near-surface Q field F1 is formed by fitting; Figure 3 As shown, the mid-deep Q field F2 is obtained by interpolating the VSP well information and the layer information, completing the data preparation of the two Q fields. Of course, the near-surface Q field F1 and the mid-deep Q field F2 in this step can also be generated by other methods in the prior art.

[0024] S2, weighted subtraction of the near-surface Q field F1 and the mid-deep Q field F2 forms a data volume F3, which includes the zero-value interface of the transition region between the near-surface Q field F1 and the mid-deep Q field F2.

[0025] The zero value interface corresponds to the best fusion surface in the existing manual picking method. This interface can be regarded as a layer data (X, Y, Depth / Time). Depth is the depth domain data type, and Time is the time domain data type. You can select Depth or Time data as needed. Most of the layer scattered points on the vertical single channel are unique, and a few points are not unique. The non-unique position is the main reason for the slow speed and poor accuracy of manual picking.

[0026] like Figure 4 As shown, it is a specific method flow chart of steps S2 to S5 of this embodiment. The above-mentioned weighted subtraction formula is: F3=F1-αF2, where the value range of the weight coefficient α is 0.5 to 1.5. When the weight coefficient α is 1, the near-surface Q field F1 and the mid-deep Q field F2 are obtained by direct subtraction. In the subsequent fusion process and Q field migration application, the weight coefficient α can be adjusted through the migration result quality control data, so that the zero-value interface of F3 obtained will move up or down within a certain range in the vertical direction, thereby performing a new round of Q field fusion until the Q field quality control meets the conditions.

[0027] like Figure 5The figure shows a single-channel display of the near-surface Q field F1 in a certain seismic data of this embodiment. The jump points in the figure are abnormal points, which are abnormal values ​​calculated by the first arrival data of the cannon. Figure 6 Shown in Figure 5 The single channel display diagram of the middle-deep Q field F2 in a corresponding seismic data. Taking this single channel data as an example, in this embodiment, α is taken as 1, and the following is obtained by processing F3=F1-F2: Figure 7 A single channel difference of the data volume F3 shown by the black dotted line, Figure 7 In the figure, the intersection of the single-channel data of data body F3 and the horizontal coordinate is the data point of this single-channel data in the zero-value interface. There are black abnormal points on the left side of the intersection (the abnormal points are points below the zero axis). The abnormal points are the places where a small number of points of the above-mentioned layer scatter points on the vertical single channel are not unique, and abnormal value processing is required in the subsequent process.

[0028] S3. Remove outliers from the data volume F3 to make the zero-value interface unique on the vertical single channel, thereby obtaining the data volume F3".

[0029] In this embodiment, the method for removing outliers is a sliding average method, specifically a five-point sliding average method. The calculation steps of the five-point sliding average method are: S31. Mark the layer data where the zero-value interface in the data body F3 is located as (x, y, Depth / Time), where Depth is the depth domain data type and Time is the time domain data type; select either depth domain Depth data or time domain Time data according to the type of seismic data. This embodiment takes Time data as an example.

[0030] S32. Calculate the boundary values ​​excluding the two ends of the data. The processing formula is as follows:

[0031] Among them, A (x, y, z n ) represents the zero-value interface sample point value in the data volume F3, A”(x, y, z n ) represents the zero-value interface sample value in the processed data body F3”; the data length of the sample value is k, where 3≤n≤k-2.

[0032] S33. Calculate the boundary values ​​at both ends of the data. The processing formula is as follows:

[0033] Thus, the unique data volume F3 of the zero-value interface on the vertical single channel is obtained.

[0034] exist Figure 7In the figure, the gray dotted line is the display diagram of the data body F3" obtained by the five-point sliding average method in this embodiment in this single channel. It can be seen that after the outlier removal process, the values ​​on the left side of the intersection of the data body F3 in this single channel are all above the zero axis.

[0035] S4, perform mask transformation on the data body F3" to obtain mask bodies F 1 / 0 、F 0 / 1 , mask body F 1 / 0 、F 0 / 1 The data boundary where the values ​​0 and 1 in are located is the ideal fusion surface.

[0036] Masking technology is a technology that selectively filters or modifies data through binary bit operations, logical operations or specific rules. The core idea is to extract certain parts of the target data through a "template". The resulting mask is the attribute mark of the data.

[0037] In this step, the logical function GE (>=) or LE (<=) (GE represents greater than or equal to, LE represents less than or equal to) is used to perform logical judgment on the data body F3". Figure 4 In the example, the logic function GE is used to set the values ​​greater than 0 to 1 and the values ​​less than 0 to 0, thus obtaining the mask body F. 1 / 0 ; In the mask body F 0 / 1 One method of obtaining the mask body F is to use the logical function LOGICALNOT (logical negation) to obtain the mask body F. 1 / 0 Make a judgment and get the same as F 1 / 0 Complementary mask body F 0 / 1 Another mask body F 0 / 1 The method of obtaining is: use the logic function GE (>=) or LE (<=) to directly perform logical judgment on the data body F3", and take the number less than 0 as 1 and the number greater than 0 as 0.

[0038] like Figure 8 As shown in the figure, the near-surface Q field mask body F is obtained after mask transformation. 1 / 0 Schematic diagrams of various forms, Figure 8 (a) is the mask body F 1 / 0 A single-lane schematic diagram, where the value is 1 on the left side of the junction and 0 on the right side of the junction; Figure 8 (b) is the mask body F 1 / 0 A schematic diagram of one of the cross sections. Figure 8 (c) is the mask body F 1 / 0 3D schematic diagram; formed by superposition of single channel data Figure 8 (b) Figure 8 (c) Masked body morphology.

[0039] like Fig. 9As shown, after mask transformation, the mid-depth Q field mask body F is obtained. 0 / 1 Schematic diagrams of various forms, Fig. 9 (a) is the mask body F 0 / 1 Single-lane schematic diagram, where the value is 0 on the left side of the junction and 1 on the right side of the junction; Fig. 9 (b) is the mask body F 0 / 1 A schematic diagram of one of the cross sections. Fig. 9 (c) is the mask body F 0 / 1 3D schematic diagram; formed by superposition of single channel data Fig. 9 (b) Fig. 9 (c) The mask shape of the mid-depth Q field mask F 0 / 1 and the near-surface Q-field mask F 1 / 0 To complement each other.

[0040] S5, near-surface Q field F1, mid-depth Q field F2 respectively multiply the mask volume F 1 / 0 、F 0 / 1 Then add them together to complete the fusion of the two Q fields.

[0041] Specifically, the calculation steps of this embodiment are as follows: Figure 4 As shown in the figure, the near-surface Q field F1 is multiplied by the mask volume F 1 / 0 , get the data body F 上 , multiply the middle and deep Q field F2 by the mask volume F 0 / 1 , and obtain the data body F 下 ; F 上 and F 下 Contains only the corresponding parts and complements along the ideal fusion surface, and then the data volume F 上 With data body F 下 The above calculation process adopts the calculation method of "multiplying and adding", and finally obtains the following Fig.10 The three-dimensional schematic diagram of the fused full-layer Q field is shown in the figure, which realizes the fusion of the two Q fields and completes the establishment of the full-layer Q field.

[0042] In the processing of three-dimensional seismic data of a certain block in a certain basin of the present embodiment, the method of directly using high-speed top surface fusion in the Q field fusion stage is adopted, resulting in poor coupling and low precision in the transition zone, while the traditional zero difference surface fusion manual picking method is time-consuming and there are manual picking errors. The single round fusion takes an average of 30 hours, and multiple rounds of fusion are often required in actual operation. The Q field fusion process lasts for half a month, which cannot meet the needs of industrial production. Using the fusion method of the present invention, the average single round takes only 20 minutes, which not only greatly shortens the processing time, but also takes into account the accuracy of zero difference surface fusion, while also avoiding the error of layer picking that varies from person to person, which is beneficial to production and quality control requirements, and achieves the expected effect of intelligent and efficient fusion of Q field.

[0043] In the subsequent migration application of the full-layer Q field after fusion of this embodiment, the quality of the Q time migration profile can be compared with the results of conventional time migration (without applying the Q field fusion technology), so as to synchronize it as a quality control method for the Q field fusion effect, thereby quality controlling the fusion Q field effect from multiple dimensions.

[0044] In the quality control process, the weight coefficient α can be adjusted by offsetting the result quality control data, so that the zero-value interface of F3 can move up or down within a certain range in the vertical direction, thereby performing a new round of Q-field fusion until the Q-field quality control meets the conditions. This method can be used multiple times in this iterative process and has high operating efficiency. Fig.11 The full-layer Q prestack time migration profile in Example 1 of the present invention, wherein: (a) is a conventional prestack time migration profile (Q field fusion technology is not applied), and the migration profile effect at the arrow is poor; (b) is a full-layer Q prestack time migration profile obtained by fusing the Q field by the manual picking method, and the migration profile effect at the arrow is good; (c) is a full-layer Q prestack time migration profile obtained by applying the fusion method of the present invention, and the migration profile effect at the arrow is good. It can be seen from this figure that the quality of the full-layer Q field time migration profile fused by the method of the present invention can reach the quality of the migration profile using the manual picking method, but the method greatly shortens the processing time, greatly improves the implementation effect of the full-layer Q field application technology, and meets the needs of industrial production. Example 2

[0045] A near-surface and mid-depth Q-field fusion system based on a mask body, such as Fig.12 As shown, including: Q field data generation module, used to generate near-surface Q field F1 and mid-deep Q field F2; The data volume F3 acquisition module is used to obtain the zero-value interface including the transition area between the near-surface Q field F1 and the mid-deep Q field F2; The data body F3 processing module is used to remove the abnormal values ​​of the zero-value interface on the vertical single channel, so that the zero-value interface is unique on the vertical single channel, and obtain the data body F3"; Data body F3" mask transformation module, used to obtain mask body F 1 / 0 、F 0 / 1 , mask body F 1 / 0 、F 0 / 1 The data boundary where the 0 and 1 values ​​are located is the ideal fusion surface; The Q field fusion module is used to fuse the two Q fields on the ideal fusion surface, and the point multiplication mask volume F corresponding to the near-surface Q field F1 and the mid-depth Q field F2 is 1 / 0 、F 0 / 1 , we get a data volume F that only contains the corresponding parts and is complementary along the ideal fusion surface 上 and data body F下 ; The data body F 上 With data body F 下 By adding them together, the fusion of the two Q fields is realized, completing the establishment of the Q field of the entire system. Example 3

[0046] A computer device includes a memory and a processor, wherein the memory stores a computer program, and the processor calls the computer program in the memory to execute the computer program to implement the method in embodiment 1. Example 4

[0047] This embodiment provides a computer-readable storage medium, in which a computer program is stored. When the computer program is executed by a processor, it is used to implement the method in Embodiment 1.

[0048] Among them, the computer-readable storage medium can be a computer storage medium or a communication medium. The communication medium includes any medium that facilitates the transmission of a computer program from one place to another. The computer storage medium can be any available medium that can be accessed by a general-purpose or special-purpose computer. For example, a computer-readable storage medium is coupled to a processor so that the processor can read information from the computer-readable storage medium and write information to the computer-readable storage medium. Of course, the computer-readable storage medium can also be a component of the processor. The processor and the computer-readable storage medium can be located in an application-specific integrated circuit (ASIC). In addition, the ASIC can be located in a user device. Of course, the processor and the computer-readable storage medium can also exist in a communication device as discrete components. Specifically, the computer-readable storage medium can be implemented by any type of volatile or non-volatile storage device or a combination thereof, such as a static random access memory (SRAM), an electrically erasable programmable read-only memory (EEPROM), an erasable programmable read-only memory (EPROM), a programmable read-only memory (PROM), a read-only memory (ROM), a magnetic memory, a flash memory, a magnetic disk or an optical disk, etc. The storage medium can be any available medium that can be accessed by a general-purpose or special-purpose computer.

[0049] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present disclosure, rather than to limit them. Although the present disclosure 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 by 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 disclosure.

Claims

1. A near-surface and mid-deep Q-field fusion method based on a mask volume, characterized in that: The method comprises the following steps: S1, generate near-surface Q field F1 and mid-deep Q field F2; S2, the near-surface Q field F1 and the mid-deep Q field F2 are weightedly subtracted to form a data volume F3 containing a zero-value interface; S3, remove outliers from the data body F3 to obtain data body F3”; S4, mask the data volume F3" to obtain the near-surface mask volume F 1 / 0 , medium-deep mask body F 0 / 1 ; S5, near-surface Q field F1, mid-depth Q field F2 respectively multiply the mask volume F 1 / 0 、F 0 / 1 Then add them together to complete the fusion of the two Q fields.

2. The method for fusion of near-surface and mid-deep Q fields based on a mask volume according to claim 1, characterized in that: The weighted subtraction formula in S2 is: F3=F1-αF2, where the value range of the weight coefficient α is 0.5 to 1.

5.

3. The method for fusion of near-surface and mid-deep Q fields based on a mask volume according to claim 1, characterized in that: The method for removing outliers in S3 is the sliding average method.

4. The method for fusion of near-surface and mid-deep Q fields based on a mask volume according to claim 3, characterized in that: The sliding average method adopts a five-point sliding average method, and the calculation steps of the five-point sliding average method are: S31, marking the layer data where the zero-value interface in the data volume F3 is located as (x, y, Depth / Time), where Depth is the depth domain data type, and Time is the time domain data type; S32. Calculate the boundary values ​​excluding the two ends of the data. The processing formula is as follows: Among them, A (x, y, z n ) represents the zero-value interface sample point value in the data volume F3, A”(x, y, z n ) represents the zero-value interface sample value in the processed data body F3”; the data length of the sample value is k, where 3≤n≤k-2; S33. Calculate the boundary values ​​at both ends of the data. The processing formula is as follows: ， ， ， ， Thus, the unique data volume F3 of the zero-value interface on the vertical single channel is obtained.

5. A mask-based near-surface and mid-deep Q field fusion method according to any one of claims 1 to 4, characterized in that: The mask body F in S4 1 / 0 The logic function GE (>=) or LE (<=) is used to perform logical judgment on the data body F3", and the numbers greater than 0 are taken as 1, and the numbers less than 0 are taken as 0; the mask body F 0 / 1 The logical function LOGICALNOT is used to judge it and get the same as F 1 / 0 Complementary mask body F 0 / 1 .

6. The method for fusion of near-surface and mid-deep Q fields based on a mask volume according to claim 5, characterized in that: The mask body F 0 / 1 It uses the logic function GE (>=) or LE (<=) to perform logical judgment on the data body F3", taking the number less than 0 as 1 and the number greater than 0 as 0.

7. The mask-based near-surface and mid-deep Q field fusion method according to claim 6, characterized in that: The near-surface Q field F1 in S1 is obtained by inverting the near-surface velocity field using the first arrival data of the cannon and fitting it with the micro-logging data of the two wells; the mid-deep Q field F2 is obtained by using the VSP well information and the interpolation of the horizon.

8. A near-surface and mid-deep Q-field fusion system based on a mask, characterized in that: include: Q field data generation module, used to generate near-surface Q field F1 and mid-deep Q field F2; The data volume F3 acquisition module is used to obtain the zero-value interface including the transition area between the near-surface Q field F1 and the mid-deep Q field F2; The data body F3 processing module is used to remove the abnormal values ​​of the zero-value interface on the vertical single channel, so that the zero-value interface is unique on the vertical single channel, and obtain the data body F3"; The data volume F3" mask transformation module is used to obtain the near-surface mask volume F 1 / 0 , medium-deep mask body F 0 / 1 , F 1 / 0 、F 0 / 1 The data boundary where the 0 and 1 values ​​are located is the ideal fusion surface; The Q field fusion module is used to fuse the two Q fields on the ideal fusion surface, multiplying the points corresponding to the near-surface Q field F1 and the mid-depth Q field F2 by the near-surface mask volume F 1 / 0 , medium-deep mask body F 0 / 1 , we get a data volume F that only contains the corresponding parts and is complementary along the ideal fusion surface 上 and data body F 下 ; The data body F 上 With data body F 下 By adding them together, the fusion of the two Q fields is realized, completing the establishment of the Q field of the entire system.

9. A computer device, characterized in that: The invention comprises a memory and a processor, wherein the memory stores a computer program, and the processor calls the computer program in the memory to execute a mask-based near-surface and mid-deep Q field fusion method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, which, when executed by a processor, is used to implement a mask-based near-surface and mid-deep Q field fusion method as described in any one of claims 1 to 7.

Images