Seismic plane attribute value domain adjusting method and device based on changed distribution curve

Abstract

The invention discloses a seismic plane attribute value domain adjustment method and device based on change of a distribution curve, and relates to the field of seismic exploration, and the key points of the technical scheme are as follows: through the analysis of a seismic plane attribute value distribution curve, marking a target information value domain, determining a seismic plane attribute value domain CAD and an adjustment value CRS which need to be highlighted, and determining an adjustment direction; a distribution matching technology is utilized, a seismic attribute distribution curve is carried out to carry out telescopic iteration, a distribution curve meeting requirements is determined, seismic plane attribute value domain distribution adjustment is carried out, and a result can effectively highlight target information concerned by exploration personnel. Meanwhile, the method can effectively support seismic plane attribute connecting piece splicing.

Description

Technical Field

[0001] The present invention relates to the field of seismic exploration technology, and more particularly to a method for adjusting a seismic plane attribute value range based on changing a distribution curve. Background Art

[0002] Before drilling operations, multiple 3D seismic layer data in the target area need to be spliced ​​into a map for attribute analysis, reservoir inversion, fracture prediction, and comprehensive evaluation, thereby improving the success rate of drilling. Due to the differences in the time and parameters of 3D seismic data acquisition and processing, the data boundaries between different work areas are obvious and the closure difference of the layers is large, which greatly increases the difficulty of cross-work area contiguous interpretation and mapping, is not conducive to conducting structural analysis and research from a regional contiguous perspective, and makes it difficult to understand the regional reservoir formation laws and determine key targets.

[0003] In this regard, Chinese patent CN112269208A proposes a method for mapping seismic attributes across three-dimensional seismic work areas, which includes the following steps: Step 1: Observe the plane characteristics of the selected attribute layers in combination with the geological task, and remove the boundary outliers caused by the boundary effect of the seismic data; Step 2: Normalize the value range of the attribute layers of each three-dimensional seismic work area to make its value range uniform; Step 3: Observe the concentrated areas of the outlier values ​​of each attribute layer, and calculate the layer value range again to make the outlier value range roughly the same; Step 4: Splice the adjusted layers together.

[0004] When adjusting the plane attribute value range, the above patent first unifies the value range of each 3D seismic work area, and then makes the abnormal value range of each attribute layer abnormal value concentration area roughly the same. It can be seen that the adjustment of the value range in this patent is an overall adjustment of the plane attribute value range, which is difficult to highlight a local plane attribute that the exploration personnel are more concerned about according to the research needs.

[0005] To this end, the present invention hopes to provide a seismic plane attribute value range adjustment technology based on changing the distribution curve, so as to achieve local adjustment of the seismic plane attribute value range and highlight the local plane attributes required by the exploration personnel. Summary of the invention

[0006] In view of the shortcomings of the prior art, the present invention aims to provide a method for adjusting the value range of seismic plane attributes based on changing the distribution curve, which analyzes the distribution curve of seismic plane attribute values, marks the value range of target information, and determines the seismic plane attribute value range C that needs to be highlighted. AD and adjustment value C RSAnd clarify the adjustment direction, use the distribution matching technology, carry out the expansion and contraction iteration of the seismic attribute distribution curve, determine the distribution curve that meets the requirements, and carry out the distribution adjustment of the seismic plane attribute value range. The results can effectively highlight the target information that the exploration personnel are concerned about. At the same time, the present invention can also effectively support the splicing of seismic plane attributes.

[0007] To achieve the above object, the present invention provides the following technical solution: a method for adjusting the range of seismic plane attributes based on changing the distribution curve, which comprises the following steps:

[0008] S1. Statistic the seismic plane attribute value H of the target layer, obtain the value range distribution curve, and determine the minimum value C of the value range distribution curve min and the maximum value C max ;

[0009] S2. Given the earthquake plane attribute value range C that needs to be highlighted AD and adjustment value C RS , C AD In C min to C max within the scope of

[0010] S3. Obtain the adjusted seismic plane attribute value H through distribution matching technology new The value range distribution of

