Stratigraphic comparison method and system for heterogeneous multilayer sandstone oil field and computer equipment

Abstract

The invention provides a stratigraphic contrast method and system for a heterogeneous multilayer sandstone oil field and computer equipment, and belongs to the field of oil-gas field development. The stratigraphic contrast method for the heterogeneous multi-layer sandstone oil field comprises the following steps: dividing different levels of standard layer boundary lines for each development well in a target area; determining the sedimentary environment of a layer section between two adjacent standard layer boundary lines of each development well according to preset different sedimentary environment stratum development characteristic evaluation standards; dividing a boundary line of a continuous deposition unit of a layer section according to a deposition environment of the layer section between two adjacent standard layer boundary lines of each development well to obtain a preliminary stratum boundary line division result; and establishing a three-dimensional space model in the stratum section based on the preliminary stratum boundary division result, and correcting the stratum boundary division result through the three-dimensional space model in the stratum section to obtain an optimal stratum boundary division result.

Description

Technical Field

[0001] The present invention relates to the field of oil and gas field development, and particularly to a stratigraphic correlation method for heterogeneous multi-layer sandstone oilfields, a stratigraphic correlation system for heterogeneous multi-layer sandstone oilfields, and a computer device. Background Art

[0002] Stratigraphic correlation is the basis of fine reservoir description. Regardless of the type of reservoir, the structural pattern of the strata is the basic content for understanding the reservoir and the key to describing the distribution of the reservoir. With the continuous deepening of oilfield development and the continuous increase in the number of development wells, it is difficult to rely solely on manual stratigraphic correlation to meet the needs of efficient oilfield development. Especially for heterogeneous multi-layer sandstone oilfields, there are nearly a hundred vertical layers, and it is difficult to accurately identify and correlate the entire formation system.

[0003] Traditional stratigraphic correlation methods are based on high-resolution sequence stratigraphy. By using the characteristics of well logging curves that can reflect the preservation degree and accumulation rate of sediments, such as the facies sequence, stratigraphic interface, and facies differentiation, the symmetry degree of the stratigraphic cycle, the pattern of cycle thickening or thinning, and the amplitude and direction of facies dislocation at the position crossing the genetic stratigraphic boundary are identified, and then the stratigraphic boundary is divided. However, traditional stratigraphic correlation methods are relatively simple and mainly deal with homogeneous thick sections of single wells. For heterogeneous multi-well sandstone oilfields, the accuracy of comparison and identification is insufficient. Summary of the Invention

[0004] To solve the above technical defects, the present invention provides a stratigraphic correlation method, system, and computer device for heterogeneous multi-layer sandstone oilfields. The stratigraphic correlation method for heterogeneous multi-layer sandstone oilfields divides the development wells into multiple segments, then determines the sedimentary environment of the segments, and uses different methods to divide the sedimentary units within the segments according to the sedimentary environment of the segments to obtain a preliminary stratigraphic division result, and then corrects the preliminary division result to obtain an optimal stratigraphic boundary division result. This method first divides the development wells into multiple segments, and then divides the sedimentary units within the segments; each segment is treated as a unit, and dividing a segment improves the accuracy of stratigraphic division.

[0005] The first aspect of the present invention provides a stratigraphic correlation method for heterogeneous multi-layer sandstone oilfields, including:

[0006] Dividing standard layer boundaries of different levels for each development well in the target area;

[0007] Determining the sedimentary environment of the segments between two adjacent standard layer boundaries of each development well according to the preset evaluation criteria for the development characteristics of strata in different sedimentary environments;

[0008] Dividing the boundaries of continuous sedimentary units within the segment according to the sedimentary environment of the segment between the two adjacent standard layer boundaries to obtain a preliminary stratigraphic boundary division result;

[0009] Establish a three-dimensional spatial model within the layer segment based on the preliminary formation boundary division result, and correct the formation boundary division result through the three-dimensional spatial model within the layer segment to obtain the optimal formation boundary division result.

[0010] In the embodiment of the present invention, the division of different levels of standard layer boundaries for each development well in the target area includes:

[0011] Construct multiple different levels of standard layer systems according to the core data of the target area;

[0012] Construct a multi-level standard layer division model using an image segmentation method based on a convolutional neural network;

[0013] Extract the standard layer characteristic logging curves of each level of standard layer system, and use the standard layer logging curves of each level of standard layer system as the training set to input into the multi-level standard layer division model for training;

[0014] Obtain the original logging curves of each development well and input them into the trained multi-level standard layer division model to obtain the different levels of standard layer boundaries of each development well.

[0015] In the embodiment of the present invention, the method further includes: correcting the different levels of standard layer boundaries of each development well.

[0016] In the embodiment of the present invention, the determination of the sedimentary environment of the layer segment between two adjacent standard layer boundaries of each development well according to the preset evaluation criteria for the formation development characteristics of different sedimentary environments includes:

[0017] Select the evaluation standard parameters for the formation development characteristics for characterizing the sedimentary environment, and construct the evaluation criteria for the formation development characteristics of different sedimentary environments according to the evaluation standard parameters for the formation development characteristics;

[0018] Obtain the logging interpretation data within each layer segment of each well in the target area;

[0019] Calculate the parameter values of the formation development characteristic evaluation parameters according to the logging interpretation data within each layer segment;

[0020] Determine the sedimentary environment of each layer segment from the evaluation criteria for the formation development characteristics of different sedimentary environments based on the calculated parameter values of the formation development characteristic evaluation parameters.

[0021] In the embodiment of the present invention, the formation development characteristic evaluation parameters include: average effective thickness, average effective permeability, sand-to-ground ratio, drilling encounter rate, layering coefficient;

[0022] Among them, the calculation formula for the average effective thickness is:

[0023] The calculation formula for the average effective permeability is as follows:

[0024] The calculation formula for the sand-to-ground ratio is as follows:

[0025] The calculation formula for the drilling encounter rate is as follows:

[0026] The calculation formula for the layering coefficient is as follows:

[0027] Wherein, i is the i-th interval of the development well, j is the j-th sedimentary unit within the i-th interval, n is the number of sedimentary units, m is the number of wells encountering intervals in the target area, h ei is the average effective thickness of the interval, K ei is the average effective permeability of the interval, D i is the drilling encounter rate, Q i is the layering coefficient, h ej is the effective thickness of the j-th sedimentary unit within the interval, h j is the formation thickness of the j-th sedimentary unit within the interval, h syj is the sandstone thickness of the j-th sedimentary unit within the interval, k j is the effective permeability of the j-th sedimentary unit within the interval, N i is the number of wells encountering the interval, T is the total number of wells in the block, E i is the number of effective sandstone layers within the interval.

[0028] In the embodiment of the present invention, the method for dividing the boundaries of continuous sedimentary units within the interval according to the sedimentary environment between the boundaries of two adjacent standard layers to obtain a preliminary formation boundary division result includes:

[0029] Select multiple types of well logging curves for formation division and correlation features, and synthesize the multiple types of well logging curves for formation division and correlation features into a fusion feature curve;

[0030] According to the sedimentary environment of each interval, apply the fusion feature curve to divide the boundaries of continuous sedimentary units within the current interval to obtain a preliminary formation boundary division result.

[0031] In the embodiment of the present invention, the method for dividing multiple sedimentary units within the current interval according to the sedimentary environment of each interval by applying the fusion feature curve to obtain a preliminary formation division result includes:

[0032] If the sedimentary environment of the layer section is the frontal facies, then according to the fusion characteristic curve, the wavelet transform method and the local INPEFA logging cycle analysis method are used to divide the boundaries of the continuous sedimentary units within the current layer section to obtain a preliminary formation division result;

[0033] If the sedimentary environment of the layer section is the plain facies, then according to the fusion characteristic curve, the sand body superposition pattern constraint method and the continuous formation correlation division method of dynamic programming are used to divide the boundaries of the continuous sedimentary units within the current layer section to obtain a preliminary formation division result.

[0034] In the embodiment of the present invention, establishing a three-dimensional space model within the layer section based on the preliminary formation boundary division result, and correcting the formation boundary division result through the three-dimensional space model within the layer section to obtain an optimal formation boundary division result, including:

[0035] According to the preliminary formation boundary division result, establish a three-dimensional space model within the layer section for characterizing the continuous formation distribution within the layer section;

[0036] Correct the local cross-stratum abnormal points on the three-dimensional space model within the layer section to obtain a preliminarily corrected formation division result;

[0037] Establish an isopach map of the plane formation thickness of each sedimentary unit within the layer section, and obtain and mark the local abnormal points of the preliminarily corrected formation division result in the isopach map of the plane formation thickness within the layer section;

[0038] According to the sedimentary environment of the sedimentary unit, based on the three-dimensional space model within the layer section and the isopach map of the plane thickness of each sedimentary unit within the layer section, correct the local abnormal points to obtain an optimal formation boundary division result.

[0039] The second aspect of the present invention provides a formation correlation system for a heterogeneous multi-layer sandstone oilfield, including:

[0040] A layer section division unit for dividing the standard layer boundaries of different levels for each development well within the target area;

[0041] A layer section sedimentary environment confirmation unit for determining the sedimentary environment of the layer section between the adjacent two standard layer boundaries of each development well according to the preset evaluation criteria for the formation development characteristics of different sedimentary environments;

[0042] A formation boundary division unit for dividing the boundaries of the continuous sedimentary units within the layer section according to the sedimentary environment of the layer section between the adjacent two standard layer boundaries of each development well to obtain a preliminary formation boundary division result;

[0043] A correction unit, configured to establish a three-dimensional spatial model within a layer segment based on the preliminary formation boundary division result, and correct the formation boundary division result through the three-dimensional spatial model within the layer segment to obtain an optimal formation boundary division result.

[0044] The third aspect of the present invention provides a computer device, including:

[0045] A memory;

[0046] A processor; and

[0047] A computer program;

[0048] Wherein, the computer program is stored in the memory and is configured to be executed by the processor to implement a formation correlation method for a heterogeneous multi-layer sandstone oilfield as described above.

[0049] The formation correlation method for a heterogeneous multi-layer sandstone oilfield proposed by the present invention divides development wells into multiple layer segments, then determines the sedimentary environment of the layer segments, and divides the sedimentary units within the layer segments using different methods according to the sedimentary environment of the layer segments to obtain a preliminary formation division result, and then corrects the preliminary division result to obtain an optimal formation boundary division result. This method first divides the development wells into multiple layer segments, and then divides the sedimentary units within the layer segments; each layer segment is treated as a unit, and a small layer segment is divided, improving the accuracy of formation division.

[0050] Other features and advantages of the technical solution of the present invention will be described in detail in the following specific implementation section. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] The drawings described herein are used to provide a further understanding of the present invention and constitute a part of the present invention. The schematic embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute an improper limitation to the present invention. In the drawings:

[0052] Figure 1 is a flowchart of a formation correlation method for a heterogeneous multi-layer sandstone oilfield provided by an embodiment of the present invention;

[0053] Figure 2 is a structural block diagram of a multi-level standard layer division model provided by an embodiment of the present invention;

[0054] Figure 3 is a diagram of the standard layer boundary division result provided by an embodiment of the present invention;

[0055] Figure 4 is a structural block diagram of a formation correlation system for a heterogeneous multi-layer sandstone oilfield provided by an embodiment of the present invention;

[0056] Figure 5It is a schematic diagram of the superimposed mode of single sand bodies in different deposition environments provided by the embodiments of the present invention. Detailed implementation manners

[0057] In order to make the technical solutions and advantages in the embodiments of the present invention clearer and more understandable, the following further describes the exemplary embodiments of the present invention in detail with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than an exhaustive list of all embodiments. It should be noted that, without conflict, the embodiments in the present invention and the features in the embodiments can be combined with each other.

[0058] In the description of the present invention, it should be understood that the orientation or positional relationship indicated by the terms "center", "longitudinal", "transverse", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", etc. is based on the orientation or positional relationship shown in the drawings, and is only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be construed as a limitation to the present invention.

[0059] In addition, the terms "first" and "second" are only used for descriptive purposes, and cannot be understood as indicating or implying relative importance or implicitly specifying the quantity of the indicated technical features. Thus, the features defined with "first" and "second" may explicitly or implicitly include one or more of such features. In the description of the present invention, the meaning of "plurality" is at least two, such as two, three, etc., unless otherwise specifically defined.