[0011] The formula for the distribution matching technique is as follows:

[0012] H=C AD When H new =C RS ;

[0013] H>C AD hour,

[0014] H<C AD hour,

[0015] The present invention is further configured as follows: C min <C RS <C max .

[0016] The present invention is further configured to include: S4: obtaining H and H new The value range distribution curve of is given with a standard value A%; calculate H and H new The difference between the value range distribution curves is N, and whether N is less than A%, if N is less than A%, then output H new Otherwise, H new As a new H, repeat S1-S4 until N is less than A%.

[0017] The present invention is further configured as follows: N is calculated by: new Collect Q points on the distribution curve of the range of the value of and calculate H and H according to the coordinates of these points new The Euclidean distance and H or H new The root mean square of H and H new Subtract the Euclidean distance of , and subtract the difference from the calculated root mean square to calculate N.

[0018] The present invention is further configured as follows: in S4, the horizontal coordinates of the Q points collected on the H value range distribution curve are respectively equal to H new The horizontal coordinates of the Q points collected on the range distribution curve.

[0019] The present invention is further configured such that: a root mean square of H is used in S4.

[0020] The present invention is further configured as follows: 2≤Q≤256 in S4.

[0021] The present invention is further configured as follows: A% in S4 is equal to 3%.

[0022] The present invention is further configured as follows: H and H new The minimum accuracy is 0.001.

[0023] The present invention is further configured as follows: the value range distribution curve is obtained based on a histogram of seismic plane attributes.

[0024] The present invention also provides a device for adjusting the value range of seismic plane attributes based on changing the distribution curve, which includes a storage medium on which a computer program is stored; and also includes a processor, which is used to implement the above-mentioned method for adjusting the value range of seismic plane attributes based on changing the distribution curve when executing the computer program.

[0025] In summary, compared with the prior art, the present invention has the following beneficial effects: the present invention analyzes the distribution curve of seismic plane attribute values, marks the target information value range, and determines the seismic plane attribute value range C that needs to be highlighted. AD and adjustment value C RS And clarify the adjustment direction, use the distribution matching technology, carry out the expansion and contraction iteration of the seismic attribute distribution curve, determine the distribution curve that meets the requirements, and carry out the distribution adjustment of the seismic plane attribute value range. The results can effectively highlight the target information that the exploration personnel are concerned about. At the same time, the present invention can also effectively support the splicing of seismic plane attributes. BRIEF DESCRIPTION OF THE DRAWINGS

[0026] Figure 1 It is a schematic diagram of the process of an embodiment;

[0027] Figure 2is the histogram of the earthquake plane attribute H in S1;

[0028] Figure 3 Histogram of seismic plane attributes adjusted for traditional methods;

[0029] Figure 4 The seismic plane attribute H after adjustment for S3 new Histogram;

[0030] Figure 5 It is the plane distribution diagram of the value range of the earthquake plane attribute H;

[0031] Figure 6 The plane distribution diagram of the seismic plane attribute value range after adjustment for the traditional method;

[0032] Figure 7 It is the plane distribution diagram of the earthquake plane attribute value range adjusted according to S3;

[0033] Figure 8 This is a layer mosaic map after adjustment using the embodiment. DETAILED DESCRIPTION

[0034] The technical solution of the present invention will be clearly described below in conjunction with the accompanying drawings. Obviously, the described embodiments are not all embodiments of the present invention, and all other embodiments obtained by ordinary technicians in the field without making creative work are within the scope of protection of the invention.

[0035] It should be noted that the directions or positional relationships indicated by terms such as “center”, “up”, “down”, “horizontal”, “left”, “right”, “front”, “back”, “lateral” and “longitudinal” are based on the directions or positional relationships shown in the accompanying drawings and are only for the convenience of describing the present invention and simplifying the description. They do not indicate or imply that the device or element referred to must have a specific direction, be constructed and operated in a specific direction, and therefore should not be understood as a limitation on the present invention.

[0036] Example