[0060] In the present invention, unless otherwise clearly defined and limited, terms such as "installation", "connection", "connection", "fixation", etc. should be understood in a broad sense. For example, it may be a fixed connection, a detachable connection, or integrated; it may be a mechanical connection, an electrical connection or communication with each other; it may be directly connected, or indirectly connected through an intermediate medium, and may be the internal communication of two elements or the interaction relationship between two elements. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific circumstances.

[0061] In the process of implementing the present invention, the inventors found that with the continuous deepening of oilfield development and the continuous increase in the number of development wells, it is difficult to meet the requirements of efficient oilfield development by relying solely on manual formation correlation. Especially for heterogeneous multi-layered sandstone oilfields, with nearly a hundred vertical layers, it is difficult to accurately identify and correlate the entire formation system. The traditional formation correlation method is based on high-resolution sequence stratigraphy. By using well logging curve characteristics that can reflect the preservation degree and accumulation rate of sediments, such as facies sequences, formation boundaries, and facies differentiation, the symmetry degree of formation cycles, the patterns of cycle thickening or thinning, and the amplitude and direction of facies dislocation at the positions crossing genetic formation boundaries are identified, and then formation boundaries are divided accordingly. However, the traditional formation correlation method is relatively simple and mainly deals with homogeneous thick sections of single wells. For heterogeneous multi-well sandstone oilfields, the accuracy of correlation and identification is insufficient.

[0062] In view of the above problems, an embodiment of the present invention provides a formation correlation method for heterogeneous multi-layered sandstone oilfields, including: dividing standard formation boundaries of different levels for each development well in the target area, and a section is formed between two adjacent standard formation boundaries of the development well; determining the sedimentary environment of each section according to the preset evaluation criteria for formation development characteristics in different sedimentary environments; dividing the boundaries of continuous sedimentary units of the section according to the sedimentary environment of each section to obtain a preliminary formation boundary division result; establishing a three-dimensional space model within the section based on the preliminary formation boundary division result, and correcting the formation boundary division result through the three-dimensional space model within the section to obtain an optimal formation boundary division result. The formation correlation method for heterogeneous multi-layered sandstone oilfields provided by the present invention divides the development wells into multiple sections, then determines the sedimentary environment of the sections, uses different methods to divide the sedimentary units within the sections according to the sedimentary environment of the sections to obtain a preliminary formation division result, and then corrects the preliminary division result to obtain an optimal formation boundary division result. This method first divides the development wells into multiple sections, and then divides the sedimentary units within the sections; each section is treated as a unit, and a small section is divided, which improves the accuracy of formation division.

[0063] Figure 1 It is a flowchart of a formation correlation method for heterogeneous multi-layered sandstone oilfields provided by an embodiment of the present invention. As Figure 1 shown, a formation correlation method for heterogeneous multi-layered sandstone oilfields provided by this embodiment includes the following steps:

[0064] S1. Divide standard formation boundaries of different levels for each development well in the target area;

[0065] S2. Determine the sedimentary environment of the section between two adjacent standard formation boundaries of each development well according to the preset evaluation criteria for formation development characteristics in different sedimentary environments;

[0066] S3. Divide the boundaries of continuous sedimentary units within the interval according to the sedimentary environment of the interval between two adjacent standard layer boundaries of each development well, and obtain a preliminary formation boundary division result;

[0067] S4. Based on the preliminary formation boundary division result, establish a three-dimensional spatial model within the interval, and correct the formation boundary division result through the three-dimensional spatial model within the interval to obtain an optimal formation boundary division result.

[0068] In step S1, the division of different-level standard layer boundaries for each development well in the target area, where a section is formed between two adjacent standard layer boundaries of a development well, includes:

[0069] S11. According to the core data of the target area, construct multiple different-level standard layer systems;

[0070] Specifically, the different-level standard layer systems include a first-level standard layer and a second-level standard layer. The stability of the first-level standard layer is more than 90%, and the stability of the second-level standard layer is 50 - 90%.

[0071] S12. Use an image segmentation method based on a convolutional neural network to construct a multi-level standard layer division model.

[0072] S13. Extract the standard layer characteristic logging curves of each level of the standard layer system, and use the standard layer logging curves of each level of the standard layer system as a training set to input into the multi-level standard layer division model to be trained for training; the multi-level standard layer division model is based on the image segmentation method of CNN, learns the standard layer characteristic logging curves, and obtains a trained multi-level standard layer division model.

[0073] S14. Obtain the original logging curves of each development well and input them into the trained multi-level standard layer division model to obtain different-level standard layer boundaries of each development well.

[0074] In step S1, the method further includes:

[0075] S15. Correct the different-level standard layer boundaries of each development well.

[0076] Specifically, according to the different-level standard layer boundaries of each development well, establish a three-dimensional mapping of the standard layer of the development well with the section as a unit, and correct the different-level standard layer boundaries of each development well through the three-dimensional mapping of the standard layer.

[0077] Figure 2 is the structural block diagram of the multi-level standard layer division model provided by the embodiment of the present invention; Figure 3 is the standard layer boundary division result diagram provided by the embodiment of the present invention;

[0078] Taking the western part of the central area with the richest heterogeneous multi-layer characteristics in the Daqing Oilfield as the research object, by establishing standard layer systems at different levels and using the image segmentation method based on CNN, a multi-level standard layer division model as shown in Figure 2 can be constructed. The standard layer boundary division as shown in Figure 3 can be carried out. The coincidence rate of the standard layer stratigraphic division is 96.3%. Using Petrel software, the three-dimensional mapping of the standard layer at multiple levels is realized, and the local distortion is corrected. The closing rate of the standard layer stratification boundary reaches 100%.

[0079] In step S2, determining the sedimentary environment of the interval between the adjacent two standard layer boundaries of each development well according to the preset evaluation criteria for the stratigraphic development characteristics in different sedimentary environments includes:

[0080] S21. Selecting the evaluation standard parameters for the stratigraphic development characteristics used to characterize the sedimentary environment, and constructing the evaluation criteria for the stratigraphic development characteristics in different sedimentary environments according to the evaluation standard parameters for the stratigraphic development characteristics;

[0081] S22. Obtaining the well logging interpretation data within each interval of each well in the target area; taking the interval 1 between standard layer 1 and standard layer 2 as an example, taking the depth of the adjacent standard layer boundary as the boundary, extracting the well logging interpretation data within the interval between the adjacent standard layers of all wells in the target area;

[0082] S23. Calculating the parameter values of the evaluation parameters for the stratigraphic development characteristics according to the well logging interpretation data within each interval;

[0083] S24. Determining the sedimentary environment of each interval from the evaluation criteria for the stratigraphic development characteristics in different sedimentary environments based on the calculated parameter values of the evaluation parameters for the stratigraphic development characteristics.

[0084] In step S21, the evaluation parameters for the stratigraphic development characteristics include at least three of the following: average effective thickness, average effective permeability, sand-to-ground ratio, drilling penetration rate, layering coefficient;

[0085] Among them, the calculation formula for the average effective thickness is:

[0086] The calculation formula for the average effective permeability is:

[0087] The calculation formula for the sand-to-ground ratio is:

[0088] The calculation formula for the drilling penetration rate is:

[0089] The calculation formula for the layering coefficient is:

[0090] wherein, i represents the i-th interval of the development well, j represents the j-th sedimentary unit within the i-th interval, n represents the number of sedimentary units, m represents the number of wells encountering intervals in the target area, h ei is the average effective thickness of the interval, K ei is the average effective permeability of the interval, D i is the encounter rate, Q i is the layering coefficient, h ej is the effective thickness of the j-th sedimentary unit within the interval, h j is the formation thickness of the j-th sedimentary unit within the interval, h syj is the sandstone thickness of the j-th sedimentary unit within the interval, k j is the effective permeability of the j-th sedimentary unit within the interval, N i is the number of wells encountering the interval, T is the total number of wells in the block, E i is the number of effective sandstone layers within the interval.

[0091] Based on the standard layer division results obtained from S1, using the method described in S2, the evaluation criteria for the development characteristics of strata in different sedimentary environments as shown in Table 1 were established. Through the calculation of relevant parameters of 779 wells, the environmental discrimination at the sandstone group level within 20 intervals was completed.

[0092] Before developing the development wells, through the exploration of the target area, the number of sedimentary units within the interval can be known, but the boundaries of the sedimentary units within the interval are not clear.

[0093] In step S3, the boundaries of the continuous sedimentary units within the interval are divided according to the sedimentary environment of the interval between the boundaries of two adjacent standard layers of each development well, and the preliminary formation boundary division result is obtained, including:

[0094] S31. Select multiple types of logging curves for formation division and correlation features, and synthesize the multiple types of logging curves for formation division and correlation features into a fusion feature curve; specifically, preferably the spontaneous potential SP, natural gamma ray GR, and microelectrode (micro-gradient RMG, micro-potential RMN) are used as the logging curves for formation division and correlation features. Using the feature curve filtering and synthesis method, all the logging curves in the target area are comprehensively processed to maximize the curve signal-to-noise energy ratio, and the fusion feature curves of the reference well and the target well in the target area are obtained.

[0095] S32. According to the sedimentary environment of each interval, apply the fusion feature curve to divide the boundaries of the continuous sedimentary units of the current interval to obtain the preliminary formation boundary division result.

[0096] In step S32, the multiple sedimentary units of the current interval are divided according to the sedimentary environment of each interval by applying the fusion feature curve to obtain the preliminary formation division result, including:

[0097] S321. If the sedimentary environment of the layer section is the frontal facies, then according to the fused characteristic curve, the wavelet transform method and the local INPEFA logging cycle analysis method are used to divide the boundaries of the continuous sedimentary units within the current layer section to obtain a preliminary formation division result;

[0098] S322. If the sedimentary environment of the layer section is the plain facies, then according to the fused characteristic curve, the sand body superposition pattern constraint method and the continuous formation correlation and division method of dynamic programming are used to divide the boundaries of the continuous sedimentary units within the current layer section to obtain a preliminary formation division result.

[0099] Specifically, if it is determined that the sedimentary environment of the layer section is the frontal facies, the wavelet transform and the local INPEFA logging cycle analysis method are used to divide the internal continuous sedimentary units;

[0100] If it is determined that the sedimentary environment of the layer section is the plain facies, the continuous formation correlation and division method based on the sand body superposition pattern constraint and dynamic programming is used to divide the internal continuous sedimentary units.

[0101] Figure 5 It is a schematic diagram of the single sand body superposition pattern in different sedimentary environments provided by the embodiment of the present invention. As Figure 5 shown, if the layering coefficient Qi < 1.0, refer to the Figure 5 sand body superposition pattern of the flood plain facies in; if the layering coefficient Qi > 1.0, refer to the Figure 5 sand body superposition pattern of the distributary plain facies in. Based on the results obtained from S1 and S2, through the logging curve processing of 779 wells, through the learning of the formation pattern shown in Figure 4 the continuous sedimentary units between the standard layer sections are divided, and the formation boundaries of 78 sedimentary units (sub-layers) in 779 wells are completed, forming a set of data for formation boundary division.

[0102] In this embodiment, the continuous sedimentary unit means that multiple sedimentary units are adjacent in sequence.

[0103] In step S4, the three-dimensional space model within the layer section is established based on the preliminary formation boundary division result, and the formation boundary division result is corrected through the three-dimensional space model within the layer section to obtain the optimal formation boundary division result, including:

[0104] S41. According to the preliminary formation boundary division result, a three-dimensional space model within the layer section is established to represent the continuous formation distribution within the layer section;

[0105] S42. The local cross-stratum abnormal points are corrected on the three-dimensional space model within the layer section to obtain a preliminarily corrected formation division result;

[0106] S43. Establish an isopach map of the planar formation thickness within each interval, obtain and mark the local abnormal points of the preliminary formation division result after initial refinement in the isopach map of the planar formation thickness within the interval;

[0107] S44. Based on the sedimentary environment of the sedimentary unit, and based on the three-dimensional spatial model within the interval and the isopach map of the planar thickness within the interval, correct the local abnormal points to obtain the optimal formation boundary division result.