[0037] like Figure 1 FIG. 1 is a flow chart of a method for adjusting the value range of seismic plane attributes based on changing the distribution curve in a preferred embodiment of the present invention, which includes the following steps:

[0038] S1. Prepare the horizons for which the seismic plane attribute value range needs to be adjusted, perform statistics on a seismic plane attribute value H of the target horizon, obtain the value range distribution curve, and determine the minimum value C of the value range distribution curve. min and the maximum value C max The horizontal axis of the range distribution curve represents the range of the seismic plane attribute, and the vertical axis represents the frequency of data within the range. min and Cmax Both are the horizontal axes of the value range distribution curve.

[0039] S2. Given the earthquake plane attribute value range C that needs to be highlighted AD and adjustment value C RS , C AD In C min to C max within the range.

[0040] S3. Calculate the adjusted seismic plane attribute value H through distribution matching technology new , and plot the earthquake plane attribute value H new The value range distribution curve of .

[0041] The formula for the distribution matching technique is as follows:

[0042] H=C AD When H new =C RS .

[0043] H>C AD hour,

[0044] H<C AD hour,

[0045] This embodiment marks the target information value range based on the seismic plane attribute value distribution curve, and determines the seismic plane attribute value range C that needs to be highlighted. AD and adjustment value C RS And clarify the adjustment direction, use distribution matching technology to carry out expansion and contraction iteration of the seismic attribute distribution curve, carry out distribution adjustment of the seismic plane attribute value range, and realize local adjustment of the seismic plane attribute value range. The results can effectively highlight the target information that the exploration personnel are concerned about.

[0046] Specifically, in this embodiment, C min <C RS <C max .

[0047] This embodiment further includes S4: obtaining H and H new The value range distribution curve of is given with a standard value A%; calculate H and H new The difference between the value range distribution curves is N, and whether N is less than A%, if N is less than A%, then output H new Otherwise, H new As the new H, S1-S4 are repeated until N is less than A%, so as to avoid excessive distortion of the adjusted seismic plane attribute value range distribution. In this embodiment, A% in S4 is equal to 3%.

[0048] Specifically, the calculation method of N is:new Collect Q points on the distribution curve of the range of the value of and calculate H and H according to the coordinates of these points new The Euclidean distance and H or H new The root mean square of H and H new The Euclidean distance of H and H is subtracted from the obtained difference and the obtained root mean square to obtain N. The Euclidean distance and root mean square can be used to calculate the new The similarity between the value range distribution curves is quantitatively evaluated without relying on domain expertise, the calculation is efficient and the results are intuitive.

[0049] Specifically, in S4, the horizontal coordinates of the Q points collected on the H value range distribution curve are equal to H new The horizontal coordinates of the Q points collected on the range distribution curve ensure that H and H new The point data collected on the value range distribution curve correspond to each other, avoiding errors in the calculated N due to misalignment of the collected point data.

[0050] Specifically, S4 adopts the root mean square as the root mean square of H, emphasizing the importance of the value range distribution curve of H, so that the subsequent value range distribution adjustment uses the value range distribution curve of H as a reference standard.

[0051] Specifically, in S4, 2≤Q≤256. In modern computer architecture, 8-bit bytes are usually used as basic storage units. When Q is not greater than the maximum value of eight-bit binary, the coordinate data of Q points can be directly mapped to the binary representation in the computer. This can simplify the data representation and calculation process, without the need for additional data conversion or processing, and improve the calculation efficiency. In this embodiment, Q is equal to 256 to provide a more comprehensive and accurate data sample to better reflect H and H new The overall characteristics and trends of the value range distribution curve can reduce the error introduced by N due to insufficient number of Q.

[0052] Specifically, in this embodiment, H and H new The minimum accuracy is 0.001.