[0108] Specifically, based on the formation correlation result obtained in S3, i.e., the preliminary formation boundary division result. Use modeling software such as Petrel or RMS to establish a horizon model, realize the three-dimensional spatial model within the interval of the continuous formation distribution between adjacent standard horizons, correct the local cross-strata abnormal points spatially, and export the formation division result;

[0109] S42: Establish an isopach map of the planar formation thickness of a single sub-layer, and mark the local abnormal parts;

[0110] S43: Through the linkage of the isopach map of the planar formation thickness of a single sub-layer in the two-dimensional plane and the three-dimensional spatial model within the interval, correct the local abnormal parts to complete the fine adjustment of the formation boundary.

[0111] S431: For the sub-layers in the frontal facies, the isopach characteristics are of a gradual change type. For the mutant distortion parts, establish cross-well and longitudinal-well profiles, and use adjacent well constraints to correct the mutant distortion;

[0112] S432: For the sub-layers in the plain facies, the isopachs show the characteristics of alternating high and low values and different widths. Considering that the formation thickness between adjacent standard horizons is similar, within the adjacent standard horizons, use the method of continuous formation image recognition constrained by the sedimentary pattern of fluvial sand bodies to circle out the areas with unfavorable constraints of each single sand layer pattern between the standard horizons, establish cross-well and longitudinal-well profiles, correct the areas with unfavorable pattern constraints, and complete the rapid correlation of the formation at the sedimentary unit level (i.e., single sand layer).

[0113] Based on the result obtained in S3, conduct logging curve processing for 779 wells, and based on the method described in S4 and the formation pattern learning as Figure 4 shown, divide the continuous sedimentary units between the standard intervals, complete the formation boundary division of 78 sedimentary units (sub-layers) in 779 wells, and the error rate of the well-zoned comparison of the sedimentary unit wells is 0.97%.

[0114] Figure 4 It is the structural block diagram of a formation correlation system for a heterogeneous multi-layer sandstone oilfield provided by an embodiment of the present invention. As Figure 4As shown in the figure, a stratigraphic correlation system for heterogeneous multi-layer sandstone oilfields provided by this embodiment includes: a layer section division unit for dividing standard layer boundaries of different levels for each development well in the target area; a layer section sedimentary environment confirmation unit for determining the sedimentary environment of the layer section between two adjacent standard layer boundaries of each development well according to the preset evaluation criteria for the development characteristics of strata in different sedimentary environments; a stratigraphic boundary division unit for dividing the boundaries of continuous sedimentary units within the layer section according to the sedimentary environment of the layer section between two adjacent standard layer boundaries of each development well to obtain a preliminary stratigraphic boundary division result; and a correction unit for establishing a three-dimensional space model within the layer section based on the preliminary stratigraphic boundary division result and correcting the stratigraphic boundary division result through the three-dimensional space model within the layer section to obtain an optimal stratigraphic boundary division result.

[0115] The above-mentioned stratigraphic correlation system for heterogeneous multi-layer sandstone oilfields is implemented by using the above-mentioned stratigraphic correlation method for heterogeneous multi-layer sandstone oilfields. The method includes:

[0116] S1. Divide standard layer boundaries of different levels for each development well in the target area;

[0117] S2. Determine the sedimentary environment of the layer section between two adjacent standard layer boundaries of each development well according to the preset evaluation criteria for the development characteristics of strata in different sedimentary environments;

[0118] S3. Divide the boundaries of continuous sedimentary units within the layer section according to the sedimentary environment of the layer section between two adjacent standard layer boundaries of each development well to obtain a preliminary stratigraphic boundary division result;

[0119] S4. Establish a three-dimensional space model within the layer section based on the preliminary stratigraphic boundary division result and correct the stratigraphic boundary division result through the three-dimensional space model within the layer section to obtain an optimal stratigraphic boundary division result.

[0120] In step S1, the division of standard layer boundaries of different levels for each development well in the target area, and a layer section is formed between two adjacent standard layer boundaries of the development well, including:

[0121] S11. Construct multiple standard layer systems of different levels according to the core data of the target area;

[0122] Specifically, the standard layer systems of different levels include a first-level standard layer and a second-level standard layer. The stability of the first-level standard layer is more than 90%, and the stability of the second-level standard layer is 50 - 90%.

[0123] S12. Use an image segmentation method based on a convolutional neural network to construct a multi-level standard layer division model.

[0124] S13. Extract the standard layer characteristic logging curves of the standard layer systems at all levels, and use the standard layer logging curves of the standard layer systems at all levels as the training set to input into the multi-level standard layer division model to be trained for training; the multi-level standard layer division model is based on the image segmentation method of CNN, learns the standard layer characteristic logging curves, and obtains the trained multi-level standard layer division model.

[0125] S14. Obtain the original logging curves of each development well and input them into the trained multi-level standard layer division model to obtain the standard layer boundaries at different levels of each development well.

[0126] In step S1, the method further includes:

[0127] S15. Correct the standard layer boundaries at different levels of each development well.

[0128] Specifically, according to the standard layer boundaries at different levels of each development well, establish a three-dimensional map of the standard layers of the development well with the layer segment as the unit, and correct the standard layer boundaries at different levels of each development well through the three-dimensional map of the standard layers.

[0129] Taking the western part of the central area with the richest heterogeneous multi-layer characteristics in the Daqing Oilfield as the research object, by establishing standard layer systems at different levels and using the image segmentation method based on CNN, construct a Figure 2 multi-level standard layer division model as shown, and can carry out the standard layer boundary division as shown Figure 3 The coincidence rate of the standard layer stratigraphic division is 96.3%. Use Petrel software to realize the three-dimensional mapping of the standard layers at multiple levels, correct the local distortion, and the closure rate of the standard layer stratification boundary reaches 100%.

[0130] In step S2, the determination of the sedimentary environment of the layer segment between the adjacent two standard layer boundaries of each development well according to the preset evaluation criteria for the stratigraphic development characteristics in different sedimentary environments includes:

[0131] S21. Select the evaluation standard parameters for the stratigraphic development characteristics used to characterize the sedimentary environment, and construct the evaluation criteria for the stratigraphic development characteristics in different sedimentary environments according to the evaluation standard parameters for the stratigraphic development characteristics; specifically, according to the historical data, set the thresholds for each evaluation standard parameter for the stratigraphic development characteristics in different sedimentary environment cases to form the evaluation criteria for the stratigraphic development characteristics in different sedimentary environments.

[0132] S22. Obtain the logging interpretation data of each well in each layer segment within the target area; taking the layer segment 1 between standard layer 1 and standard layer 2 as an example, extract the logging interpretation data of the layer segment between the adjacent standard layers of all wells in the target area with the adjacent standard layer boundary depth as the boundary;

[0133] S23. Calculate the parameter values of the formation development characteristic evaluation parameters according to the logging interpretation data within each interval;

[0134] S24. Determine the sedimentary environment of each interval from the formation development characteristic evaluation criteria for different sedimentary environments based on the calculated parameter values of the formation development characteristic evaluation parameters.

[0135] In step S21, the formation development characteristic evaluation parameters include at least three of the following: average effective thickness, average effective permeability, sand-to-shale ratio, encounter rate, layering coefficient;

[0136] Among them, the calculation formula for the average effective thickness is:

[0137] The calculation formula for the average effective permeability is:

[0138] The calculation formula for the sand-to-shale ratio is:

[0139] The calculation formula for the encounter rate is:

[0140] The calculation formula for the layering coefficient is:

[0141] Among them, i is the i-th interval of the development well, j is the j-th sedimentary unit within the i-th interval, n is the number of sedimentary units, m is the number of wells encountering the intervals in the target area, h ei is the average effective thickness of the interval, K ei is the average effective permeability of the interval, D i is the encounter rate, Q i is the layering coefficient, h ej is the effective thickness of the j-th sedimentary unit within the interval, h j is the formation thickness of the j-th sedimentary unit within the interval, h syj is the sandstone thickness of the j-th sedimentary unit within the interval, k j is the effective permeability of the j-th sedimentary unit within the interval, N i is the number of wells encountering the intervals, T is the total number of wells in the block, E i is the number of effective sandstone layers within the interval. This embodiment also provides the formation development characteristic evaluation criteria for different sedimentary environments as shown in Table 1 below: Formation development characteristic evaluation criteria for different sedimentary environments

[0142]

[0143]

[0144] Based on the standard layer division results obtained from S1, using the method described in S2, an evaluation standard for the development characteristics of strata in different sedimentary environments as shown in Table 1 was established. Through the calculation of relevant parameters of 779 wells, the environmental discrimination at the sandstone group level within 20 intervals was completed.

[0145] Before the development of development wells, through the exploration of the target area, the number of sedimentary units within the interval can be known, but the boundaries of the sedimentary units within the interval are not clear.

[0146] In step S3, the boundaries of the continuous sedimentary units of the interval are divided according to the sedimentary environment of the interval between the boundaries of two adjacent standard layers of each development well, and a preliminary formation boundary division result is obtained, including:

[0147] S31. Select multiple types of logging curves for formation division and correlation characteristics, and synthesize the multiple types of logging curves for formation division and correlation characteristics into a fused characteristic curve; specifically, preferably select spontaneous potential SP, natural gamma ray GR, and microelectrodes (micro-gradient RMG, micro-potential RMN) as the logging curves for formation division and correlation characteristics, and use the characteristic curve filtering synthesis method to comprehensively process all the logging curves in the target area to maximize the signal-to-noise energy ratio of the curves, and obtain the fused characteristic curves of the reference well and the target well in the target area.

[0148] S32. According to the sedimentary environment of each interval, use the fused characteristic curve to divide the boundaries of the continuous sedimentary units within the current interval to obtain a preliminary formation boundary division result.

[0149] In step S32, the boundaries of the continuous sedimentary units within the current interval are divided according to the sedimentary environment of each interval by using the fused characteristic curve to obtain a preliminary formation division result, including:

[0150] S321. If the sedimentary environment of the interval is the frontal facies, then according to the fused characteristic curve, use the wavelet transform method and the local INPEFA logging cycle analysis method to divide the multiple sedimentary units of the current interval to obtain a preliminary formation division result;

[0151] S322. If the sedimentary environment of the interval is the plain facies, then according to the fused characteristic curve, use the sand body stacking pattern constraint method and the continuous formation correlation and division method of dynamic programming to divide the multiple sedimentary units of the current interval to obtain a preliminary formation division result.

[0152] Specifically, if it is determined that the sedimentary environment of the interval is the frontal facies, use the wavelet transform and the local INPEFA logging cycle analysis method to divide the internal continuous sedimentary units;

[0153] If it is determined that the sedimentary environment of the interval is the plain facies, use the continuous formation correlation and division method based on sand body stacking pattern constraint and dynamic programming to divide the internal continuous sedimentary units.

[0154] Figure 5 It is a schematic diagram of the superposition pattern of single sand bodies in different deposition environments provided by the embodiments of the present invention. As Figure 5 shown, if the layering coefficient Qi < 1.0, refer to the Figure 5 superposition pattern of sand bodies in the floodplain facies in Figure 5 ; if the layering coefficient Qi > 1.0, refer to the

[0155] superposition pattern of sand bodies in the distributary plain facies in Figure 4 . Based on the results obtained from S1 and S2, by processing the logging curves of 779 wells, through learning the

[0156] stratigraphic pattern shown, continuous sedimentary units between standard stratigraphic sections are divided, and the stratigraphic boundaries of 78 sedimentary units (sub-layers) of 779 wells are completed, forming a set of data for stratigraphic boundary division.

[0157] In step S4, the three-dimensional spatial model within the layer segment is established based on the preliminary stratigraphic boundary division result, and the stratigraphic boundary division result is corrected through the three-dimensional spatial model within the layer segment to obtain the optimal stratigraphic boundary division result, including:

[0158] S41. According to the preliminary stratigraphic boundary division result, establish a three-dimensional spatial model within the layer segment for characterizing the continuous stratigraphic distribution within the layer segment;

[0159] S42. Correct local cross-stratum abnormal points on the three-dimensional spatial model within the layer segment to obtain a preliminarily corrected stratigraphic division result;

[0160] S43. Establish an isopach map of the planar stratigraphic thickness within each layer segment, and obtain and mark the local abnormal points of the preliminarily corrected stratigraphic division result in the isopach map of the planar stratigraphic thickness within the layer segment;