[0053] Specifically, the range distribution curve is obtained based on the histogram of seismic plane attributes. Figure 2 That is the histogram of the earthquake plane attribute H in S1. Figure 4 As shown in FIG. 1 , the earthquake plane attribute H finally outputted in this embodiment is new Histogram, with Figure 3 Compared with the histogram adjusted by the traditional method, it can be seen that the above range distribution matching technology can effectively highlight the target information that the exploration personnel are concerned about. Figure 5-7 As shown in Figure 2, this effect is more intuitive in the range plane distribution diagram.

[0054] This embodiment also provides a device for adjusting the value range of seismic plane attributes based on changing the distribution curve, which includes a storage medium on which a computer program is stored; and also includes a processor, which is used to implement the above-mentioned method for adjusting the value range of seismic plane attributes based on changing the distribution curve when executing the computer program.

[0055] In summary, this embodiment analyzes the distribution curve of seismic plane attribute values, marks the target information value range, and determines the seismic plane attribute value range C that needs to be highlighted. AD and adjustment value C RS And clarify the adjustment direction, use distribution matching technology to carry out expansion and contraction iteration of seismic attribute distribution curve, determine the distribution curve that meets the requirements, and carry out seismic plane attribute value range distribution adjustment. The results can effectively highlight the target information that exploration personnel are concerned about. Figure 8 As shown, the present invention can also effectively support the splicing of seismic plane attributes.

[0056] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. A method for adjusting the range of seismic plane attributes based on changing the distribution curve, characterized in that: The steps include: S1. Statistic the seismic plane attribute value H of the target layer, obtain the value range distribution curve, and determine the minimum value C of the value range distribution curve min and the maximum value C max ; S2. Given the earthquake plane attribute value range C that needs to be highlighted AD and adjustment value C RS ; S3. Obtain the adjusted seismic plane attribute value H through distribution matching technology new The value range distribution of The formula for the distribution matching technique is as follows: H=C AD When H new =C RS ; H>C AD hour, H <C AD hour, 2. The method for adjusting the range of seismic plane attributes based on changing the distribution curve according to claim 1, characterized in that: C min <C RS <C max 。 3. The method for adjusting the range of seismic plane attributes based on changing the distribution curve according to claim 2, characterized in that: Also includes S4: Get H and H new The value range distribution curve of is given with a standard value A%; calculate H and H new The difference between the value range distribution curves is N, and whether N is less than A%, if N is less than A%, then output H new Otherwise, H new As a new H, repeat S1-S4 until N is less than A%.

4. The method for adjusting the range of seismic plane attributes based on changing the distribution curve according to claim 3, characterized in that: The calculation method of N is: new Collect Q points on the distribution curve of the range of the value of and calculate H and H according to the coordinates of these points new The Euclidean distance and H or H new The root mean square of H and H new Subtract the Euclidean distance of , and subtract the difference from the calculated root mean square to calculate N.

5. The method for adjusting the range of seismic plane attributes based on changing the distribution curve according to claim 4 is characterized in that: In S4, the horizontal coordinates of the Q points collected on the H value range distribution curve are equal to H new The horizontal coordinates of the Q points collected on the range distribution curve.

6. The method for adjusting the range of seismic plane attributes based on changing the distribution curve according to claim 4, characterized in that: In S4, the root mean square is H.

7. The method for adjusting the range of seismic plane attributes based on changing the distribution curve according to claim 4, characterized in that: 2≤Q≤256 in S4.

8. The method for adjusting the range of seismic plane attributes based on changing the distribution curve according to any one of claims 3 to 7, characterized in that: In S4, A% is equal to 3%.

9. The method for adjusting the range of seismic plane attributes based on changing the distribution curve according to claim 1, characterized in that: H and H new The minimum accuracy is 0.

001.

10. The method for adjusting the range of seismic plane attributes based on changing the distribution curve according to claim 1, characterized in that: The range distribution curve is obtained based on the histogram of seismic plane attributes.

11. A device for adjusting the range of earthquake plane attributes based on changing the distribution curve, characterized in that: It includes a storage medium on which a computer program is stored; and it also includes a processor, which is used to implement the method for adjusting the value range of seismic plane attributes based on changing the distribution curve as described in any one of claims 1-10 when executing the computer program.

Images