[0161] S44. According to the deposition environment of the sedimentary unit, based on the three-dimensional spatial model within the layer segment and the isopach map of the planar thickness within the layer segment, correct the local abnormal points to obtain the optimal stratigraphic boundary division result.

[0162] Specifically, based on the stratigraphic correlation result obtained from S3, that is, the preliminary stratigraphic boundary division result. Use modeling software such as Petrel or RMS to establish a horizon model to realize the three-dimensional spatial model within the layer segment of the continuous stratigraphic distribution between adjacent standard horizons, correct local cross-stratum abnormal points in space, and export the stratigraphic division result;

[0163] S42: Establish an isopach map of the planar stratigraphic thickness of a single sub-layer and mark the local abnormal parts;

[0164] S43: Through the linkage between the isopach map of a single small layer plane in a two-dimensional plane and the three-dimensional space model within the layer section, local abnormal areas are corrected to complete the fine adjustment of the stratigraphic boundary.

[0165] S431: For the small layers in the frontal facies, the contour feature is a gradual change type. For the mutated and distorted areas with abrupt change type, transverse and longitudinal cross-well profiles are established, and the mutation and distortion are corrected using the constraints of adjacent wells.

[0166] S432: For the small layers in the plain facies, the contour lines show the characteristics of alternating high and low values and varying widths. Considering that the stratigraphic thickness between adjacent standard layers is similar, within the adjacent standard layers, using the continuous stratigraphic image recognition method constrained by the sedimentary pattern of fluvial sand bodies, the unfavorable areas constrained by each single sand layer pattern between the standard layers are circled, transverse and longitudinal cross-well profiles are established, the unfavorable areas constrained by the pattern are corrected, and the rapid correlation of the sedimentary unit-level strata (i.e., single sand layers) is completed.

[0167] Based on the results obtained from S3, well logging curve processing of 779 wells is carried out. Based on the method described in S4 and the stratigraphic pattern learning as Figure 4 shown, the continuous sedimentary units between the standard layer sections are divided, and the stratigraphic boundaries of 78 sedimentary units (small layers) in 779 wells are completed. The error rate of the well-layer comparison quality of the sedimentary unit around the well is 0.97%.

[0168] The embodiment of the present invention also provides a computer device, including: a memory, a processor, and a computer program. The computer program is stored in the memory and is configured to be executed by the processor to implement the above-mentioned stratigraphic correlation method for heterogeneous multi-layer sandstone oilfields.

[0169] The embodiment of the present invention also provides a machine-readable storage medium, on which computer program instructions are stored. When the computer program instructions are executed by the processor, the above-mentioned stratigraphic correlation method for heterogeneous multi-layer sandstone oilfields is implemented.

[0170] Based on the characteristic curves of standard layers at different levels, this application divides the boundaries of standard layers of all wells in the target area; selects the evaluation parameters of formation development characteristics in the interval between adjacent standard layers to quickly determine the sedimentary environment of the interval; for the characteristics of stratigraphic superposition patterns in different sedimentary environments, different stratigraphic correlation methods are adopted to vertically divide the continuous stratigraphic units between adjacent standard layers; through the interactive comparison of the three-dimensional spatial representation of continuous stratigraphic boundaries and the isopach map of single sand layer formation thickness, the abnormal information of local single-well stratigraphic boundaries between adjacent standard layers is corrected, and the rapid stratigraphic correlation at the sedimentary unit level is completed. This application realizes the rapid identification and correlation of multi-level standard layers and constructs a rapid judgment method for different sedimentary environments around the characteristics of different formation development patterns and different stratigraphic division and correlation methods in multi-layer heterogeneous sandstone oilfields, thus quickly completing the rapid stratigraphic correlation from sandstone group level to small layer level in different sedimentary environments of heterogeneous multi-layer sandstone oilfields, greatly improving the accuracy and efficiency of stratigraphic correlation work, and can be widely applied in the field of basic geological research with high timeliness requirements.

[0171] Those skilled in the art should understand that the embodiments of the present invention can be provided as a method, a system, or a computer program product. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present invention can 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. The solutions in the embodiments of the present invention can be implemented in various computer languages, for example, object-oriented programming languages such as Java and interpreted scripting languages such as JavaScript.

[0172] The present invention is described with reference to the flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to the 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 the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate a device for realizing the functions specified in Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks.

[0173] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, so that the instructions stored in the computer-readable memory generate a manufactured article including an instruction device, and the instruction device realizes the functions in the processFigure 1 one process or multiple processes and / or boxes Figure 1 the functions specified in one box or multiple boxes.

[0174] These computer program instructions can also be loaded onto a computer or other programmable data processing device, so that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process, and thus the instructions executed on the computer or other programmable device provide for implementing the process Figure 1 one process or multiple processes and / or boxes Figure 1 the steps of the functions specified in one box or multiple boxes.

[0175] Although the preferred embodiments of the present invention have been described, those skilled in the art can make additional changes and modifications once they learn the basic creative concept. Therefore, the appended claims are intended to be construed to include the preferred embodiments as well as all changes and modifications falling within the scope of the present invention.

[0176] Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalent technologies, the present invention also intends to include these modifications and variations.

Claims

1. A formation correlation method for heterogeneous multi-layer sandstone oilfields, characterized in that, Including: Dividing the standard layer boundaries of different levels for each development well in the target area; Determining the sedimentary environment of the interval between two adjacent standard layer boundaries of each development well according to the preset evaluation criteria for the stratigraphic development characteristics of different sedimentary environments; Dividing the boundaries of continuous sedimentary units within the interval according to the sedimentary environment of the interval between the two adjacent standard layer boundaries to obtain a preliminary result of stratigraphic boundary division; Establishing a three-dimensional spatial model within the interval based on the preliminary result of stratigraphic boundary division, and correcting the stratigraphic boundary division result through the three-dimensional spatial model within the interval to obtain an optimal stratigraphic boundary division result.

2. The formation correlation method for heterogeneous multi-layer sandstone oilfields according to claim 1, characterized in that, The dividing the standard layer boundaries of different levels for each development well in the target area includes: Constructing multiple standard layer systems of different levels according to the core data of the target area; Constructing a multi-level standard layer division model using an image segmentation method based on a convolutional neural network; Extracting the standard layer characteristic logging curves of each level of the standard layer system, and using the standard layer logging curves of each level of the standard layer system as a training set to input into the multi-level standard layer division model for training; Obtaining the original logging curves of each development well and inputting them into the trained multi-level standard layer division model to obtain the standard layer boundaries of different levels for each development well.

3. The formation correlation method for heterogeneous multi-layer sandstone oilfields according to claim 2, characterized in that, The method further includes: correcting the standard layer boundaries of different levels for each development well.

4. The formation correlation method for heterogeneous multi-layer sandstone oilfields according to claim 1, characterized in that, The determining the sedimentary environment of the interval between two adjacent standard layer boundaries of each development well according to the preset evaluation criteria for the stratigraphic development characteristics of different sedimentary environments includes: Selecting the evaluation standard parameters for the stratigraphic development characteristics representing the sedimentary environment, and constructing the evaluation criteria for the stratigraphic development characteristics of different sedimentary environments according to the evaluation standard parameters for the stratigraphic development characteristics; Obtaining the logging interpretation data within each interval of each well in the target area; Calculating the parameter values of the evaluation parameters for the stratigraphic development characteristics according to the logging interpretation data within each interval; Determining the sedimentary environment of each interval from the evaluation criteria for the stratigraphic development characteristics of different sedimentary environments based on the calculated parameter values of the evaluation parameters for the stratigraphic development characteristics.

5. The formation correlation method for heterogeneous multi-layer sandstone oilfields according to claim 4, characterized in that, The evaluation parameters for the stratigraphic development characteristics include: average effective thickness, average effective permeability, sand-to-shale ratio, drilling encounter rate, layering coefficient; Among them, the calculation formula for the average effective thickness is as follows: The calculation formula for the average effective permeability is as follows: The calculation formula for the sand ratio is as follows: The calculation formula for the drilling encounter rate is as follows: The calculation formula for the layering coefficient is as follows: where, i is the i-th interval of the development well, j is the j-th sedimentary unit within the i-th interval, n is the number of sedimentary units, m is the number of wells drilled through the intervals in the target area, h ei is the average effective thickness of the interval, K ei is the average effective permeability of the interval, D i is the encounter rate, Q i is the layering coefficient, h ej is the effective thickness of the j-th sedimentary unit within the interval, h j is the formation thickness of the j-th sedimentary unit within the interval, h syj is the sandstone thickness of the j-th sedimentary unit within the interval, k j is the effective permeability of the j-th sedimentary unit within the interval, N i is the number of wells drilled through the intervals, T is the total number of wells in the block, E i is the number of effective sandstone layers within the interval.

6. The formation correlation method for heterogeneous multi-layer sandstone oilfields according to claim 1, characterized in that, The dividing the boundaries of continuous sedimentary units within the interval according to the sedimentary environment of the interval between the two adjacent standard layer boundaries to obtain a preliminary result of stratigraphic boundary division includes: Selecting multiple types of stratigraphic division and correlation characteristic logging curves, and synthesizing the multiple types of stratigraphic division and correlation characteristic logging curves into a fusion characteristic curve; According to the sedimentary environment of each interval, applying the fusion characteristic curve to divide the boundaries of continuous sedimentary units within the current interval to obtain a preliminary result of stratigraphic boundary division.

7. The formation correlation method for heterogeneous multi-layer sandstone oilfields according to claim 6, characterized in that, The dividing the multiple sedimentary units of the current interval according to the sedimentary environment of each interval by applying the fusion characteristic curve to obtain a preliminary result of stratigraphic division includes: If the sedimentary environment of the interval is the frontal facies, then according to the fusion characteristic curve, using the wavelet transform method and the local INPEFA logging cycle analysis method to divide the boundaries of continuous sedimentary units within the current interval to obtain a preliminary result of stratigraphic division; If the sedimentary environment of the layer section is a plain facies, then according to the fusion characteristic curve, the sand body stacking mode constraint method and the continuous stratigraphic correlation division method of dynamic programming are adopted to divide the boundaries of the continuous sedimentary units within the current layer section, and a preliminary stratigraphic division result is obtained.

8. The formation correlation method for heterogeneous multi-layer sandstone oilfields according to claim 1, characterized in that, Based on the preliminary stratigraphic boundary division result, a three-dimensional spatial model within the layer section is established, and the stratigraphic boundary division result is corrected through the three-dimensional spatial model within the layer section to obtain an optimal stratigraphic boundary division result, including: According to the preliminary stratigraphic boundary division result, a three-dimensional spatial model within the layer section for characterizing the continuous stratigraphic distribution within the layer section is established; Local cross-stratum abnormal points are corrected on the three-dimensional spatial model within the layer section to obtain a preliminarily corrected stratigraphic division result; A planar stratigraphic thickness isoline map of each sedimentary unit within the layer section is established, and the local abnormal points of the preliminarily corrected stratigraphic division result in the planar stratigraphic thickness isoline map within the layer section are obtained and marked; According to the sedimentary environment of the sedimentary unit, based on the three-dimensional spatial model within the layer section and the planar thickness isoline map of each sedimentary unit within the layer section, the local abnormal points are corrected to obtain an optimal stratigraphic boundary division result.

9. A formation correlation system for heterogeneous multi-layer sandstone oilfields, characterized in that, Including: A layer section division unit for dividing standard layer boundaries of different levels for each development well within the target area; A layer section sedimentary environment confirmation unit for determining the sedimentary environment of the layer section between two adjacent standard layer boundaries of each development well according to the preset evaluation criteria for the stratigraphic development characteristics of different sedimentary environments; A stratigraphic boundary division unit for dividing the boundaries of the continuous sedimentary units within the layer section according to the sedimentary environment of the layer section between two adjacent standard layer boundaries of each development well to obtain a preliminary stratigraphic boundary division result; A correction unit for establishing a three-dimensional spatial model within the layer section based on the preliminary stratigraphic boundary division result and correcting the stratigraphic boundary division result through the three-dimensional spatial model within the layer section to obtain an optimal stratigraphic boundary division result.

10. A computer device, characterized in that, Including: A memory; A processor; And A computer program; Wherein, the computer program is stored in the memory and is configured to be executed by the processor to implement a stratigraphic correlation method for a heterogeneous multi-layer sandstone oilfield as described in any one of claims 1-8.

Images